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ABSTRACT 


The Normalized Radar Cross-Section (NRCS), the fundamental measurement made 
by radar scatterometers, was obtained as part of the Water-Air Vertical Exchanges 1987 
(WAVES87) experiment. The experiment was designed to evaluate the effects of 
environmental parameters on the NRCS and was performed from a research tower located 
in Lake Ontario, on which two microwave scatterometers operating at 14.0 and 5.0GHz 
were installed for six weeks in the autumn of 1987. The novel aspect of this experiment 
was that the 14.0GHz radar automatically rotated through 300° in azimuth angle at six 
different incidence angles to the water surface, accompanied by simultaneous 
measurements of wind stress and high resolution directional wave spectra. Therefore, the 
incidence and azimuthal angle behavior of the NRCS was examined as a function of wind 
speed, friction velocity, wind direction, wave direction and atmospheric stability. 

The dependence of the NRCS on wind speed for various incidence angles is similar to 
previous results. However, the slope exponents of the NRCS vs. 19.5m wind speed 
curves at intermediate incidence angles are higher than the corresponding open ocean 
measurements. Scaling the lake neutral wind speed data by the ratio of lake to ocean drag 
coefficients reduces the slopes of the curves and suggests the drag coefficient has a sea 
state dependence. The correlation between NRCS and neutral wind speed at 1m 1s higher 
(0.91) than between the NRCS and friction velocity (0.73 at 40°). The minima in the 
sinusoidal modulation of the NRCS as a function of relative wind angle (the angle between 
the wind and antenna directions) are often shifted (by as much as 45°) such that the minima 
do not always occur at cross-wind angles. Instead, the angular distance between the 
NRCS minima in the case of a wind-wave sea appears to approximate the directional spread 
of the waves about the upwind direction, generally rather less than 180°. The degree of 
sinusoidal modulation of the NRCS with relative wind angle is highly correlated with 
significant slope and inverse wave age at 20° incidence angle (0.90) and moderately 
' correlated at 40" (0.75); 1.e., increased azimuthal modulation at 20° is associated with a 
steeper wave field. The dependence of the NRCS on atmospheric stability shows the 
NRCS to decrease by about 5dB between air-water temperature differences of about -16 to 
+10°C. This stability effect is removed by parameterization of the NRCS in terms of either 
the friction velocity or neutral wind speed at 1m, with the neutral wind speed providing the 
best normalization of the data. The results show that radar scatterometers are an especially 
sensitive means by which to study the air-sea interface: the magnitudes of the 5.0GHz and 
14.0GHz NRCS respond nearly instantaneously to changes in the near-surface neutral 
wind speed, but the directionality of the (Ku-band) NRCS is the result of complicated 
interrelationships among the influencing environmental variables. 
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I. INTRODUCTION 


A. REMOTE SENSING OF THE OCEAN 

In the three decades since the launch of the first meteorological satellite, 
the enormous potential of remote sensing of the Earth's environment has 
been proved. Products of visible and infrared imagery are routinely used to 
monitor the evolution of weather systems, oceanic fronts and eddies. Passive 
infrared and microwave radiometers provide vertical temperature soundings 
of the atmosphere and information on its moisture content. Active 
microwave radars serve the dual purpose of measuring sea surface elevation 
and waves, and inferring sea surface winds over areal expanses, both with 
accuracies never before realized. To capitalize on these now proven 
techniques, space agencies have initiated ambitious programs to build and 
launch environment observing satellites in the early 1990's. 

At the heart of these programs are intriguing and unresolved 
questions, answers to which will provide the basis for next generation 
instrumentation and delimit the applicability of remotely sensed data. 
Questions raised are 1) theoretical: How are the spectral emissions/reflections 
from the earth related to geophysical variables of interest?, 2) technical: How 
can advances in computer and engineering sciences be incorporated into 
remote sensing instruments?, and 3) operational: How are the huge volumes 
of data to be received processed, distributed and analyzed in near real-time? 
The issue addressed in this study is of the first general type. In particular, the 


purpose of this study is to evaluate environmental effects on the functional 


form vvhich relates radar echoes from the sea surface to the sea surface vvind 


vector. 


B. THE RESEARCH PROBLEM 
1. Statement of the Problem 

Radar backscatter at intermediate incidence angles (20 to 70 deg) can 
be related by a "model function" to the sea surface wind speed and direction 
due to backscatter which is produced primarily by "Bragg resonant” reflections 
of the incident radiation from successive crests of short, wind generated 
capillary-gravity waves (0.2 to 20 cm wavelength). For the signal to be 
scattered to the radar, these waves must have components traveling along the 
direction of the antenna line-of-sight. Therefore, the backscattered power 
displays a strong angular dependence: maxima occur when the antenna is 
pointed either upwind or downwind, minima occur when the antenna is 
pointed crosswind. Furthermore, the upwind maximum is usually larger 
than the downwind maximum; this difference is enhanced for horizontally 
polarized radiation. The idea, then, behind scatterometer wind vector 
determination, is that a particular area on the ocean surface viewed from 
several different directions; ie., azimuth angles, yields different wind 
components per azimuth that can be used to determine a single wind vector. 

The above information has been incorporated into previous model 
functions by writing the normalized radar cross-section (NRCS) of the sea as a 
function of radar incidence angle, radar polarization, and azimuth angle 
relative to the true wind vector (Jones, et al., 1977; Moore and Fung, 1979; 
Schroeder, et al., 1982). However, Daley, et al., (1984) identify the deficiencies 


in this simple type of formulation and list the additional dependencies which 


may need to be included in the scatterometer model function. Based on 


recent research, these dependencies are: 


(1) cross-section dependence on wave slope and wave direction relative to 
the wind vector (Keller, et al., 1985, 1989; Li, et al., 1989), 


(2) cross-section dependence on air-sea temperature difference, and thus 
on Stability of the marine surface layer (e.g., Keller, et al., 1989), 


(3) a lack of self-consistency of the model between vertical and horizontal 
polarizations (Wentz, 1984; Woiceshyn, et al., 1986), 


(4) decreased sensitivity of the model function at low wind speeds (below 7 
m/s), (Woiceshyn, et al., 1986), 


(5) cross-section dependence on surface contaminants and films 
(Huhnerfuss, et al., 1983), and 


(6) cross-section dependence on viscosity of the water as a function of 
temperature (Woiceshyn, 1986; Donelan and Pierson, 1987). 


For the scatterometer wind velocity to meet the stringent operational 
requirements of the Navy (*2m/s, +2.5°) , an updated model function is 
required. 
2. Research Objective 

The effects of the various environmental parameters listed above are 
represented schematically in Figure 1. An increase in wind speed produces an 
associated increase in the NRCS through amplification of the capillary-gravity 
waves from which the electromagnetic waves are reflected. Observations 
show that the wind speed dependence can be represented as a power law 
(Figure 2a). Referring to Figure 2b, the cross-section varies with wind 
direction as a result of the anisotropy of the capillary-gravity waves with the 
relative maxima occurring when the radar beam is directed perpendicular to 
the wavecrest (upwind/downwind) and the minima occurring when the 
radar beam is parallel to the wavecrests (crosswind). This behavior is rather 


successfully modeled as a truncated Fourier cosine series. 
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Figure 1. Environmental Effects on the Normalized Radar Cross-Section 

If the short waves are tilted towards the radar and modulated in 
amplitude and frequency by the slope of the longer waves on which they ride, 
the NRCS will again be increased (Figure 2c). Observational evidence for the 
long wave influence has been reported by Keller, et al., (1985,1989) who show 
that the cross-sections increase with wave slope only when the atmosphere is 
stably stratified. Measurements during neutral and unstable conditions 
yielded no discernable dependence of the NRCS on long wave slope; it is 
suspected that the influence of stability is a consequence of the coupling 


between the surface waves and wind turbulence (see Figure 2d). 
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Figure 2. Functional Dependence of Normalized Radar Cross-Section on a) 
Wind Speed for Various Angles of Incidence, 0 (Schroeder, et al., 1982); b) 
Wind Direction (Jones, et al. 1977); c) Long Wave Slope (Keller, et al., 1985); 
d) Atmospheric Stratification (Keller, et al., 1989) 


Unstable atmospheric stratification, as indicated by a negative air-sea 
temperature difference (Figure 2d), results in 1) a slightly higher near-surface 
wind speed, 2) a wind profile which deviates from the neutral stratification, 
logarithmic form, and 3) an increase in the surface drag coefficient. The 
NRCS is seen to decrease in response to these changes in the surface layer as 
the stratification changes from unstable to stable. | 

The dissipation and growth rates of the capillary-gravity waves have 
been shown to be proportional to water viscosity (Lamb, 1932) which is an 
inverse function of water temperature; therefore, this mechanism has been 
proposed to explain observations which show an increase in the radar 
backscatter for increasing surface temperature, particularly at low wind speeds 
(Woiceshyn, et al., 1986; Donelan and Pierson, 1987; Kahma and Donelan, 
1988). Keller and Plant (1988) present recent wavetank measurements which 
do not show a pronounced temperature dependence. Additional 
observations are needed to reach a definite conclusion. 

The physical processes which influence radar backscatter all fall into 
the realm of describing, modeling and measuring the microphysics of the 
air/water surface layers. Among other institutions which study these 
processes, the Canada Centre for Inland Waters has supported air-water 
interaction experiments since 1976 under its WAVES (Water-Air Vertical 
Exchange Studies) program (Tsanis and Donelan, 1987). The primary 
platforms for the CCIW program are laboratory wavetanks and a tower at the 
western end of Lake Ontario, positioned 1.1km offshore from Hamilton, 
Ontario. In the autumn of three consecutive years (October-December, 1985- 


1987), the tower was outfitted with instruments to study wave spectra, wind 


stress, near-surface vvind profiles and vvater turbulence. The Naval Research 
Laboratory added tvvo radar scatterometers (14Ghz and 5 Ghz) as elements of 
the 1986 and 1987 experiments in order to take advantage of the opportunity 
to obtain collocated radar and environmental measurements. In 1986, both 
radars vvere mounted on the tovver railing at fixed incidence and azimuth 
angles. In 1987, the 14Ghz radar was mounted on an antenna rotator and 
extended on a boom away from the tower. This configuration allowed the 
radar to be rotated in azimuth at different incidence angle settings. 

The objective of this study is to examine the effects of the parameters 
shown in Figure 1 on measured scatterometer returns from an analysis of the 
WAVES data and to determine if these effects are successfully predicted by 
currently available model functions. Explanations are offered for the observed 
and predicted behavior of the radar cross-section. 

3. Limitations of the Study 

It should be noted that the above list of parameters which affect the 
cross-section is not complete (Masu'.o, et al., 1986). In effect, any oceanic 
phenomenon which has a surface expression may alter the cross-section; e.g., 
surface slicks caused by collection of surfactants by internal waves or 
Langmuir circulations (Huhnerfuss, 1983; Hughes and Gower, 1983), surface 
currents (Phillips, 1981), and the surface adjustment due to topographic 
interactions in shallow coastal areas (Valenzuela, 1985). Similarly, 
atmospheric phenomena such as planetary boundary layer rolls and or coastal 
fronts can be equally important in generating spatially varying surface 


signatures (Geernaert, 1990). 


Due to the location of the study site in Lake Ontario, the atmospheric 
and oceanic variations described in the above paragraph vvill not be addressed 
in this study. The maximum wind speed observed, 15 m/s, is not as great or 
sustained as would be possible in the open ocean. On Lake Ontario, the 
longest wave period observed is 8s and "wind rows" aligned with the wind, 
believed to be due to Langmuir circulations, are the most common form of 
surface slick. The fetch limits are the geographic bounds of the lake: 1 to 
300km. The wind direction is primarily from the west; i.e., the short fetch 
direction. Finally, the water is fresh versus saline. At the radar frequencies 
used in the experiment, the dielectric constants of fresh and saltwater are very 
nearly the same, therefore, the cross-section will not be affected. However, 
residence time of foam in a breaking wave regime is longer in saline water, 
which could cause the backscatter observed in the ocean to differ from that on 
the lake, especially at higher incidence angles and wind speeds (Gucinski, 


1986). 


C. NAVAL RELEVANCE: REQUIREMENTS FOR THE WIND VELOCITY 
AT THE SEA SURFACE 


Marine wind is the only environmental parameter listed as a 
requirement in all Naval Warfare areas (Brown and McCandless, 1988). In 
terms of criticality, only sea-state is more important although the correct 
specification or prediction of the sea-state has accurate marine winds as a 
necessary first requirement (Janssen, et al., 1989). Both strategic and tactical 
naval operations are affected by marine winds. ` 


The surface wind is an important parameter in strategic battle 


management. For example, the general weather picture (e.g., winds, storms, 


rain) is essential for decision making, resource management and optimum 
track ship routing. The systematic observation of ocean areas for effective sea- 
control and surveillance is hindered by high winds and waves which degrade 
sonar array performance by increasing ambient noise and reduce radar 
effectiveness by increasing sea Clutter. 

All operations at the tactical level can be suspended on account of 
extreme wind and sea-state conditions. In the traditional sense, these tactical 
Operations or decisions require only a local measure of environmental 
conditions. However, today's battlegroup size is O(1000 km) which extends 
the battle horizon to ocean basin scales when the deployment range of 
modern weaponry is considered. Therefore, next-generation tactical decision 
aids will have expanded areas of coverage. This concept is already being 
implemented in the Tactical Environmental Support System (TESS). 
Shipboard users of this system will use direct readout from overhead satellite 
passes and conventional data to produce local data products. These products 
will be supplemented by those produced at central site locations using all 
available satellite data, model outputs and conventional data. 

The pervasive requirement for surface wind coverage at spatial scales of 
tens to thousands of kilometers argues for global satellite measurements of 
surface winds. An Operational Requirement (OR) for this capability, which 
listed the performance goals shown in Table I (Daley, et al., 1984), was issued 
by the Chief of Naval Operations in 1977. This OR has since been rescinded in 
order to consolidate requirements for environmental Һс but the 


performance goals and requirements are still valid. The 


TABLE I. PERFORMANCE GOALS FOR SATELLITE MEASUREMENT OF 
SURFACE WINDS 


NOTE: The double entries indicate the most desirable values on the left and the minimum 
acceptable thresholds on the right of the slashes. 





Surface Data Receipt Rate Horizontal Measurement Absolute Range 
Wind Freq Timeliness Resolution Precision Accuracy 


_ Speed 3//12hrs 0.25//3hrs 10//25km 5//20% 2mps/ /Amps 1-75mps 
/ /3-25 mps 


Direction 3//12 hrs 0.25/ /3hrs 10//50km +5°//10° +250 0-360” 


most recent compilation of OR"s shovvs at least 24 requirements vvhich 


specify surface vvind as an important or critical element. 


D. SUMMARY 

İn this study, radar backscatter as a function of vvind speed and direction, 
long-wave slope, atmospheric stratification and water temperature is 
examined using data acquired from an instrumented tower in Lake Ontario. 
A historical review of scatterometry is presented in the next chapter along 
with background material on wave and air-sea interaction theory which is 
essential for the interpretation of the results. A review of current model 
function theory and supporting observations is given in the final section of 
Chapter II. A description of the experiment and data processing follows in 
Chapters HI and IV. The major results of the study are presented in Chapter 
V. These findings are summarized and the major conclusions stated in 


Chapter VI. 
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IL LITERATURE REVIEW 


A. INTRODUCTION 

To understand the physics of scatterometry, one must couple 
electromagnetic, ocean wave and air-sea interaction theories. It is the purpose 
of this chapter to summarize the aspects of each of these areas which are 
pertinent to scatterometry. The intent is to make clear how certain 
environmental parameters can (theoretically) influence the radar 
measurements and how currently proposed model functions incorporate 
these parameters. Each of the theory sections (C,D,E) derive, in outline form, 
the expressions used in the model function section (Section F). A synopsis of 
the history of scatterometry is presented first in order to provide a larger 


context in which to set the model function problem. 


B. HISTORICAL REVIEW OF SCATTEROMETRY 

In his experiment to test Maxwell's electromagnetic theory, Hertz, in 1886 
demonstrated that radio waves were reflected from various solid objects. In 
1904 Hulsmeyer (Ulaby et al., 1981) applied this result to the tactical problem 
of ship detection and in the war years to follow, ground-based radars were 
developed to detect aircraft and ships. These radars operated at meter or 
decimeter wavelengths, much longer that those of present-day radar 
scatterometers. To preserve the secrecy of the early radar development 
program, arbitrary letter designations, which are still in common usage, were 


assigned to the microwave frequency bands (Table IT). 


D 


TABLE II. MICROWAVE FREQUENCY BAND DESIGNATIONS 


Frequency Band GHz Central Wavelength Cm 


27100 
..- 5 
1.55 - 4.20 
4.20 575 
575 10.90 
10.9 


en 


Je 
L 
S 
C 
X 
K 
K 
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фә 


z «o 





By World War II, airborne radars were used to detect other aircraft 
and ships-at-sea, and some in the 3 and 10 GHz bands were producing images 
of the ground. As the instruments were refined, it became clear that radar 
reflections from ocean waves, which obscured the echoes from the actual 
targets (sea clutter), would have to be accounted for to improve data 
interpretation. Thus, the radars were directed at the sea surface itself with the 
purposes of improving radar performance, developing techniques to study 
the sea surface, and testing new theories which attempted to explain scatter 
from the sea (Stewart, 1985). 

In one cf the first controlled radar measurement experiments, 
Goldstein (1946) introduced the standard unit of measurement for radar 
backscatter from the sea, 0°, the (area) normalized radar cross-section (NRCS). 
Goldstein’s results suggested that o° varied with radar frequency, polarization 


and incidence angle (Huebner, et al., 1975). From 1959 to 1975, the Naval 
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Research Laboratory (NRL) continued the efforts to establish the character of 
the scattering cross-section via an extensive measurements program vvhich 
used airborne multifrequency radars to examine the dependence of o” on sea 
state and vvind speed (Guinard, et al., 1971, Daley, 1973). The determination 
that o” was in some way proportional to wind speed led to the proposal for a 
satellite-borne radar, by then known as a scatterometer, to obtain oceanic 
wind and wave information (Moore and Pierson, 1966). However, at that 
time, a successful correlation of o° with wind speed was yet to be established 
and the technology with which to build a stable, accurate (£2dB absolute) 
radar was not yet available (Jones et al., 1982). 

iuringaethe 1970's, “the National Aecronautiemand Space 
İD inistraion (NASA) sponsored scatterometer research and development 
which led to a technically improved aircraft scatterometer (Advance 
Applications Flight Experiment Radiometer-Scatterometer, AAFE- 
RADSCAT) and a proof-of-concept spacecraft scatterometer (Skylab, 5-193). 
With the RADSCAT instrument, Jones et al., (1977) demonstrated the 
azimuthal dependence of 6^; i.e., the cosine-like dependence of o with wind 
direction relative to the antenna. Using the radar data and corresponding 
in-situ data, empirical algorithms ("model functions") were derived to relate 
the radar backscatter to the ocean surface wind vector. In June 1978, as part of 
a suite of microwave instruments onboard Seasat, The Seasat-A Satellite 
Scatterometer (SASS) was launched. 

Three months into its planned three year mission Seasat suffered a 
failure in its power subsystem and communications with the satellite ceased; 


in November 1978 the mission was officially terminated (Pounder; 1980). Ten 


15 


years later the scatterometer data from this abbreviated mission are still being 
analyzed to confirm the instrument as a breakthrough in maritime 
meteorology and oceanography, but also to reveal the shortcomings in the 
scatterometer model function (VVoiceshyn, et al., 1986). 

Follow-on missions to the "proof-of-concept" Seasat were already 
planned at the time of its failure. NASA had proposed a Seasat-B mission in 
the summer of 1978. Early in 1979 there was an interagency (NASA, Dept of 
Defense, Dept of Commerce) proposal for a National Oceanic Satellite System 
(NOSS, Joint Effort by NASA, NOAA, Navy, 1979 internal report). Both of 
these systems were to be "high-heritage" and "limited-operational;" i.e., they 
would carry microwave instruments nearly identical to those on Seasat and 
no new, high risk (developmental) instruments would be flown. There were 
some changes proposed for NOSS: the scatterometer was re-designed to have 
three instead of two antennas, no synthetic aperture radar would be flown 
and a large-aperture multichannel microwave radiometer was to be 
developed. NOSS grew beyond its original charter but still survived internal 
review processes until 1981 when it was cancelled. 

Reiterating the need for an improved, all-weather oceanic data base, 
the Navy proposed N-ROSS (Navy Remote Ocean Sensing System), a 
supposedly "scaled-down" version of NOSS, in April, 1981 (Honhart, 1984). 
This satellite was to be the centerpiece of a constellation of ocean observing 
satellites to be launched in the early 1990's by the US, کت‎ and Japan. 
Funding and responsibility for the program was primarily slated by Navy 
and NASA; NOAA and the Air Force were to provide data acquisition and 


distribution support. Citing an unacceptable cost growth, the Navy cancelled 
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the program in December, 1986 (Graham, 1987, Matthevvs, 1987). The 
replacement launch vehicle for the N-ROSS NASA-Scatterometer (NSCAT) 
is currently planned to be the Japanese Advanced Earth Observing Satellite 
(ADEOS), tentatively scheduled for launch in 1995. Therefore, work still 
continues on the data processing algorithms including the model function. 
The European Space Agency will launch a scatterometer on its Earth 
Resources Satellite-1 in 1990. Looking to the next century, NASA scientists are 
proposing an alternate design for the Earth Observing Satellite system 
scatterometers which will use scanning instead of fixed antennas. A scanning 
system is more difficult to stabilize on orbit, but these antennas are more cost- 
effective and energy-efficient. Furthermore, they are able to operate as 
radiometers as well as radars so that it is possible to obtain corresponding 
passive microwave measurements such as atmospheric moisture parameters. 


With these changes the scatterometer will become a far more attractive 


‘instrument in a programmatic sense (Brown and McCandless, 1988), thereby 


increasing the chance that the US will have a scatterometer in space 


sometime during the next decade. 


C ELECTROMAGNETIC SCATTERING THEORY 

While scatterometer programs have been proposed and cancelled in the 
recent past, advances in knowledge of the scattering physics have been steady 
(Rice, 1951; Barrick, 1968; Wright, 1966, 1968; Valenzuela, 1968, 1978; Fung and 
Chan, 1969; Keller and Wright, 1975; Brown, 1978; Bahar, 1981; Durden and 
Vesecky, 1985; Plant, 1986; Holliday, 1986; Donelan and Pierson, 1987). To 


date, the three most common approaches to scattering problems are physical 
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optics, small perturbation and tvvo-scale composite surface methods!. Details 
of the derivations for these methods can be found in recent publications 
(Ulaby et al., 1982; Durden, 1986; Plant, 1988). Here, only the assumptions, 
results and experimental support for these theories are summarized. 
1. Fundamentals of Surface Scattering 

An expression for the NRCS (see Appendix A) is in terms of the 
electric field incident on the surface, Ej, and the field scattered from the 
surface, Es, for a given incidence angle, 8, azimuth angle, x, polarization, p, 


illumination area, A, and range, R: 





4лк? |Е(р)Г 
A |E (pJ) 





o°(0,x,p) = (1) 
i.e., the computation of the cross-section requires specification of the scattered 
field. Fundamentally, this means that the wave equation for E; is solved 
using Green's second vector theorem (See Appendix A, Beckmann and 
Spizzichino, 1987; Ulaby et al., 1982). This theorem states that the field at any 
point within a source-free domain bounded by a closed surface S can be 
expressed in terms of the tangential fields on the surface. 

For a surface rough in two-dimensions, the "Stratton-Chu" 


equations result as follows: 


E. FREA Har A(nAH)+(na E) Ts (2a) 


TAn alternative method is the generalized "full-wave" approach 
employed by Bahar (1981) and Brown (1978). Due to its generality, the 
approach and results are considerably more complicated and less easily 
interpreted in terms of the ocean wave spectrum. Therefore, this method has 
not been a widely adopted solution method. 


16 


FISS KE اوہ‎ K(i aH) + (ñ aE) parras (2b) 
0 


where a time factor of e!®t is understood; E, H are the total electric and 
magnetic fields on the surface, K is the scalar Green's function; R is the range 
to the surface; rs is the unit vector in the scattered direction; No is the intrinsic 
impedence of the scattering medium, and kem is the electromagnetc (em) 
wavenumber (Ulaby, et al, 1982). As required by Green's theorem, the 
scattered field is given in terms of the tangential fields, n^ Eand n AH. 

These equations simplify to the "Helmholtz Integral" for a surface 
rough in one dimension from which it is seen that the scattered field at 
Observation point r above the surface requires the field value and its normal 


derivative at the scattering surface 


1 OG OE 
F r= EG Hs (3) 
ит) 21 ön on 


In equation (3) G is the two-dimensional Green's function and the magnetic 
field H replaces E in the case of vertical polarized incident radiation. 

The essence of the scattering problem is the specification of the total 
field on the scattering surface. If the surface is irregular, then no general 
solution for the fields exist (Ulaby, et al., 1982) and three solution methods are 
possible: 1) the field is expanded as a series and the exact boundary conditions 
are applied (Brown, 1978; Bahar, 1981); 2) the integral equations are solved 


numerically (Durden, 1986); or 3) approximations of the boundary conditions 


on 


(EZ aae) are made, allowing closed form solutions to be derived as 
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outlined belovv. The simplifying assumptions generally include one or more 


of the follovving (Beckmann and Spizichinno,1987): 
1) (often) the surface is perfectly conducting, 


2) the dimensions of the scattering elements are much larger or smaller 
than the electromagnetic wavelength, 


3) the radius of curvature of the scattering elements is much larger than 
the electromagnetic wavelength, 


4) shadowing effects from larger elements are neglected, 

5) only the far-field values are computed, 

6) multiple scattering is neglected, 

7) the density of the irregularities (per unit area or volume) is ignored, and 


8) treatment is restricted to a certain model of surface roughness. 

These assumptions al! relate to the reflection properties of the 
scattering surface which depend on its electrical properties and roughness. 
The roughness of a random surface is radar wavelength dependent and is 
characterized by the rms surface elevation, h, and the surface correlation 
length, |. The appropriate approximations to use for a particular problem are 
based on certain limiting values of these two parameters. 

To be "smooth," the rms height must satisfy the Fraunhofer 


condition: 


A 


ҺГ 
00+00 (4) 


where A is the electromagnetic wavelength and @ is the incidence angle. This 
condition guarantees that there will be no more than a x/8 difference in 
phase between two reflected waves. For wavenumber k - 2x / 2, the limit for 
a smooth surface is kh < 0.3 (Ulaby et al., 1982) and the correlation length, I, is 


infinite since every point on the surface is correlated with every other point. 
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Reflections from a smooth surface are "specular" or mirror-like in that the 
angle of reflection equals the angle of incidence and very little of the radiation 
is scattered (Figure 3). Therefore, a monostatic (single antenna) radar must 
transmit at normal incidence to receive backscattered povver from a smooth 


surface. 


Incident 
Wave 


Scattering 
Pattern 


Medium Rough 





Figure 3. Examples of Surface Scattering Patterns 
2. The Physical Optics Approximation 
The near normal condition can be met for the sea surface at 
incidence angles less than about 20° since it is possible to have wave slopes of 
similar degree to maintain the required normal incidence. In the physical 
optics, or Kirchoff approximation, the roughened sea surface is approximated 


by small planar "facets," each locally tangent to the sea surface as shown in 


19 


Figure 4a. Each facet specularly reflects the incident radiation so that the total 
surface field is the sum of the incident and reflected fields 

E = (14 R)E, (5a) 
and 

dE 

r (5b) 
where R is the Fresnel reflection coefficient of the tangent plane at a 
particular point and r is the unit vector in the scattered direction. These are 
the boundary condition approximations for the physical optics method. 
These approximations are valid when the radius of curvature of the surface 
elements is much greater than the electromagnetic wavelength and the 
surface roughness is comparable to or larger than the electromagnetic 
wavelength (see Figure 4a). 

The cross-section is then computed as the sum of the reflected energy 
from all the facets oriented normal to the radar in the illuminated area using 
the probability distribution of surface slopes. Assuming an isotropic rough 
surface and Gaussian distributed surface slopes, Barrick (1968) determined the 
cross-section to be 
bt 


s” 


-[tan? в/5?) (6) 





0 
o ser" Oe 


where R(0) is the Fresnel reflection coefficient at normal incidence and s? is 
the total slope variance of the wave slopes in the x and y directions. If the 
surface is non-isotropic, then the cross-section in the upwind direction is 


written 


20 


2 
Ops = OY (uc en ers) (7) 


where sy and s; correspond to the slopes in the upwind and crosswind 
directions (Valenzuela, 1978). The slopes of importance here correspond to 
ocean wavelengths much longer than the microwavelength (to satisfy the 
criterion of small radius of curvature). The cross-section in (7) is termed 
quasispecular since it arises from only those parts of the sea surface which are 
normal to the incident wave vector. 
3. The Small Perturbation Approximation 

For a slightly rough surface and beyond about 30° incidence angle, 
comparable wave slopes are not physically realizable so that specular 
scattering is no longer valid, although scatter is still observed. At these 
greater angles it has been shown (Wright, 1966) that the response is from a 
resonance phenomenon between the em waves and the surface height 
displacements of the surface waves. This phenomenon, known as 
Bragg-scatter or Bragg-resonance, occurs when the em and ocean wavelengths 
(4, L) satisfy the equality: 


na 


= ДИН, 8 
2sin 0 (8) 





As shown in Figure 4b, when the excess distance from the source to each 


successive crest is nA/2, the round trip phase difference is 360” and the 
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(a) Tangent Plane--the Rough Surface as a Collection of Tangents 
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In-phase addition for Bragg scattering when AR = nA/2. 
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(b) Small Perturbation--Flat Mean Surface and Small Perturbation Imposed 


` FA OO" 
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NOTE: Separation of scales wavenumber is critical 


(c) Two-Scale--Small Perturbations Imposed on a Larger Scale Rough Surface 


Figure 4. Schematic Drawings of the Sea Surface for a) the Kirchoff (Tangent 
Plane or "Facet") Approximation, b) the Bragg-Resonance Model, and c) the 
Two-Scale Composite Surface Model Where Acm, Kem: Electromagnetic 
Wavelength and Wave number; h, Height above Mean Level; l, Correlation 
Length; L, Ag: Ocean and Bragg Wavelengths; AL Longer Wave of Two-Scale 
Surface 
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reflected signals add in phase. For example, for a radar wavelength of 2.18cm 
at 0 = 40°, the Bragg wavelength is 1.7cm corresponding to an ocean capillary- 
gravity wave. 

For such small waves it would seem that the resonant signal would 
be swamped by the return from the much larger waves. Average power 
return from the larger, randomly spaced waves is proportional to the number 
of the scatterers in an illuminated area. However, in the case of resonance, 
the average power is proportional to the square of the number of scatterers. 
Therefore, the resonant effect produces a significant power return, as well as 
providing a wavelength selective tool for examining the sea surface (Ulaby, et 
al., 1982). 

Assuming that the backscatter is from very short, small amplitude 
waves, the cross-section is derived using a perturbation approach in which 
the mean scattering surface is flat and small wave perturbations are 
superimposed (Figure 4b). "Small" perturbations are defined to have rms 
surface slopes which are sufficiently gentle (V2h/l < 0.3) and surface 
displacements small compared to the em wavelength (kh < 0.3).2 The total 
field on the surface is written as the sum of the incident and reflected fields 
(specular reflection from the flat mean surface) plus a small scattered 


component from the small amplitude roughness. 


?Chen and Fung (1988) suggest "that a more appropriate description for 
the range of validity of this model is a two-dimensional space with kh and kl 
as the axes"; ie. in addition to small height and slope the horizontal 
correlation length must also be large compared to the electromagnetic 
wavenumber. 
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The rough surface is considered random vvithin a particular L by L 
area but this random surface repeats itself in distance, L, i.e., is periodic. 
Flouquet's theorem is then invoked which states that a scattered field from a 
periodic surface is itself periodic. Therefore, it is appropriate to write the 
scattered field as a harmonic series expansion: 

p Y B, E(p +m,n;z) (9) 
.2x 2F 
where E(p+m,n;z) = S cae р = 2r /2 sin0 and kz is the 
vvavenumber component in the z direction. That is, the scattered field is 
periodic in x and y and propagates upvvard. 

The relationship between the coefficients, Bmn, for the different field 
components is found from Maxwell's divergence equation V- E = 0; .i.e., there 
are no sources in the region. The boundary condition of perfect conductivity 
(the tangential fields on the surface are zero) provides the additional two 
equations needed to solve for the three coefficients. Since it was assumed that 
kh was small all exponentials involving z are expanded in Taylor's series of 
kh about h = 0 and the coefficients are expanded as increasing powers of the 
surface elevation. For first-order Bragg scattering only the first terms are 
retained. The dominant effect of the first order expansion is to give rise to a 
z-component in the scattered field. Higher order expansions (Valenzuela, 
1968) produce cross-polarization effects. 

The radar cross-section for resonant scatter from a dielectric 


(nonconducting) surface is derived to be (Wright, 1968) 


o°(8) = 162k* cos* d Tat sin 9,0) (10) 
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where 
5 
SO تھی‎ (11a) 
Ë 0 4 (e, — sin ep 


15 (le * sin? 8) - sin? 6| 


Svv = (11b) 


2 
E cos 0 + (e, - sin? oy) 


and V (ky, ky) is the two-dimensional ocean wavenumber spectrum. The 
complex relative permittivity of the ocean, £r, is given by ве, = €/€p - 1O/€o0),in 
which the real part represents the capacity of the medium to store electric 
energy (c.f. a vacuum for which £592 8.85X10-l?Farad/m). The imaginary part 
represents losses by the medium through the conductivity, o, at frequency, o. 
The subscripts refer to the transmit/receive polarizations of the signal: h- 
pol/ h-pol and v-pol/v-pol. For a perfect conductor,6 e cə, £ç m -İ co, ghhe land 
gvv = (1 + sin”0) / cos?0. Therefore, resonant scattering of h-pol radiation 
from a good conductor is more sensitive to incidence angle than is v-pol. 

In his classic vvavetank experiment, VVright (1966) demonstrated that 
some of the observed characteristics of backscatter can be explained by this 
"first-order" Bragg theory. As shovvn in Figure 5, the measured backscattering 
and the polarization ratio (Oyy/Ohh) as functions of incidence angle agree 
rather well with the theoretically computed curves. Wright hypothesized 
that the differences from the theoretical cross-sections could be due to 
oversimplifying assumptions regarding the electric fields near the surface 


when no waves were present. 
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Backscattering cross section as a function of incidence angle for capillary waves on 
tap water. Left: the scattering section calculated from theory compared with observations. Right: 
the ratio of vertically to horizontally polarized backscatter, also compered with theory (from 
Wright. 1966). 





Figure 5. Evidence of First-Order Bragg Scattering 
(Stewart, 1985 after Wright, 1966) 


Further evidence for Bragg-scattering is found via analysis of the 
Doppler spectrum of the scattered signal. The Doppler spectrum is plotted as 
the spectral density of radar returns against Doppler frequency shifts "induced 
in backscattered microwaves due to a surface moving with a line-of-sight 
velocity toward or away from the radar" (Plant and Keller, 1989). If there are 
no background currents and Bragg scattering is the ad mechanism, 


then the Doppler spectra vvill contain peaks at those frequency shifts 


corresponding to plus and minus the Bragg vvave frequency. 
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Plant and Keller (1989) provide evidence of the Bragg scattering 
mechanism in field data. Doppler spectra computed from an L-band (16cm) 
scatterometer positioned at the end of the Coastal Engineering Research 
Center pier near Duck, North Carolina are shown in Figure 6. As predicted, 
spectral peaks appear at the Bragg wave frequency of 3.35Hz for the 
microwave parameters used in this experiment (Figure 6a). The wind 
direction for this spectrum is 210°, i.e., the offshore direction, and the wind 
speed is only 4.3m/s; hence, the dominant ocean wavelength will be 
correspondingly short. On the other hand, the peaks disappear for a spectrum 
corresponding to the onshore direction (Figure 6b). The longer waves have 
orbital velocities which contribute to the Doppler shift, thereby broadening 
the Doppler spectrum. It appears then that the first-order theory is valid until 
the dominant wind wavelength slightly exceeds the Bragg wavelength (Plant 
and Wright, 1977). At this point, the koh &41 and the mean surface is no 
| longer flat. 

4. The Two-Scale Composite Surface Approximation 

These limitations of the pure first-order theory led Wright (1968) and 
Bass et al., (1968) to propose independently a composite surface (two-scale) 
theory in which the surface is envisioned to consist of a large scale part which 
has small curvature and a small-scale part which backscatters according to 
first-order theory. The sea surface is broken up into an infinite number of 
slightly rough patches and the net backscattered power is computed as the 
summation of the power from a single patch distributed over’ the slopes of 


the dominant wind waves of the ocean (Figure 4c). The wavenumber 
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spectrum is decomposed into tvvo regions corresponding to the longer tilting 


waves and the shorter Bragg-scattering waves. 


Field Data Two-Scale Model 
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(a) Offshore Wind Direction 


30 Hz | 0 0 0 7 HE 





(b) Onshore VVind Direction 


Figure 6. Evidence for First-Order Bragg Scattering and Validity of Tvvo-Scale 
Model (a) (b) Left, Doppler Spectra from Field Data, (a) (b) Right, Doppler 
Spectra from Tvvo-Scale model 
(adapted from Plant and Keller, 1989) 


The radar cross-section is the sum of tvvo terms: the quasispecular 


contribution from the tilted surfaces and the Bragg contribution from the 


small amplitude vvaves. 
CHG) nie + J o,(8,y,6)p(tan y,tan 6)d tan 6 (12) 


where 
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oş(6,y,ö))  16nk" cos“ dioù. y, Ó) 





y,(2ksin(@ + y),2kcos(@+ y)sinó) (13) 


is the Bragg-scatter cross-section for a slightly rough surface tilted by yin the 
x-direction and ӧ 1п the y-direction. 

In (13), Ys is the short wavelength portion of the wave spectrum. 
The quasispecular term in (12) is now reduced by a factor exp(-4k2<h2 >) due 
to the small waves which make the specular "facets" less smooth, and the 
Bragg contribution is altered by the presence of the longer waves which tilt 
the Bragg waves in and out of the plane of incidence (Valenzuela, 1968). A 
separation wavenumber must be defined between the two assumed scales for 
the composite surface to be applied correctly. If such a "spectral gap" exists, 
then the scales are effectively decorrelated, allowing the linear summation 
of the two types of backscattered power. 

The relationships between backscatter and incidence angle 
corresponding to the three types of scattering models discussed are illustrated 
in Figure 7. Composite theory compares well with observed values of o° for 
vertical polarization, with an assumed ocean wave spectral form proportional 
to E but the curves diverge from the observed values for h-pol at higher 
incidence angles. Assuming higher, somevvhat unrealistic, background slopes 
(10° and 20°) helps to fit the data but the remaining differences imply that 
other scattering mechanisms may become important for h-pol scattering at 
grazing angles (Valenzuela, 1978; Donelan and Pierson, 1987). Better 
agreement (not shown) between observations and theory are found at lower 


microwave frequencies (.438 Ghz, L-band) than higher frequencies (4.4 Ghz, 
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Figure 7. Backscatter Versus Incidence Angle for Three Scattering Regimes 
(from Stewart, 1985 after Valenzuela, 1978) 
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C-band) (Daley, 1973). Finally, composite surface theory explains the widened 
spectral peaks observed in the Doppler spectra as being due to the advection of 
the small-scale structure by the orbital velocity of the longer waves (Figure 


6b) (Keller and Wright, 1977). 


D. SPECIFICATION OF THE ROUGH SEA SURFACE 

The expressions for the radar cross-section in the previous section show 
that the scattering cross-section is linearly dependent on the two-dimensional 
ocean wavenumber spectrum. To examine pure Bragg scattering with the 
first order theory, the spectral density for a small range of wavenumbers near 
the Bragg wavenumber is required (the high frequency tail of the full 
spectrum). Two scale composite surface theory also requires the joint 
probability density function of the upwind and crosswind slopes or, 
alternatively, the lower frequency gravity wave part of the wavenumber 
spectrum. Although the correspondence is good between the observed and 
computed backscatter, the absolute accuracies of the Bragg and two-scale 
formulations cannot be unambiguously established since a universally 
acccepted functional form of the two-dimensional wavenumber spectrum has 
not yet been derived. Thus, present model function research is necessarily 
semi-empirical and semi-theoretical: wave and radar data are used in a 
feedback loop with wave theory to derive new forms for the wavenumber 
spectrum and to determine the physical variables of importance to the 
development of the ocean waves, thence backscatter. It is at this point that 
the transition is made from a pure electromagnetic scattering problem to an 


ocean wave dynamics problem. 
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1. Ocean Wave Spectra 

The ocean surface is an evolving, random rough surface composed 
of waves which are generally categorized by their generating or restorative 
forces. The surface waves visible to the observer are most often the result of 
the driving force of the wind and the restorative forces of surface tension 
(capillary waves) or gravity (wind generated gravity waves). These waves 
span wavelengths of millimeters to hundreds of meters (wavenumbers of 
1000 to .001 m`!) and have periods of less than O(10s) (intrinsic frequencies of 
greater than 0.1 s!) Wave triads or quartets (capillary, gravity waves, 
respectively) can interact nonlinearly, with the result of an energy transfer 
among the waves, provided the wavenumbers of the individual waves satisfy 
certain resonance conditions (Hasselmann, 1962). Energy is removed from 
the wave system by wavebreaking, the generation of "parasitic" capillaries on 
the steeper gravity waves, or frictional dissipation at the highest frequencies 
(Hasselmann, et al., 1973). 

All of these competing generation and dissipation mechanisms lead 
to a random sea surface which is best described by its overall statistics, as 
formalized by spectral (Fourier) methods; i.e., as a linear superposition of an 
infinite number of sinusoidal components (Pierson, et al., p.24, 1955). This 
requirement is satisfied if the components are of small amplitude and slope 
and can be approximated by a sinusoidal profile. Under these circumstances 
the Fourier transform of the surface elevation covariance is computed as a 
function of wavenumber and frequency, with position E time as 
parameters. This "power spectrum" is the spectral "signature" of the sea 


surface. Donelan, et al., (1985) point out that the higher-order effects of 
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deviation of the vvave profile from a pure sinusoid, nonlinear interactions 
and viscous dissipation vvill prohibit the signature from being permanent, but 
to first order it is an accurate representation of the vvave field. 

Assuming the wave field is spatially homogeneous and statistically 
stationary, then the ergodic theory applies: space and time averages from a 
single time series (realization) are equal to ensemble averages across all 
possible realizations (Bendat and Piersol, 1986, p. 144). This random ocean 
surface is then characterized by the covariance, p(r, t), of the surface 


displacement, G, at points separated by a distance r(x,y) and time, t 


p(r,t) ex e(x,ty)e(x * r,ty * t) » (14) 


where the angle brackets denote an ensemble average. The covariance 
quantifies the decorrelation scales of the sea surface in space and time. There 
are then three spectra of varying dimension that can be derived from this 
‘covariance (Kinsman, 1965). 

The three dimensional (ky, ky, œ) wavenumber-frequency spectrum 
is defined as the Fourier transform of the covariance (Wiener-Kinchin 


theorem, LeBlond and Mysak, 1978, p. 311) 
1 | 
5(К,о)= —— ДИ (15) 
and inversely, 


p(r,t)= || [S(k, ae" ۶۹ (16) 


At r=0 and t=0 equation (15) gives the variance < C2 > of the surface elevation, 
or mean-square elevation. The calculation of S(k,@) requires simultaneous 


measurements of surface displacement in time and space, an infrequently 
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performed task vvhich requires an array of vvave staffs and intensive data 
processing (Donelan, et al., 1985). 
The two-dimensional wavenumber spectrum required for the NRCS 


calculation is the integral of the three-dimensional spectrum over all 


frequencies 
¥(k) = [S(k, @)do (17) 
At t=0, the covariance is 
p(r) = p(r,0) = [f ¥(K)e tak (18) 
then 
ے۔‎ ۷۲٢ رین‎ 


This spectrum is also difficult to measure since it requires surface elevation as 
“a function of space, an inherently more difficult in-situ measurement to 
make which involves the use of either stereophotography (Shemdin, et al., 
1988), laser optic techniques (Tang and Shemdin, 1983; Keller and Gotwals, 
1983) or a spatial array of sensors (Donelan, et al., 1985). 

It is much more common to have only a single wave staff acquiring 
data in time. In this case, a one-dimensional frequency spectrum results 


which is independent of wavenumber magnitude or direction 


O(a)- İİS(k, o)dk (20) 
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p(t) = p(0,t) = (f G(o)e "do (21) 
and 
p(t)e" dt (22) 


Or 


О(о) = = { p(t)cos cotdt (23) 


0 


since Q(w) is symmetric about w=0. 

An approximate wavenumber spectrum can be determined from the 
frequency spectrum using the dispersion relation, @=@(k). The wavenumber 
spectrum is often expressed in polar coordinates as the product of a 


wavenumber magnitude and direction spectrum 


W(k) = (Кр) = 5(Е)Ф(Ф) (24) 





where k = (kcoso, ksin o), k = |k|, and q is the wavenumber bearing (Pierson, 
1955). The wavenumber spectrum is obtained from the frequency spectrum 


using the dispersion relation from which 


1 dw Q(w) 
sü ee (25) 
[S CC wn 
(LeBlond and Mysak, 1978). This procedure cannot be used to determine the 
large wavenumber portion of the spectrum since the intrinsic frequency of 
the shortest (capillary) waves, Wo, is Doppler shifted by the orbital velocity of 


the longer waves on which they ride. The orbital velocity is essentially a 


зә 


slowly varying current of velocity U(x,t) so that the frequency of the short 
waves is given by 


w(k) = @)(k)+k-U(x,t) (26) 


That is, the relationship between frequency and wavenumber is no longer 
unique (Atakturk and Katsaros, 1987). 

Higher order statistics can be obtained from the wavenumber 
spectrum. For example, the mean-square slope of the waves is given by the 


second moment of the wavenumber spectrum (Cox and Munk, 1954) 

s - [IV (k)dk (27) 
The slope spectrum is often used in the study of the high frequency portion of 
the wave spectrum since higher wave numbers are weighted more heavily. 


Cox and Munk (1954) found that the mean-square slope is wind speed 


dependent. For a clean surface and windspeed at 12.5m in m/s, 


-2٦ )28(‏ ٭8 


and 


6 = 10?(84 1.56U) (29) 


for a slick-covered surface. In the Cox and Munk study, the slick was 
artificially generated and composed of diesel fuel and fish oil; the observed 


Slopes may differ with typical concentrations of natural films. 
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In polar coordinates, with the x-axis pointed upwind, the 


mean-square slopes in the upwind/downwind and crosswind directions are 
52 = {J (kcose) w(k, o)dkdo (30) 


and 
s? - İİ (ksin g)”v(k, p)dkdg (31) 
y 7 
These quantities are of particular interest in the present study since they are 
required for the computation of the quasispecular and Bragg scattering terms 
of the two-scale model and for the expressions describing the hydrodynamic 
modulation of the Bragg waves by the longer waves. 
2. The Spectral Transport Equation 
The adoption of spectral forms in wave analyses is a statistical means 
of revealing basic structures in the random wave field. Underlying these 
Statistics are, of course, the physical mechanisms responsible for the waves. 
The individual wave components generate, propagate and decay at different 
rates so that the spectral shape changes as the wave field develops. LeBlond 
and Mysak (1978, p. 322) then pose the question 
...aS to why spectral (and other statistical properties) take the shapes they 
do at various stages in the life of a wave field. Is it possible to derive 
predictive (as opposed to purely empirical or phenomemological) 
theories for the statistical properties of varying wave fields? More 
specifically, given the statistics of the forcing functions (such as wind, the 
atmospheric pressure), the boundary and initial conditions and the 
coefficients of the governing differential equations, each of which will, 


in general, have a random component, can we compute the solution and 
its statistics in terms of the statistics of the forcing? 
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The "radiative transport theory" of the evolution of spectra, vvhich 
arises in the area of stochastic fluid dynamics (Monin and Yaglom, 1965,1967), 
provides the mathematical formalism to ansvver these questions. The 
application of the theory (Hasselmann, et al., 1973) rests on the assumption 
that the time and space scales of spectral evolution are much longer than the 
longest resolvable period and vvavelength in a data record. This assumption 
allovvs the evolution process to be vveakly nonstationary. İt is found in the 
absence of currents that the ocean vvavenumber spectral density is conserved 


along a vvavenumber path except for certain exchange processes, denoted by X. 


a 
E 
di (32) 
or, 
2555575: 
ee 77٤ )33( 


These exchange processes include the energy input by the wind Qj; 
Qn], the redistribution of energy due to nonlinear interactions and the 
generation of parasitic capillaries on the steeper gravity waves; and Qg, the 
dissipation of energy by wave breaking and frictional dissipation,. The group 
velocity of the waves, cg = dw/dk, is the velocity at which the wave energy 
propagates. 

A more general expression which allows the waves to be 
superimposed on a space- and time-dependent mean flow, U(x,t) is in terms 
of the wave action density,N(k;x,t) 


dN _y 


Fis (34) 
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or, 


кии 


1 
E? “1 ok C+ U) VN =a (O tOu +Q) (35) 


where N = ¥(k) / wo. The frequency is now defined by equation (26) where 
Oo is the intrinsic frequency of the wave in a coordinate system which is 
moving with the current, U. In the presence of such a current, the wave 
action is conserved, whereas the spectral density is not. 

A visual summary of these equations is given in Figure 8 which 
shows a typical frequency spectrum and the associated source terms 
(Hasselmann, et al., 1973). Wind energy is primarily input at and above the 
frequency of the spectral peak and removed by dissipation at frequencies 
greater than approximately three times the peak frequency. The nonlinear 
terms transfer energy from the spectral peak to both higher and lower 
frequencies, but at the higher frequencies it is cancelled by the dissipative 
effects. Different stages of the wave spectrum evolution can be evaluated by 
examining different balances among the terms in the radiative transfer 
equations (32) and (34). When dV /dt 50, the wave spectrum is in a growth 
stage, a fully arisen sea is represented by d‘¥/dt=0 and the spectrum is in a 
decaying state when d'¥/dt<0. 

3. The Fully Arisen Sea, d‘¥/dt = 0 

The wave spectrum cannot grow without limit. Instead, for a given 
constant wind speed, it will reach some maximum state dictated by the length 
of time the wind blows (duration) and the distance over which the wind 
blows (fetch). At this point dV/dt-0, which implies that the source and 


dissipation terms are in equilibrium at all frequencies 
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— JONSWAP 


Figure 8. The Mean JONSWAP Spectrum and Source Functions 2’() 


together with the Computed Nonlinear Energy Transfer Z, | 
(from LeBlond and Mysak, 1978, after Hasselmann, et al., 1973) 


Q) + О, + Q, اد‎ )36( 


At wavenumbers below the wavenumber of the spectral peak, k < kp, the 
phase speed of the long waves is greater than the wind speed (for a fully 
developed sea) so that the waves lose some of their energy to the atmosphere 
only to have it restored by nonlinear transfer of energy from the spectral peak. 
These long waves are still in a growth stage and can accept more energy 
without breaking. Near the spectral peak, k = kp, wind input is transferred to 
higher and lower frequencies by nonlinear interactions. Above the spectral 
peak but below the wavenumber where viscous effects are predominant, 


l5kp«k« 3kp, (Leykin and Rozenberg, 1984; Donelan, et al., 1985), input 
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from the vvind and nonlinear interactions are balanced by dissipation 
processes, including wave breaking and the formation of "parasitic" 
capillaries ahead of the steep wave crests of the longer primary waves. 
Phillips (1958) identified this part of the equilibrium spectrum as the 
"saturation" range: the spectral density of these wavenumbers is limited by 
the onset of dissipation processes. | 

Observed spectra, appropriately normalized, exhibit remarkable 
similarity with the following characteristics. The wavenumber of the spectral 
peak decreases with increasing fetch and increasing wind speed and the 


spectral density increases with increasing fetch and wind speed (Figure 9a) 
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Figure 9a. The Evolution of Surface Wave Spectra with Fetch, from the 
JONSWAP Observations. The Fetch Increases with the Number labelling the 
Spectral Peak. 

(From LeBlond and Mysak, 1978, after Hasselmann, et al., 1973) 
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(Kitaigorodskii, 1962, Hasselmann, et al., 1973). There is very little energy 
belovv the vvell-defined spectral peak, vvhile above the peak, the spectral 
density decreases more slowly (Figure 9b) (Donelan, et al., 1985). Slope spectra 
(Cox, 1958; Tang and Shemdin, 1983; Shemdin and Hwang, 1988) also display 
characteristic features. These similarities have led researchers to compute 


parametric fits to the spectral shape. 





Figure 9b. Normalized Frequency Spectra Grouped into Classes by Ucp. The 
Vertical Bars at the Top of the Figure are an Estimate of the 90% Confidence 
Limits Based on the Standard Error of the Mean. 

(from Donelan, et al., 1985) 


The analytic fitting functions of the earlier models (before 1977) are 
based on Phillips' concept of a saturation range defined above. Phillips (1958) 
argued that wavebreaking in the saturation range was a dynamic instability in 
which the local acceleration of the fluid parcel exceeded the gravitational 
acceleration, g. Therefore, the spectral density was defined as a function only 
of frequency (or wavenumber) and gravity. Since the wavenumber spectral 


density has units of L4, by dimensional considerations we have 
SEET (37) 
or, in terms of frequency 


О(о)= аеја” (38) 


where «a and f are universal constants. The spectral slope is either k-4or c7? 
and independent of wind speed in the saturation range. Equation (38) is used 
as the "kernel" of the spectral model proposed by Pierson and Moskowitz 
(1964) and the Joint North Sea Waves Project JONSWAP) model proposed by 
Hasselmann, et al., (1973). 

Phillips (1977) modified the form of the saturation range to include 
the possibility of "microscale" breaking of the short waves (O(10cm)) in the 
presence of a significant wind drift layer, produced by the tangential 
component of the wind stress (Banner and Phillips, 1974). Banner and 
Phillips proposed that the drift current augments the phase speed so that 


waves whose phase speeds are a small multiple of the Surface nit preak 
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before reaching their maximum amplitude’. The criterion upon which the 


surface drift is assumed unimportant (after Wu, 1975) is 


2 
ku: <<2 (39) 
8 


where ux = 7 /p is the "friction velocity" of the surface drift layer for an 


applied wind stress tx. The spectral form is then 


S(k) = UR (40) 
8 
or 
— o 


Therefore, in the presence of a vvind drift layer, the slopes of the spectrum are 


k? or o” 


, assuming f1 and f» are linear. 

Phillips (1985) amended his concept of a hard saturation limit, based on 
observations that the constants, a and B, and the slopes in equations (40) and 
(41) varied and indicated a wind speed dependence (Leykin and Rozenberg, 
1984; Donelan, et al., 1985; Geernaert, et al., 1986). With the additional 


dependence of wind speed included, the functions in equation (40) and (41) 


are evaluated as 


3 However, Plant and Wright (1977) present empirical results which 
show the maximum wave height attained by a given wave depends only on 
its wavelength, i.e., augmented wind drift is not the limiting factor for the 
growth of the short gravity waves. 


-1 -7 


Noe 
S(k) 2 B|cos p|2u.g ? k? (42) 
Or, 


Q(@) = au.go™ (43) 


Kitaigorodskii (1987) suggests that the equilibrium spectrum changes from the 
wind speed dependent form (k-3-) to the saturated form (k-4) at about twice 
the spectral peak. However, scatterometer returns at Bragg wavenumbers 
above 2kp still display wind sensitivity. 

The spectral forms based on the equilibrium range are valid in 
wavenumber up to the wavenumber where surface tension, T, becomes 
important k<<Vpg/T (Kinsman, 1965). Beyond this wavenumber and before 
that where viscous dissipation becomes dominant, the governing variables 
are the surface tension, gravity, water density, frequency and wavenumber. 
Radar backscatter is sensitive to wind speed in this wavenumber range so that 
the spectral form should also display a similar dependence. The waves at the 
very highest frequencies are damped by molecular viscosity. As stated before, 
direct observations of frequency spectra in this range are broadened by the 
Doppler shifts in frequency induced by background currents. 

Pierson and Stacy (1973) used field and laboratory data to develop a 
spectral model which is broken into wavenumber ranges of applicability as 
specified above. More recently, Shemdin and Hwang (1988) made 
measurements of the high frequency slope spectra with a 55- sensor 
in order to evaluate the Pierson and Stacy (P&S) model (and others) as well as 


produce a new empirical model. Their results suggest that there are five 


spectral regions with frequency bounds as shown in Figure 10. The regions 


and their associated frequency and slope spectra (Ф(о)), s(œ)) are: 
1) the inertial range (region I) where 
Q(o) - au.go* (44a) 
or in terms of the slope spectrum, 
so) » ag u. (44b) 


2) the gravity range (region II) where 


О(о)= Ве?о” (45а) 


slæ) = Bo! (45b) 
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to Peak Energy in the Wave Height Spectrum) 
(from Shemdin and Hwang, 1988) 


Figure 10. 
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3) the surface drift range (region III) 


SIE E (46a) 
Ыо) = iə (46b) 
$ 
4) the capillary range (region IV) 
0 a 
Q(w) = «yo г (47а) 


S "ds (470) 
5 
where v is the kinematic viscosity, y is the ratio of surface tension to water 
density, and f» is a function defined by Lleonart and Blackmann (1980), and 5) 
the viscous range (region V) for which a functional form is not yet known. 
An alternative method for determining the high frequency portion 
of the spectrum is to invert the two-scale radar backscatter equation for the 
functional values of the short wavenumber spectrum (Lawner and Moore, 
1981, Chaudrey and Moore, 1983). It is al^^ possible to suggest a model for the 
spectrum and use the radar measurements to adjust its coefficients (Donelan 
and Pierson, 1987). This inversion method is predicated on the the existence 
of an accurate electromagnetic scattering model. 
Ihe ere spectral forms described above suggest that 
comparison of radar cross-section model functions reduces to a comparison of 


wave spectra. The rationale for this statement is shown in Figure 11 which 
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compares the (significantly different) frequency normalized slope spectra for 
the Pierson and Stacy,(P&S) Shemdin and Hwang (S&H), and Donelan and 
Pierson models (D&P) (Shemdin and Hwang, 1988). Shemdin attributes the 
differences between the P&S model and the S&H model to inadequacies in 
the data set used to generate the P&S model. Differences between the D&P 


and S&H models in the higher frequency region are attributed to the 


inversion method used by D&P. 
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Comparison of the Proposed Comparison of the Proposed Model 
Model with the Donel:n and Pierson with the Pierson and Stacey [1973] 
Pierson [1987] Model Model 


Figure 11. Slope Spectra Model Intercomparisons from Shemdin and Hwang, 
1988. 


4. The Source Terms 
The examination of the spectral forms in the previous section 
treated the entire right-hand side of the transport equation only qualitatively 


as a single entity (see a further discussion in Phillips, 1987). Theoretical 
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treatments of the balance of terms are provided by Kitaigorodskii (1983) and 
Phillips (1985). Explicit forms of the individual terms are novv briefly 
considered. A review of the literature pertinent to each term is given by 
Donelan and Hui (1990). Here, only a statement of the functional form for 
each term is given in order to recognize these quantities in the later 
discussion of the model functions. 


The wind input term is presently thought of as a sum of two terms 


Q, = a(k)+ B(k)¥(k) (48) 


where a(k) is a linear growth term predicted by Phillips (1957) wave 
generation theory and B(k)(k) is an exponential growth term due to Miles 
(1957, 1959). Phillips resonance theory applies only for the initial stage of 
wave generation when the surface is unperturbed. It is based on the 
hypothesis that fluctuations in air pressure due to atmospheric turbulence 
cause waves to be initiated. If the waves so generated move with a phase 
speed equal to the advective speed of the turbulent eddy, then the waves and 
pressure fluctuations will be in resonance. The forcing will continue and the 
results show that the waves will grow linearly in time. Once the waves have 
@eveloped sufficiently to perturb the air flow above them, then Miles’ theory 
of wave generation by shear flow applies. In this theory, the waves grow 
exponentially in proportion to the wave energy already present. 

Most wave generation studies show that the wave growth is 
exponential (Plant, 1982) except at very light winds (Kahma and Donelan, 


1988). Therefore, attention has been focussed on the specification of the 
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exponential grovvth rate B(k). VVave tank and field data have been used to 
detemine the growth rate in terms of friction velocity as (Plant, 1982) 
2 
Е il) cos @ 08<%<3 (49) 
o C C 


or wind speed at 10m (Hsaio and Shemdin, 1983), 


2 
B رر‎ d dı 
o p. 


C, 
or wind speed at one-half wavelength (Donelan and Hui, 1990) 


B в.Гид./2)сово | | 
—= (1 t0.2) | "HEY | (51) 


w 


where Cp is the phase speed of the wave component, 6 is the relative angle 
between wind and waves and pa/py is the air-water density ratio. Using (49), 


the wind input term ۶۶) ٢] 


` 


2 
Q, = 0.04cos pal) N(k) (52) 
C 
7 


The nonlinear interaction term, Qn}, in the transport equation (35) 
represents the redistribution of wave energy among wavenumbers. These 
interactions occur when several wavenumbers combine to form a new 
component. Phillips (1960) developed the resonance condition for gravity 
waves on deep water which occurs at third order (three waves combine to 


make a fourth) 


ki + Ka = ks + K4 (53) 
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In a continuous wave spectrum, an isolated resonant interaction produces a 
wavenumber which may interact with yet another set of wavenumbers. 
Therefore, the redistribution of wave energy becomes a conservative property 


of the entire spectrum; i.e., 


|a 


— *(44ЈАК, = 0 (54) 


Hasselmann (1962, 1967) derived the equation for the radiation balance to be 





d 
o 2 Ј//с18(4, + k, — k, — k,)öto, *0,—0,— c», ) Hk, dk;dk, (55) 


where G is a homogeneous cubic function of Ny, Ҹә, Ҹә, N4 and a 
homogeneous quadratic function of kj, kz, k3, and ky. Phillips (1985) then 
showed through dimensional analysis that the nonlinear term is 


proportional to the cube of the action density 


Qu = G”N?(k )k“or” (56) 


The final term in the transport balance, Qa, represents dissipation of 
wave energy, which results from viscous damping, wave breaking, 
turbulence and dissipation of parasitic capillaries at short gravity wavelengths 
produced by nonlinear interactions (Donelan and Hui, 1990). The form of the 
viscous dissipation rate of decay is known from classical wave theory (Lamb, 
1932) to be By = 4vk? where v is the kinematic water viscosity. The viscosity is 
an inverse function of temperature; therefore, explicit inclusion of the 
viscous damping will produce a temperature dependent spectrum. However, 
viscosity is only important at the highest frequencies when compared to wind 


forcing because of the quadratic dependence on wavenumber. Theoretical 


öl 


forms for the breaking of vvaves and generation of parasitic capillaries are not 
known. Phillips (1985) suggests that it is not necessary to know an explicit 
form for these processes, rather it is only required that the dissipation term 
balance the other source terms. Therefore, a dimensional analysis is written 
for the dissipation rate (Plant, 1986; Donelan and Pierson, 1987). Following 


Donelan and Pierson (1987) 





By _ 6, 
e 2 (Kk) E, yg] 7 (57) 


That is, the dissipation rate is a function of the wavenumber spectral density, 


wavenumber, surface tension and gravity, plus viscous dissipation. Then, 


b, 4 X^ | B. 
k 
– dd ¥(k); Ch - (58) 
Finallv, the dissipation term is written 
yk” 
Q, iron ç bəy (59) 


VVith the specification of the source terms, it is then possible to use 
the transport equation (35) to examine the interdependence of these terms 
and to model the wavenumber spectrum. The parameterizations of the 
source terms identify the most plausible environmental variables on which 
the wave spectrum could depend. These are the wind velocity or wind stress, 
surface layer stability, wind duration and constancy, fetch, water depth, 
surface tension, viscosity and water temperature. The selection of slightly 


different functional forms for these terms will produce correspondingly 


сл 
bə 


different radar cross-sections vvhen the resulting spectra are used in radar 


cross-section models. 


E. FUNDAMENTALS OF AIR-SEA INTERACTION THEORY 
The previous two sections showed that the ocean wavenumber spectrum 
was needed for the calculation of the radar cross-section and that this 
wavenumber spectrum was for the most part, dependent on the wind stress. 
The wind stress over (relatively flat) land can be directly measured using a 
drag plate of surface roughness comparable to the ground in which it is 
installed. No such direct measure of stress is possible over the undulating 
water surface; therefore, the stress is inferred from quantities which are 
directly measureable. The purpose of this section is to provide the standard 
mathematical framework which connects the wind stress to these measured 
quantities. This summary is a synthesis of similar treatments provided by 
Byers (1985), Arya (1988), Donelan (1989) and Geernaert (1990). 
1. Basic Theory 

The atmospheric planetary boundary layer (APBL) consists of the 
lowest 1000m or so of the atmosphere and is formed as a result of the 
interactions between the atmosphere and underlying surface on time scales of 
hours (Arva, 1988). In this laver, the vertical variation of the mean stresses is 
large compared to the horizontal gradients, therefore, turbulent shear eddies 
act to transfer heat, momentum and moisture to and from the surface at a 
rate much faster than simple molecular diffusion (Holton, 1979). 

The mean vertical fluxes of momentum, Uw Uw; sensible heat 
flux, Hs; and latent heat, EL, in the APBL are given by (using the notation of 


Geernaert, et al., 1988) 


(1 
Gə 





T= p|-i'wi = vwj (60) 
H. = pCpw' 6, (61) 


where 1 is the surface stress, pw is the water density, Cpis the specific heat of 
air at constant pressure, Ly is the latent heat of vaporization, OT is the 
potential temperature, q the specific humidity and u,v,w, are the 
longitudinal, latitudinal and vertical wind velocity components. Primed 
quantities represent fluctuations about a mean value where the mean, as 
denoted by an overbar, is taken over a time long enough for the mean 
vertical velocity to approach zero. It is assumed that the process is stationary 
and the ergodic theorem applies. These equations show that the stress can be 
computed from the covariance of the wind velocity components. 

The covariance measurements require a stable platform for accurate 
measurement of the vertical velocity and fast response sensors. Other 
methods for determining the wind stress (log-profile, dissipation) are 
generally neither suitable nor available for use on a standard basis. Therefore, 
bulk aerodynamic methods are often used to express the stress in terms of 


more easily obtained parameters; i.e., 


t, ٦ (63) 
H. 2 ۶ (64) 
E, = pL,C;Ul, (q, — 9) (65) 


where the numeric subscripts indicate the measurement height in meters 
above the surface. The bulk exchange coefficients are Cp for wind drag (the 
drag coefficient), Cy for the sensible heat flux (the Stanton number) and Cr for 
the latent heat flux (the Dalton number). The surface humidity is assumed 
Saturated with respect to the sea surface temperature. All of the coefficients 
are of magnitude O(10°). 

By convention, the wind stress is positive downward and the heat 
fluxes are positive upward. It is typical to define a "streamwise" coordinate 
system such that v and vu is assumed to be considerably less than "IT 
hence the stress is approximated by the x-component of (60) only. This 
approximation requires that turning of the flow by the Coriolis force can be 
neglected, which implies that the measurement be made close to the surface 
where friction has a greater influence. 

2. Surface Layer Scaling for Neutral Atmospheric Stratification 

The surface laver is defined as the lowest ca. 10% of the APBL in 
which the Coriolis force can be neglected. Furthermore, in this region the 
vertical gradients of the atmospheric parameters (u, O, gq) are significantly 
greater than the corresponding gradients of the fluxes. Therefore, the vertical 
gradients of the atmospheric parameters are scaled with the (assumed) 
constant surface fluxes (Businger, 1973). If viscous effects near the surface are 
ignored and the density stratification is negligible, then a similarity 
hypothesis is proposed that the velocity gradient depends on height above the 


surface, the surface stress, and air density; i.e., 


54 = Ы, тәр) (66) 


(л 
LA) 


where the surface stress accounts for all other influential parameters such as 


surface roughness. The only velocity scale which can be made from the above 


variables is the so-called friction velocity, u,, where 


n i (67) 
p 


The characteristic length scale has to be z. Therefore, from dimensional 


analysis 
2 (t (68) 
öz |) 
or 
ou lu 
r (69) 


where x is the von Karman's constant which is found from experiment to be 
about 0.4. Integrating equation (69) with respect to z (and now employing the 
constant flux approximation) vields the logarithmic wind profile, or as 
known in fluid mechanics literature, the "law of the wall" 
u(z) = CZ (70) 
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SCH , but zç is generally much less than z and (70) is used as 


EU 


Actually, it is In| 


an approximation. Here, the integration constant, the roughness length, zo, is 
a measure of the effective roughness of the surface and is introduced to keep 


the integral convergent. 
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A surface is considered aerodynamically smooth if its roughness 
elements? are small enough to allow a viscous sublayer to form near the 
surface and the elements are completely submerged in this layer. At the top 
of the viscous sublayer, z is not much larger than Zp viscosity dominates the 


stress and 


1 
ee. apr (71) 





where v is the kinematic viscosity of the air and g is the gravitational constant 
(Donelan, 1989). The thickness of the viscous sublayer for atmospheric flows 
is only 1mmy; therefore, the ocean surface must be "ultrasmooth" for a 
viscous sublayer to form. This may occur for very light winds (U « 3m/s). 
The typical sea surface is rarely so smooth and for wind speeds 
greater than about 7.5m/s, the roughness elements interact directly with the 
turbulent flow. In this case, form drag on the roughness elements accounts 
for nearly all the stress and the roughness length is proportional to the height 
of the roughness elements, as well as their areal densitv, characteristic shapes 
and dynamic response characteristics (flexibility and mobility) (Arya, 1988; 
Donelan, 1989). Charnock (1955) proposed that the roughness elements for 
fully rough flow were the short gravity waves, so that by dimensional 


analysis 


“ 
L 
4 “Roughness elements are surface features with sufficient steepness to 
cause flow separation on their lee faces and hence form drag. The flow may 
remain attached to quite large features having gentle slopes, but separate from 
small, abrupt roughnesses,” (Donelan, 1989). 


Un 
N 


5 7 mu. (72) 





where m is of the order of 10:7. Cardone (1969) also derived an empirical 
expression for Zo purely in terms of u», 


2. s 0.684 , 4.28x107u? — 0.443 (73) 
u. 





for u, in cm/s (See also Garrett, 1977). 


Donelan (1989) and Geernaert et al. (1986) revievved studies vvhich 
show that Charnock's relation is essentially an asymptotic form for Zo, 
appropriate only for a fully developed wave field. If the waves are actively 
growing, then zo appears to be more closely related to the wave parameters. 


Geernaert et al. derive a time dependent form 


ии 
21 el cut gi (74) 
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where X is the upwind fetch and t is time and zo 1s in meters. Donelan (1989) 


suggests that for fully rough flow and wind speed parameter, uj 


b 
Z H, 
" ox 5 (Z5) 


That is, the ratio of roughness length to the rms height of the vvaves is 
proportional to the inverse wave age, U/cp, where cp is the phase speed. This 
parameterization makes use of the fact that the longer (faster) waves are less 
steep. The coefficients, a and b, differ for differing wind speed parameters; i.e., 


wind speed at different levels or friction velocity, and field vs laboratory data. 


3. The Non-Neutral Surface Layer 

The previous section did not consider non-neutral stratification in 
the surface layer. Monin-Obukhov (1954) proposed a similarity hypothesis 
which includes buoyancy effects for a non-neutral, surface layer. The 
hypothesis is that in a horizontally homogeneous, quasistationary, "constant 
flux" surface layer, the mean flow characteristics depend only on the height 
above the surface, Kz, the friction velocity, us, the surface kinematic heat flux, 
w'T', and buoyancy flux, represented here as g/T, The simplifying 
assumptions are again that the Coriolis force can be ignored, viscous effects 
are negligible compared to the turbulent fluxes and the influence of surface 
roughness, boundary layer height and geostrophic winds are completely 
accounted for by u*. The result of the dimensional analysis is that the vertical 
profiles of wind, temperature and humidity are functions of a stability 


Parameter 7/L; e.¢., 





oz “527” E 

vvhere 
5 kzw' T, (77) 
i dti 


and L is the Monin-Obukhov length. 
The turbulent length scale, L, represents the relative importance of 
buovancy to shear effects. It is negative (positive) for unstable (stable) 


conditions and zero for neutral conditions. The nondimensional bulk 


9 


Richardson number has an equivalent interpretation as seen from its 


definition 





)78( 5 8م 
Tu‏ 
where ATy is the difference between the virtual temperatures between two‏ 
given heights (the air-water temperature difference in the present case). The‏ 
virtual temperature accounts for a moist atmosphere and is related to the‏ 
specific humidity by‏ 


T, 2 T(14 0.614) (79) 


The quantity z/L is a ratio of fluxes, hence it is difficult to measure. 
When measured directly it is typically "noisy" since the fluxes are second 
order moments which have large error variances associated with them (see 
Chapter III, Section C). Therefore, z/L is generally determined from the (first- 


order) bulk Richardson number and exchange coefficients by 


ENCE (80) 


(see Donelan, 1989; Hsu, 1989). 
The universal stability functions in equations (76) for the wind 
profile were measured over land (Businger, 1971) and found to be 
- al 
E (1— сети ии (81) 
L t 1 19 
where a and B were 15 and 47, respectively for x20.37. Panofsky and Dutton 
(1984) suggest 16 and 5 for these constants (for k=0.4). Kitaigorodskii and 


Donelan (1984) propose 17 and 5.4 (for x=0.4). 
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The diabatic log profile for the vvind is obtained from (81) by 


integrating from C5 = 29/L to Ç = z/L after the change of coordinates z = CL: 


mer “(| (82) 


where 
9 
4 
GEI e dD (83) 
SO 


Paulson (1970) solved the integral for unstable conditions, (<0 


1 22 5 
vde E aae (84) 
2 D 2 


where x = (1-15¢). For stable conditions €>0, and 


Hİ DE (85) 


A functional form for the bulk drag coefficient for height z is 


obtained from the combination of equations (63) and (67) 


"s 2 i 
с dis (86) 
и, = 


Here, the surface drift velocity, uo, was assumed to be negligible. The drag 





coefficient can then be interpreted physically as a function of measurement 


height, Z, wave state, Zo, and stratification, z/L (Geernaert, et al., 1988a). 


?»Geernaert (1986) notes that this omission of the wind drift amounts to 
only a 2 to 3% error in the wind stress. 


6] 


Many studies in the recent years have investigated the behavior of 
the drag coefficient with changes in environmental parameters (see 
Geernaert, et al., 1988a, for references). The results are usually normalized 


for stratification and reported in terms of a neutral drag coefficient, Cox, 
ku 
Con = Ë ə. z (87) 


where (87) is derived by assuming Zo to be constant for a given wind speed. 

The drag coefficient has a wind speed dependence since it depends 
on the roughness length (wave state). Therefore, most investigators plot 
their measurements against wind speed and parameterize the drag coefficient 
in terms of the wind speed, although there is a clear dimensional 
inconsistency. 

Some of the more frequently quoted drag coefficient 
parameterizations are illustrated in Figure 12. The physically based model 
functions discussed in the next section utilize the Large and Pond drag 
coefficient (curve 1) which specifies the neutral drag coefficient for the open 


ocean aS 


_ {0.114 LI € 10m / s 
Cox = C Си ири AG (88) 
since Lake Ontario typically has a shorter, steeper wave field, the 
drag coefficient is comparably higher as shown in the Donelan formulation 


for Lake Ontario (similar slope as Curve 7, Figure 12): 


Сры = .001(0.37 + 0.137) (89) 








Figure 12. Distribution of the Neutral Drag Coefficient vvith VVind Speed 
from the Following Investigations: 1, Large and Pond [1981]; 2, Smith [1980]; 
3, Smith and Banke [1975]; 4. Geernaert et al. [1987b]; 5, Geerneart et al. [1986]; 
6, Sheppard et al. [1972]; 7, Donelan [1982]; 8, Kondo [1975]. 
(from Geernaert, et al., 1988a) 


Therefore, for a given wind speed the wind stress over the lake will be higher 
than the same wind speed over the ocean. A proper comparison of the lake 
cross-section data to the scatterometer model function predictions as a 
function of wind speed, or stratified on wind speed, then requires that the 
neutral wind speeds over the lake first be "converted" to an equivalent 
wind speed over the ocean through an open ocean drag coefficient. 
4. The Thermal-Mechanical Internal Boundary Layer 

The Monin-Obukhov surface layer scaling described above strictly 
applies to an equilibrium boundary layer over horizontally homogeneous 
terrain. However, when such a flow over a uniform terrain encounters a step 


change in surface conditions such as roughness, an internal boundary layer 
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(IBL) develops downstream of the discontinuity as illustrated in Figure 13a 


(Rao, et al., 1974). 
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Figure 13. Schematic Diagram of the Internal Boundary Layer (IBL) over a 
step change of a) surface roughness (from Rao, et al., 1974) and b) temperature 
(Hsu, 1989). 


The research tower used in the present study is located 1.1km from 
the western shore of Lake Ontario, hence an IBL is formed when the wind 
blows offshore from the short fetch directions (approximately 180° to 360° 
True). The continental air masses blowing from these directions flow over 
the comparatively rough terrain of a midsize industrial city. The associated 
air temperatures often contrast significantly with the water temperature of 
the lake since the sources are either warmer air from the Southwest or cold, 


polar air from the North. Therefore, the IBL formed offshore can he 
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described as both "mechanical," due to a change in surface roughness, and 
"thermal," due to nonzero land-water temperature differences (Figure 13b). 
The height of this thermal-mechanical internal boundary layer (TMIBL) 
increases with distance from shore and its growth rate depends on the degree 
of contrast between the two regimes. 

The existence of an TMIBL is of concern when making wind stress 
measurements since the air flow above the TMIBL is representative of 
conditions upstream of the discontinuity, and the air flow within the TMIBL 
is representative of the downstream conditions. If a true representation of 
the wind stress over the water surface is the desired quantity, then the 
measuring device must be far enough from the land/water discontinuity and 
at a height which is within the TMIBL, in order to not be dominated by the 
upstream land conditions. Furthermore, the measuring device must be 
within the new "constant flux" layer (the height above the surface in which 
the fluxes are within 10% of their surface values) for Monin-Obukhov scaling 
to apply. 

Various parameterizations for the height and growth of the IBL have 
been derived for a change in surface roughness (Rao, et al., 1974; Claussen, 
1987), and a change in temperature (Mulhearn, 1981; Garrett, 1987; Hsu, 1989). 
For a change in surface roughness from rough-to-smooth, Rao et al. (1974) 
find that the height-to fetch ratio of the IBL height is about 1/10 and the new 
equilibrium laver ("constant flux") layer is about 0.5476 of the fetch; i.e., in the 
present case of a 1.1km fetch to the tower (shortest fetch encountered) the IBL 
is 110m high and the constant flux layer is only 5.5m. In the case of a thermal 


IBL, Garrett (1987) shows that the growth rate of the layer is slower for stable 


conditions since the turbulence is buoyantly damped. Hsu (1989) recasts 
Garrett's results to show that under stable conditions the IBL height is 16 
times the offshore fetch (in km); therefore, at the Lake Ontario tower for very 
stable conditions the IBL height is about 17.6m. 

These rough comparisons with the IBL parameterizations have 
several implications: 1) the anemometer used in the present experiment, 
which is 11.5m above the surface, is within the internal boundary layer but 
not necessarily within the "constant flux” layer as defined in Monin- 
Obukhov scaling. Therefore, the stress measurements may differ by more 
than 10% of their surface values, and 2) both the atmospheric boundary layer 
and the wave field are in "growth" stages when the wind blows from the 
short fetch directions in this experiment. The comparisons made above are 
for the worst case scenario; i.e., stable flow from the shortest fetch direction. 
A wide range of wind speed, fetch and stability conditions were encountered 
during this experiment, and alternative scaling arguments for the 
computation of the wind stress are not available, therefore, it was assumed 
that for the most part the anemometer was within the constant flux layer and 
Monin-Obukhov scaling was emploved to facilitate comparison with 


previous scatterometer studies. 


F. SYNOPSIS OF RECENT MODEL FUNCTIONS 

The previous sections in this chapter have provided a brief background 
on electromagnetic scattering, ocean wave and air-sea interaction theory 
required for the development of a scatterometer model function. The path 
through these intertwining theories led to the identification of 


environmental parameters which may affect the radar cross-section. These 
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parameters include: the vvind speed (or vvind speed parameter), vvind 
direction, long vvave field (long vvith respect to the Bragg vvave), atmospheric 
stratification (as it affects vvind speed) and vvater temperature (as it affects 
viscosity). Inasmuch as these quantities all relate to the ultimate form of the 
ocean wavenumber spectrum used in the two-scale scattering model, 
differences among the "physically based" model functions amount to 
differences in the specification of the spectrum. The purpose of this section is 
to outline the forms and predictions of some currently proposed model 
functions in terms of their treatment of the ocean wavenumber spectrum and 
their methods for including/excluding the environmental parameters above. 
This study is intended to be a demonstration of environmental effects on 
the cross-section, not a comprehensive test and evaluation of model 
functions. Therefore, the models are only discussed at a level of detail 
sufficient to give the reader the justification for comparisons with the field 
data and the conclusions to be drawn from these comparisons. The attributes 
of the different models are summarized in Table III at the end of this section. 
1. The Seasat Scatterometer SASS1 Model Function (1980) 

The genesis of the SASS1 model function is chronicled in Schroeder, 
et al. (1982). Three candidate model functions produced prior to the Seasat 
launch were based on Statistical fits of cross-section to the neutral wind speed 
at 19.5m using scatterometer data from aircraft circle flights (Jones, et al., 1977, 
1982). After the demise of Seasat, the three model functions went through 
refinements via a series of workshops in which the predicted wind vectors 
were tuned to additional in-situ data from the Joint Air-Sea Interaction 


(ASIN) experiment and compared to in-situ winds from the Gulf of Alaska 


SEASAT experiment, GOASEX (see Brown, et al., 1983 for a description of the 
in-situ data). The models attempted to reproduce the observed increase of the 
cross-section as some power of the wind speed and the approximately cos2x 
dependence of the cross-section on relative wind angle, x (recall Figure 2b). 
The result of these workshops was SASSI, a "power law" model 
function which relates the NRCS to the neutral wind speed at 19.5m, 


incidence angle and azimuth angle: 


o°(U, 621 ES 10İG(U, 0 2 +(U, 0, x)logio ew (90) 


where the wind speed is in m/s. The G and H coefficients are obtained from 
lookup tables which tabulate the coefficients as functions of incidence angle 
(0^ to 70^ in 2 increments) and azimuth angle (0 to 180 in 10^ increments. 
Two separate tables were produced for vertical and horizontal polarization. 

The entries in the tables were a synthesis of the three models and 
were tuned to eliminate errors between the observed and predicted winds 
(Wentz, 1984). Of the three models, the Wentz model, which used the aircraft 
and in-situ. data to determine the coefficients of a two-scale model was used 
for table entries between D and 12 incidence. Between 24 and 64,a 
combination of the three models was used (Schroeder, et al; 71982» 
Interpolation of the table filled the gap between 12° and 24° and extrapolation 
extended the table to 70° incidence (Boggs, 1981). 

The SASS1 model was used to process the three-months of SEASAT 
SASS data. Many preliminary studies successfully used the SASS1 wind 


vectors to demonstrate the utility of scatterometer winds (JGR Special Issue, 


Feb, 1983), even though problems in the model function vvere later 
uncovered, as described belovv. 
2. The Wentz SASS2 Model Function (1984, 1986) 
Wentz, et al. (1984) point out systematic errors in the SASS1 derived 
wind vectors, which are attributed to tuning the scatterometer cross-sections 


to a very limited in-situ. data set. As stated in Wentz (1986) the errors are: 
1) a positive bias between h-pol and v-pol derived winds, 


2) an artificial wind speed gradient across the satellite swath due to 
incorrect specification of the incidence angle dependence, 


3) a positive bias of wind speeds below 6m/s due to the exclusion of low 
signal-to-noise ratio NRCS values (the logarithmic form of the model 
function is undefined for zero and negative values of NRCS), and 


4) a positive 1m/s bias of the SASS wind speeds. 

Wentz (1984) notes that additional physical variables such as the 
viscosity, wind fetch and duration and long gravity waves which may be 
important as second order effects, will not significantly alter the power law 
relationship. As listed in the introduction here, Woicheshyn, et al. (1986) 
pointed out the above failings in the SASS] model function. However, they 
do blame incomplete specification of the physics involved in using the cross- 
section measurement to predict wind when a power law relationship is used; 
particularly the neglect of the effects of atmospheric stratification and the 
water temperature. Wentz (1986) counters that systematic departures from 
the power law have not been observed and that the errors in the SASSI 
algorithm were not due to the form of the power law but . the original 


coefficients used in the model. ( 
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Based on the observed povver lavv dependence of the cross-section to 
wind speed and the azimuth dependence with wind direction, Wentz (1984) 


fits the cross-section data with a truncated Fourier series as follows: 


o°= AQ + A, cos 7 ++ ظ4‎ 7 (91a) 
n Le (91b) 
A, = (a, + a, log U)A, (91c) 
A, 7 (a; * a, logU )A, (91d) 


where ag is the power law coefficient. The coefficients an and Gp are 
dependent on the polarization and incidence angle. Instead of fitting the 
coefficients to any particular data set, the coefficients were computed by 
assuming some Statistics about the global wind field: the wind speed as 
deduced from the SASS measurements are assumed to have a Rayleigh 
distribution, and a uniform distribution of wind direction®. The assumptions 
are validated by comparing wind speeds computed from the SASS1 model 
and the aircraft data used to produce the SASS1 model. The comparisons 
show better agreement between the SASS2 predictions and the aircraft data 
than similar comparisons with the SASSI predictions. There is still a 
systematic 0.5m/s bias between the h-pol and v-pol wind speeds which 


becomes greater with increasing incidence angle. 


The assumption of a uniform distribution of wind directions seems 
intuitively less defensible since the global winds have a definite latitudinal 
dependence. However, Wentz states that the predictions of the model 
function are not particularly sensitive to the actual distribution used; 
therefore, the uniform distribution is employed for the sake of simplicity. 
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3. The Durden and Vesecky Model (1985, 1986) 

The SASSI and SASS2 model functions accept the two-scale 
scattering theory as the physical basis for the observed cross-sections but the 
theory is not explicitly used to calculate the cross-sections. Furthermore, they 
are em frequency specific. In the Durden and Vesecky model (D&V) model, 
patterned after the Fung and Lee (1982) model, the radar cross-section is 
computed using the two-scale theory, as described in Section C of this chapter. 
The separation wavenumber is chosen such that both the curvature of the 
large-scale surface and the mean-square height of the surface are small in 
order to minimize the rms error in the two-scale model, as defined by Durden 
(1986). The use of the two-scale model requires a form of the ocean 
wavenumber spectrum, the definition of the separation wavenumber, long 
wave slopes and a probability distribution function of the slopes, and an 
atmospheric surface layer theory which connects the wind stress to the wind 
speed. The D&V model also requires water temperature and salinity for the 
computation of the dielectric constant in the polarization constants 
(equations (12), Section C). 


2 
The wavenumber spectrum is split at the wavenumber, K.=g/U gs, 


into a low wavenumber portion and a high wavenumber portion 
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whens 1+ = (93) 


is a modified gravitational parameter vvhich includes surface tension, T, at 
high wavenumbers. The low wavenumber portion is the Pierson- 
Moskowitz similarity spectrum (1964) which has a VI wavenumber 
dependence. Durden and Vesecky (1985) extend the Pierson-Moskowitz 
spectrum to include the effect of the wind drift layer on microscale breaking at 
the higher wavenumbers by making the wavenumber spectrum a function of 
u*^/g as in equation (39) of Section D. The functional dependence is assumed 


to be a power law in which the exponent is a function of wavenumber. 


2 r(k) 
S(k) = 2: (94) 


The exponent is modeled as a logarithm of the wavenumber, r(k) = log(k/K}), 
based on observations that the cross-section power law with wind stress 
varies with the radar frequency (Jones and Schroeder, 1978; Fung and Lee, 
2 

1982). The separation wavenumber, Kt, is the wavenumber at which سس‎ S 
is no longer valid; here K, is taken as 2. The wavenumber magnitude 
spectrum is illustrated in Figure 14 for U19.5 = 10m/s and upwind(o=0) and 
crosswind (0=90) conditions (Verma, et al., 1989). 

The directional dependence of the spectrum is the cos2@ term 
enclosed in square brackets of equation (92) and is displayed in Figure 15. The 


multiplicative exponential factor in equation (92) was introduced by Fung and 
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Figure 15. Spectral Density Image of Directional Dependence of Durden and 
Vesecky (Two-Sided) Wavenumber Spectrum 


for Wind Speed = 10m/s, phi = Wind Direction 
and k = wavenumber. 
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Lee (1982) to model the observed increase in cross-section azimuthal 
modulation vvith radar frequeney, the coefficient c must be determined from 
measured quantities. Here, c is determined by forcing the cross-vvind to 
upvvind slope ratio to be equal to the Cox and Munk (1954) value for a clean 
sea surface. The azimuthal modulation of the cross-section produced by this 
directional spectrum is symmetric in the up and dovvnvvind direction, vvith 
minima exactly at the crosswind directions.Hydrodynamic modulation of the 
small-scale waves by the large-scale waves is accounted for by introducing a 


modulation parameter, m, 


S 


Y. (k) 2 (14 ms, JV" (K) k > K, (95) 


where su is the large-scale slope in the upwind direction (Jones, et al., 1978) 
and v/ş(k) is the high wavenumber portion of the wavenumber spectrum 
(k>2). The parameter, m, can be adjusted to "calibrate" the model to a 
particular data set; i.e., hydrodynamic modulation produces the 
upwind/downwind asymmetry in the azimuthal modulation and shifts the 
minima downwind from the crosswind directions (Li, et al., 1989). 

The wind stress is related to the wind speed at 19.5m by the surface 
layer model of Large and Pond (1981) where the neutral drag coefficient at 
10m defined by Large and Pond is first converted to a value at 19.5m. The 
corresponding wind stress is computed from equation (82) with the air-sea 
temperature difference used to compute the stability correction function, V. 


The wind stress is then converted to a neutral wind speed at 19.5m by 
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assuming Zo has the functional dependence as derived by Cardone (1969) (73).7 
Fung and Lee (1982) point out that this equation is valid for U,95 >3.5m/s. 
The D&V model produces erroneous results for wind speeds below this 
value. 

Predictions of the cross-section with this model are reproduced in 
Figure 16a from Durden and Vesecky (1985). For a frequency of 13.9GHz and 
wind speed equal to 15m/s, the h-pol and v-pol calculated values are within 
3dB of the observations obtained from Guinard (1971). However, the 
polarization ratio Oyy / Oy is larger than observed. The predicted upwind 
wind speed "power law" exponents are compared to the observations in 
Figure 16b. The predicted exponents for both polarizations show reasonable 
agreement with the observations with the greatest discrepancy at 20°. 

4. The Donelan and Pierson model function (1987) 

The Donelan and Pierson (D&P) (1987) model function also uses the 
quasispecular two-scale scattering theory as its basis; therefore, the same 
quantities as specified for the D&V model are required here. Some 
differences from the D&V model should be noted about the computation of 
€, the dielectric constant, and R(0) the Fresnel reflection coefficient at normal 
incidence. The dielectric constant in the D&V model is computed directly, 
here it is assigned for a particular frequency, temperature and salinity from 


standard tables (Saxton and Lane, 1952). The polarization constants are then 


“This expression for Z, uses a different form for Cpx versus windspeed 
than that of Large and Pond. This discrepancy can be reconciled by 
recomputing the numeric constants in (73) using the data of Large and Pond. 
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Figure 16 (a) Cross-Section vs. Incidence Angle and (b) “Power Law” Wind 
Speed Exponent as a Function of Incidence Angle from 
the Durden and Vesecky Model (1985) 


computed at one degree incidence angle increments using the standard 
dielectric constant to build a lookup table which is then interpolated for 
particular conditions. The Fresnel coefficient at normal incidence is 
multiplied in the D&P model by 0.65, as suggested by Valenzuela (1978), while 
in the D&V model it maintains its full value. 

The full wavenumber spectrum in the D&P model is constructed by 


patching” together two independently derived spectra for the low and high 


wavenumbers, separated at a wavenumber kr. The low wavenumber end of 


the spectrum (for vvavenumbers less than ten times the spectral peak) is 
proportional to the mean vvind speed at 10m and vvas derived from an 


empirical analysis of vvave data (Donelan, et al., 1985) 
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(96) 
TT 0.5 2 
where F(U (10), k) = -—.- E | | 
3 
The peak of the spectrum, kp for fully developed seas is given by 
: g (98) 


ütü) 


The mean-square long wave slope is found by integrating this part of the 
spectrum using (27). Donelan and Pierson (1987) note that they were unable 
to reconcile the slopes determined from this spectrum to the Cox and Munk 
slope values; therefore, for internal consistency they used their integrated 
values. 

The high wavenumber spectrum (at the Bragg wavelengths) is based 
upon a proposed balance between the energy input to the short waves by the 
wind to that dissipated through wave breaking and viscosity. The growth rate 


(wind input) is modeled by 
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where U(n/k) is the mean wind speed just above the waves at a height of one- 
half the Bragg wavelength and c,(x/k) is the phase speed of the short waves. 
The normal stress is claimed to be more important in the amplification of the 
capillary waves (Kinsman, 1965) and is represented in (100) by the wind speed 
to phase speed ratio in the wind input term. The dissipation term is 
determined by the dimensional analysis outlined in Section D. The balance 


between wind input and dissipation then yields the high wavenumber 








spectrum: 
= : = 
0.194 u{ =) 4 vk 
V (k)o k* ——— £s و ارد‎ (100) 
ар, 42 ac, (k) 
P| k 


where a-f] (TR PRE) and n=f2(Tk*/ pw). These coefficients are used to match 
the two pieces of the spectrum together and they are obtained from gravity 
wave spectra and radar data. The model spectra, computed for water 
temperatures of 0° and 30°, and various wind speeds are shown in Figure 17. 
The directional dependence of the spectrum (Figure 18) is modeled 


as 


p(Q) = sech’[h,(g - o) (101) 


where @ is the radar look direction, @ is the peak direction of the downwind 
traveling waves and hj is an empirically derived value that controls the 
spreading of the waves about the peak direction. The value of h; decreases 


for increasing wind speed; ie., the waves spread out more about the wind 
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direction. Also, hı keeps the two-dimensional spectrum from going to zero at 
(р-ф1 =л/2. 

The short waves are modulated by the longer waves by writing the 
modulation as a function of the downwind slope of the tilting waves, sq, after 
Reece (1978): 

(1- ms, )¥(k) |ms,| < ; 


U = : (102) 
(1 — 0. 5sen(ms, It) İms,l > 5 


where m is determined empirically. The modulation will cause an 
upwind/downwind asymmetry in the computed cross-section values. Here, 
D&P set m=1 to fit the model output to observed upwind/crosswind 
differences in the NRCS. 

The remaining quantity to be specified is the wind speed parameter 
at height x/k, U(x/k)., where the value of m/k or A/2 is O(2cm). Citing 
theoretical and empirical work, D&P maintain that the logarithmic profile of 
the mean wind extends to this small height above the surface and that at 
these heights atmospheric stratification effects will be negligible since the air- 
sea temperature difference approaches zero*. Therefore, the law of the wall is 


applied to obtain 


95However, calculating 012: from U(10) should still require stability 


correction if (To - Tio) is nonzero, while 7 = A approaches Zero. 
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a(z) = U (10) "iid 


ËO —In b (103) 
vvhere the neutral drag coefficient can be either of the empirical expressions as 
found by Large and Pond for the open ocean (88) or Donelan, for the lake, (89). 

The radar cross-sections predicted using the D&P model increase 
with wind speed and are strongly influenced by water temperature at low 
wind speeds due to viscous dissipation of the short waves. In fact, a threshold 
wind speed dependence upon incidence angle and water temperature occurs 
at each wavenumber when dissipation balances input and the bracketed term 
in equation (100) equals zero. Below this threshold the viscosity inhibits the 
generation of the corresponding Bragg wave, the spectrum ceases to exist, and 
no backscatter is produced (see abrupt drop of the spectra in the lower part of 
Figure 17; see <2m/s Figure 19). Just above the threshold and up to 
approximately 10m/s for Ku band (Figure 19), the cross-sections over warmer 
water are several dB higher than over cooler water. Beyond 10m/s the 
influence of water temperature is not as apparent except for high incidence 
angles. At very high wind speeds (>30m/s) the cross-sections level off and 
"rollover," displaying a negative correlation with increasing wind speed. The 


n 
authors propose that the decay is due to the decrease of Uk) as U10 increases 


because of the increase in Cpn with Ujo. 
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Figure 17. Wavenumber Spectra from Donelan and Pierson Model (1987). 
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Figure 18. Spectral Density Image of Directional Dependence of Donelan and 
Pierson (One-Sided) Wavenumber Spectrum (from Verma, et al., 1989) 
for Wind Speed = 10m/s, phi = Wind Direction, and 

k = Wavenumber. 
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Vertically polarized backscatter for the two-scale model 
for K, band at upwind with log,, O(19.5) in meters per second on the 
abscissa and a,,.° in decibels on the ordinate. The solid curves are for 
a water temperature of 30°C and the dashed curves are for a water 
temperature of O°C. @ 1s the incidence angle. The auxiliary scale is the 
wind speed at 19.5 m in meters per second. 


Figure 19. Predicted Cross-Sections Against Wind Speed for Different 
Temperatures (from Donelan and Pierson, 1987). 


5. The Plant Model (1986) 

Plant (1986) noted that the existing model functions did not include 
an explicit accounting of the nonlinear interactions between the short Bragg 
waves and the longer dominant waves. To examine the possible dependence 
of o on the long-wave, short-wave interactions, Plant (1986) applies the 


principle of conservation of wave action to develop a spectral model for the 


short, wind generated waves. Defining short wave action as 


N(k,x,t) = pc(k)y(k,x,t) for the short wavenumber k, the governing equation is 


oN oN d oN 
۱ KL? 
(c, U) Ox ox l ) 


—- — = pN + D 104 
3 SC (104) 


where c, is the group velocity, U the horizontal current due to the orbital 
velocity of longer waves and D is a decay term which includes dissipation, 
short wave-wave interactions and wind-wave interactions. The form of D is 
found by dimensional analysis to be a function of k*'W(k) as in equation (59), 
Section D. 

The long wave modulation of the short waves is then modeled by 
writing the wave action as N = <N>(1+M) where <N> is the average action 
density and M is a term representing the modulation transfer function (MTF) 
of the long wave-short wave interactions. Equation (104) is ensemble 
averaged to yield an equation for the average action balance among long 
wave-short wave interactions, wind input and dissipation. Expanding <N> 
in a power series of s*, the mean square long wave slope, solving for the 
zeroeth and first order, and substituting the spectral definition for N yields 
the following form for the short wavenumber spectrum: 


HE S a 24 os + Y + Y,)cos(6, - d (105) 


where k, K are the short (long) wavenumbers; w, Q are the short (long) wave 


frequencies; m, is the imaginary part of the modulation transfer function; s? is 


the long wave slope; 0, 0, are the short and long wave directions and 


Yo = f(k, K, No) and yı = f(k, K, m). 
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To evaluate ¥(k) further requires a form for the wave growth f, 

chosen here to be (Plant, 1982) 
и? x 

B= —... , ت6‎ « 7 (106) 
0 , elsewhere 
The vvavenumber spectrum is functionally dependent on friction velocity, 
long vvave slope and the angular variation of the short vvave grovvth. The 
azimuth variation of the cross-section arises primarily from this cosine 
dependence in the grovvth term vvith a smaller contribution from the 
modulation term. 

Plant (1986) assumes neutral conditions and uses the neutral drag 
coefficient of Large and Pond (1981) to determine u, from equation (86) using 
a range of assumed wind speeds at 10m. The mean square slope of the long 
waves in equation (105) is estimated from the Cox and Munk (1954) values 
for a slick-covered surface, equation (29). The difference in wind speed 
between 12.5m and 10m is ignored?. 

The spectrum is used in the composite surface expression for O to 
yield cross-sections with incidence angle, azimuth angle, polarization and 
wind speed dependences. For the purposes of this preliminary model, Plant 
considers only the pure Bragg scattering regime; therefore, the model applies 


only at incidence angles between 20° and 60 incidence angles. The predictions 


9 The slope equation published in Plant (1986) and that used in the 
model program are inconsistent. Since the code uses the Cox and Munk 
values (equation 29) and the diagrams in the paper are produced from the 
code, it is assumed here and verified by Plant (pers. comm.) that (29) was the 
correct equation. 


67 


of Plant's model compare most favorably with SASS1, SASS2 and aircraft data 
at incidence angles between 30° and 50°. The predicted h-pol values are 
typically low and the crosswind values fall well below the data at 25° incidence 
angle and at low wind speeds (Figure 20a). The cross-section was found to 
increase with long wave slope but not as much as observed (Keller, et al., 
1985). Plant concludes that the most important variables to be understood for 
further progress to be made in scatterometry algorithms are short-wave 
growth rates and the surface wind stress. Since these are derived quantities, 
the fundamental variables implicated in this statement are viscosity, 


atmospheric stability and long wave slopes (Geernaert, et al., 1986). 


G. SUMMARY 

The theoretical understanding of the physical basis for scatterometry 
requires working knowledge in the interrelated theories of electromagnetic 
scattering, Ocean wave dynamics and air-sea interactions. The scattering of 
electromagnetic radiation from evolving, random surfaces leads directly to 
the most poorly understood, but currently, most actively researched problems 
in ocean wave dynamics: wave generation by the wind, nonlinear 
interactions, wavebreaking and long wave-short wave modulations. The 
presence of the moving waves affects the air flow above them and also 
precludes direct measure of the stress on the surface. Therefore, the 


connection of the winds to the waves is accomplished indirectly in air-sea 
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Figure 20a. 25° Incidence Angle. Comparison of Predictions of Plant Model 
with K, Band Data Given by Schroeder et al. [1982] and the SASS 1 Model 
Function. Data Points are those Seasat Cross Sections which Could Be 


. Associated with Good Surface Measurements of Wind. Solid Lines Represent 


the SASS 1 Model Function, While Dashed Lines Represent Predictions of the 
Plant Model (from Plant, 1986). 
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Figure 20b. 40° Incidence Angle. Comparison of Predictions of Plant Model 
with K, Band Data Given by Schroeder et al. [1982] and the SASS 1 Model 
Function. Data Points are those Seasat Cross Sections which Could Be 
Associated with Good Surface Measurements of Wind. Solid Lines Represent 
the SASS 1 Model Function, While Dashed Lines Represent Predictions of the 
Plant Model (from Plant, 1986). 
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interaction theory by the parameterization of the surface roughness in terms 
of a drag coefficient or roughness length. These quantities are themselves 
functions of the wave state, wind speed and atmospheric stratification. Since 
the wind stress is of central importance to ocean circulation studies, wave 
dynamics and remote sensing an entire body of literature which documents 
attempts to define appropriate parameterizations has sprung up in the past 
decade. 

The authors of the model functions described in this chapter attempt 
to consolidate into a scatterometer model function current information 
regarding the intertwining theories. In the process, certain environmental 
parameters are identified which may affect the radar cross-section. It now 
remains to be seen if the proposed dependencies are supported by 


observations. 


SN 


TABLE III. MATRIX OF MODEL FUNCTION CHARACTERISTICS 


MODEL 
(Author) 


SASS] 
SASS 
Team 
(1980) 


55 
VVentz, 
et al. 
(1984, 
1986) 


Durden 
& 
Vesecky 
. (1966) 


SCATTERING 
THEORY 


“Power 
law" /cosine 
statistical fit 
of aircraft 
scatteromet 
er and in- 
situ mea- 
surements 
(RADSCAT 
, JASIN 
GOASEX) 


“Power 
law” /cosine 
Statistical fit 
to SEASAT 
data 


Quasispecul 
ar Tvvo- 
Scale Model 
Scattered 
Fields com- 
puted nu- 
merically 


INTRINSIC 
RADAR 
PARAMETER 


07-70” 
İncidence 
Ku-band 
Oniy 
Symmetric 
in azimuth, 
minima at 
70777 


0°-70° inci- 
dence Ku- 
band only 


Symmetrical 


in azimuth 


07-70” inci- 
dence 

All fre- 
quency 
Symmetric 
in azimuth 
Requires 

e 
Salinity for 
dielectric 
constant 


WAVE THEORY 


None 


None 


v(k)zS(k)dbio) 

S(k): Pierson- 
Moskowitz spectrum 
low wavenumbers; 
stress dependence at 
high wavenumbers 
Q0): cos29 depen- 
dence, effective only 
near Bragg wavenum- 
ber 

Linear modulation of 
high wavenumbers 
using Cox & Munk 
(1954) clean sfc slopes 


GE: 


SURFACE LAYER 
MODEL 


Kondo (1975) 
model enters 
implicitly in 
the conver- 
sion of in-situ 
windspeed to 
Ux (19.5) 

See Boggs 
(1981) 


Assumes 
global wind 
speed has a 
Ravleigh dis- 
tribution, 
wind direc- 
tion is uni- 
formly dis- 
tributed 


Large & Pond 
drag coeffi- 
cient 
Cardone zo 
corrects for 
stability 

çiz/ L) needs 
Tam Isea 


COMMENTS 


See problems 
noted in text, 
3آ‎ 


Produces er- 
roneous re- 
sults for 
U<3.5 mps 
due tO Zo. 
Extension of 
basic wave 
mode! by 
Fung & Lee 
(1982) 


TABLE III. MATRIX OF MODEL FUNCTION CHARACTERISTICS (CONT. 


Plant Two-Scale 

(1986) Bragg-scat- 
tering do- 
main only 


Donelan Quasispec- 
& ular Two- 
Pierson Scale 


(1987) 


20°-60° 
Incidence 
All fre- 
quency 
Symmetric 
in azimuth 
Constant 
dielectric 
constant = 81 


0°-70 
Incidence 
All fre- 
quency 
Symmetric 
in Azimuth 
Look up 
table for di- 
electric con- 
stant, re- 
quires lepa, 
Salinity. 


v/(k) derived from bal- 
ance of terms in wave 
action equation; mag- 
nitude and directional 
dependence arises 
from cosine form of 
wave growth parame- 


terization. 


Uses Cox & Munk 
(1954) slick covered 
slope equation to es- 
timate long-wave 


slope 


Wk)=Sk) PM) 


S(k): Donelan (1985) 
spectrum for low 
wavenumbers; bal- 
ance of wind input 
with wave breaking 
and viscous dissi- 


pation at high 
wavenumbers 


p(o): Seck? (wind- 
wave angle) Linear 
modulation of high 
wavenumbers using 
slopes determined 
from integration of 
empirical/theoretical 


spectrum 


23 


Neutral 
Conditions, 
Large & Pond 
drag coeffi- 
cient 


Neutral 
conditions 
Large and 
Pond or 
Donelan drag 
coefficient 
uses Neutral 


wind at ; 


Bragg 
vvavelength 


Predicts o 
has a long 
wave slope 
dependence 


Predicts water 
temperature, 
wind speed, 
dependent 
threshold 
wind speed 


for G°. 


III. THE LAKE ONTARIO WAVES '87 EXPERIMENT 


A. INTRODUCTION 

The hypothesis of this study is that the radar cross-section varies with the 
wind vector and other environmental parameters; in particular, atmospheric 
stratification, long wave slope and water temperature. Experiments were 
conducted in 1986 and 1987 from an instrumented, moored tower in Lake 
Ontario in order to obtain the data necessary to test the hypothesis above. The 
fundamental measurements included wind speed and direction, air-water 
temperatures, relative humidity, radar cross-section and wave spectra. The 
purpose of this chapter is to describe the experiment and instruments by 


which these measurements were obtained. 


B. THE EXPERIMENTAL SITE 

The data analyzed in this study were acquired between 11 November and 
12 December 1987 from an instrumented, bottom-mounted tower at the 
western end of Lake Ontario, near Hamilton, Ontario (Figure 21).! The tower, 
constructed by the Canada Centre for Inland Waters, has been in nearly 
continual use since 1976 for wave and air-water interaction research. 

As described by Donelan et al. (1985, p. 515), the tower is positioned 1.1 km 


offshore at a location with the following characteristics : 


‘Selected data from the WAVES 86 experiment are also used in some of 
the analyses which require a longer averaging interval or to augment the 
range of environmental conditions. 
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Elevation ¢ 100m 





Figure 21. Map Showing the Location of the Research Tower in Lake Ontario 
and the Shore-Normal Profile in the Vicinity of the Tower 
(Donelan, et al., 1985) 


...the bottom slopes relatively rapidly (about 11 m/km) from the shore to 
the location of the tower at 12 m depth; whereas at and beyond the tower, 
the bottom slope is gentle (1.5 m / km). In addition, the shoreline is very 
straight and the bottom contours are parallel for 3 km in both directions. 
The annual variation in water level is less than 0.5 m; tides, seiches and 
wind set-up change the water level by, at most 0.1 m; there are no 
significant tidal or seiche currents, and other less organized currents are 
typically less than 10 cm/s. 


The structure is of a bilevel design (Figure 22) with an upper deck (100 m^) 
constructed of open grating to minimize flow distortion, and a lower 
walkway four meters above the water surface. At the center of the upper deck 
is a mast on which the meteorological sensors are mounted; at the foot of the 
mast is the data acquisition system (DAS) housing. Power to the tower (6 KW 
at 600 V) is supplied from an onshore source via an underwater cable. 

During the WAVES 87 experiment environmental conditions at the 
tower exhibited the ranges summarized in Table IV. Comparison of these 
values with climatology shows that the environmental conditions were 
typical of late autumn in Lake Ontario (Great Lakes Climatological Atlas, 
1986). To give an overall perspective to the experiment, the average values 
of air and water temperatures, air/water temperature difference, relative 
humidity, wind speed and wind direction for each run are plotted as time 
series in Figure 22. An unusually warm weather period which occurred 
during the week of 11 November 1987 (Days 5-10, Figure 23a) greatly expanded 
the range of air-water temperature differences beyond the climatological 
norm. The air - water temperature difference was mostly negative during the 
next two weeks of the experiment and the final week had several days of 
relatively warm air temperatures. The relative humidity was quite variable, 
ranging from 40 to 100%; precipitation is noted for those days which show a 
constant 90% value (see Figure 23b, Day 14). Only moderate wind speeds less 


than 15m/s were measured in the part of the data set currently available to 
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Figure 22. Photograph of the CCIW Tower Configured for WAVES'87. 


— 


TABLEIV. SUMMARY OF ENVIRONMENTAL CONDITIONS DURING 
THE WAVES'87 EXPERIMENT 


Air Temperature “lo ٦ 
Water temiperature” EE 900706 
Air-water temperature difference ....-13.8 to 9.6° C 
Relative humidity 40 to 100% 
Wind speed 


Wind direction 
Significant Wave Height 


Wind speed <5 m s ño Zn >7.5 m/s 
No. of data runs 39 40 47 


Fetch Short Long 
No. oí data runs ` 100 26 


(as determined from wind direction: long fetch cases occur when 
the wind direction is between 57 and 83 deg T) 


Stability Unstable Neutral Stable 
AT-0.5 |AT 1 <0.5 
No. of dataruns 56 8 
(here, an indication of stability is given by the air-water temperature 
difference) 
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Figure 23. Average Environmental Conditions During VVAVES 87 vs. Time 
(a) Air-VVater Temperature, (b) Relative Humidity 
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Figure 23. (Cont.) (c) VVind Speed and Air-VVater Temperature Difference, 


and (d) VVind Direction. 
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the author (Figure 23c)*. A feature to note is that there were no winds from 


the southeast quadrant of the compass (90° to 180°). 


C MEASUREMENTS AND INSTRUMENTATION 
1. Normalized Radar Cross-Section from the Microwave 

Scatterometers 

The intrinsic and — characteristics of the NRL scatterometers 
are described in Table V. There were two dual-pol, continuous wave, 
coherent microwave radars used in the experiment, one operating at 14 GHz 
frequency (Ku-band) and the other at 5 GHz (C-band). The received signal is 
beat down by a local oscillator to a 60 MHz intermediate frequency (IF) offset. 
This signal is then run into an IF amplifying circuit which passes the full 
bandwidth and beats down the signal to zero audio frequency (baseband). The 
final filter on the system passes frequencies between 3 dB power points at 1 Hz 
and 1.7 Khz. The filter removes dc antenna feedthrough at the low 
frequencies and noise at the high frequency end while allowing passage of the 
full Doppler bandwidth in-between the half-power points. Because there was 
a limit of 48 data communication channels available to telemeter data from 
the tower, only the amplitude modulated signal from the scatterometers was 
retained. This signal is proportional to the fields received by the antennas 
ES Tice Over 0.6 msec due {to the filtering. Thus, upon calibration and 


squaring, the signal provides the radar cross-sections as desired, but it does 


ду 


? A storm event occurred during the last vveek of the experiment for 
which the wind speeds exceeded 15m/s; however, the radar data from these 
last few runs have not yet been processed. 


ILON 


not allovv discrimination of the Bragg scatterers as can be obtained through 
analysis of Doppler frequency shifts in the frequency modulated signal. 

The Ku-band scatterometer was mounted on an antenna rotator and 
then extended on a boom 7.6 m from the tower and 7.6m above the water 
surface (Figure 24). The rotator varies the azimuth angle of the radar in a 
windshield wiper fashion, omitting those angles at which the beam intersects 
the tower. One 300 degree sweep is made every ten minutes by stepping five 
degrees in azimuth approximately every ten seconds. The incidence angle of 
the radar was also varied in the sequence summarized in Table V in order to 
cover the entire range of reflection angles from specular to grazing. The 
C-band radar was mounted on the tower railing 8.2m above the water surface 
at 45° incidence angle and 250° T (from true North) (Figure 25)3. 

The two microwave systems are illustrated schematically in Figure 
26. The Ku-band radar consisted of a dual-pol, 15 cm circular horn transmitter 
and two rectangular horn receiving antennas, one for each polarization. The 
longer wavelength C-band antenna requires a larger antenna, hence this 
system was composed of a 61 cm dual-pol parabolic reflector, and two 


4 


rectangular horn receiving antennas’. A 0.8 cm lucite septum was inserted at 


3 The original intent was to have the 5GHz radar on the rotator 
alongside the 14GHz radar. Due to weight limitations on the rotator, the ideal 
arrangement of two rotating scatterometers was abandoned and the 5GHz 
radar was fixed on the railing. 


4 The rectangular horns were used as receivers instead of a second 
parabolic dish in order to reduce the physical size of the system and to 
eliminate cross-talk between the antennas. This was an important 
consideration for the tower-based radars since they were not continuously 
monitored throughout the experiment. 
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Ku-Band (14 GHz) Radar System Installed on the Boom 


TABLE V. CHARACTERISTICS OF THE NRL RADARS 


FREQ MEASURED WAVE BRAGG WAVE H or VPOL DYNAMIC 
GHz VARIABLE LENGTH 20 40 60 80 RANGE 


Ku-band NRCS 2.14 36 1 dual 
(14.0) all cm 


C-band NRCS 6.00 4.2 cm at 45 deg dual 
(5.0) 


AZIMUTH ANGLE VARIATION 


Ku-band scatterometer sweeps through approx 300 deg every ten minutes, 
from 95 through 0 to 30 deg T. 





C-band scatterometer is fixed at 250 deg T. 


INCIDENCE ANGLE VARIATIONS 








Recording Period Angle Comment 


5 minutes 0 Specular reflections 
10 minutes 10 "Calibration" angle 
no wind speed dependence 
Ku-band radar 20 minutes 20 Spec/diffuse (soft) boundary 
40 Diffuse (Bragg) scatter 
60 Diff/grazing (soft) boundary 
80 Grazing angle reflections 
(system noise evaluation) 


— — — — — — — — — — 


95 min continuous run 


C-band radar Fixed on tower railing at 45. 
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Figure 25. C-Band (5 GHz) Radar System Mounted on the Tower Railing 
Facing Southwest 
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8/10 cm 
lucite septum 


fiberglass 
Hadome 
VPOL Dual-Pol HPOL HPOL Dual-Pol VPOL 
Receiver Transmitter Receiver Receiver Transmitter Receiver 
Rectangular Circular Rectangular| Rectangular Parabolic Rectangular 
Horn Horn Horn Horn Reflector Horn 
y y 
Coordinates: b Coordinates: | Bə 
X X 
Ф, 12.62 9.17 ФИЛМИ: GER 
Фу2 2 Фу: 6.6” 
b, 10.62 7.7” b, 21° 6.6? 
Ф,, O22 Ф,. 5-37 
Ф, SE 10.6” İ o, 5. 21 
d. 6.27 d. 5.0? 
Ф, 9.17 12517 7.5? 21° 
b. 7.39 d. 7.09 
11.0 GHZ Ku-Band System 5.0 GHz C-Band System 


Figure 26. Configuration and Half-Power Full Beam-Widths of the 
Microwave System Antennas. Py is the Vertical Beamwidth and ®, is the 
Horizontal Beamwidth 
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the midpoint of the h-pol receiving antenna to match its antenna pattern to 
that of the v-pol receiving antenna. 

The mainbeam and first sidelobes of the C-band transmitting 
antenna and the h-pol receiving antenna are shovvn in Figure 27. The half 
power beamwidth (HPBW) of the transmitting antenna (width of beam 3 dB 
below peak) is 5.5 while the receiving antenna has a considerably wider 
beamwidth of 21. The wider pattern is used in conjunction with the narrow 
transmitting pattern since it allows for some alignment error when 
attempting to focus the individual beams on the same illumination area. 

The first sidelobes of the transmitting antenna are about 15 dB below 
the peak and 9° off the main beam direction. The receiving antenna sidelobes 
are approximately 19 dB below the main peak and 56 off the mainbeam. 
These two patterns are combined into a (h-pol) two-way pattern using the 
principle of multiplication of antenna patterns to yield the solid line shown 
in the figure. Here it is seen that the disadvantage of the wider receiving 
pattern is that it does not suppress the sidelobes as much as would a narrow 
main beam with sidelobes matched in position to those of the transmitting 
antenna. In the two-way pattern the HPBW is 5.5" and the first sidelobes are 
9” off and approximately 17 dB below the main lobe. This implies that the 
radar cross-section must vary by 17 dB in only 9 in order for the sidelobes to 
significantly contaminate the total cross-section measurement. The v-pol 
two-way pattern is the same as the h-pol pattern since the two receiving 
antennas have the same one-way pattern. A similar analysis of the Ku-band 
system gives a vertical by horizontal HPBVV of 7.3” by 6.2” for the mainlobe, 


and first sidelobes 26 dB belovv and 9” off the main peak. 
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Figure 27. Antenna Patterns for C-Band (5GHz) Radar System 
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The elliptical surface area, A, illuminated by the radars is calculated 


from the tvvo-vvay antenna beamvvidths as 





7 aD, = 


cos 6 


Jee.) 5). (107) 


cos” 9 


where a is the antenna height, Фу2, Фх2, are the vertical and horizontal two- 
way beamwidths at half power, and @ is the incidence angle (Plant, 1990). 
Since the one-way beamwidths at half power for the transmitting and 
receiving antennas are unequal, the two-way beamwidth is calculated as the 


sum of the inverses of the individual one-way beamwidths (Silver, 1949) 


q ə ہے‎ (108) 


where DIR and OT are the one-way horizontal receiving and transmitting 
beamwidths of the receiving and transmitting antennas. A similar 
expression is used to compute the vertical two-way beamwidth, ®y2. Using 
the values given in Figure 25 for the one-way beamwidths and given that a is 
equal to 7.6m, the illumination areas (1, by ly) for the different incidence angle 
settings of the Ku-band radar are .82 by .97m* (at 0 incidence angle), .84 by 
1.0m? (10°), .88 by 1.1m” (207), 1.1 by 1.6m” или, by 3.9m? (60°), 4.8 by 
32.2m7 (80^). The C-band area of illumination is about .75 by 1.3m* for 0 and a 
equal to 45^ and 8.2m. 

It is important that the illuminated area be in the far-field of the 
transmitting antenna. That is, far enough away from the source that the 
components of the total field due to the source itself can be neglected and the 
spherical wavefronts of the impinging radiation are essentially planar, 


thereby allowing uniform illumination of the area. For a single antenna, the 
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limiting distance which defines the far-field to within a x/8 error in phase is 
determined from 


2D? 
Rr = X (109) 


where D is the largest linear dimension of the antenna and A is the 
electromagnetic wavelength. The one-way beamwidth di can De 
approximated by A/D; therefore, by substitution, 


22 


Rs 
- (110) 


In the present case of separate transmitting and receiving antennas, the two- 


way beamwidth must be used. From equation (110), with equal receiving 


2 1 2 
раНегп5 Ф, = ən and 


Rep = — (111) 


Finally, with unequal vertical and horizontal beamwidths, the far-field 


distance is 





ob SE (112) 


In the case of the Ku-band system, the far-field distance is 1.5m for A 
equal to 2.14cm. The range of the Ku-band radar varied from 7.6m at nadir to 
44m at 80 incidence: therefore, all measurements were made in the far-field. 
similarly for the C-band antenna and à equal to 6cm, the far-field distance is 


5.6m and all measurements were made in the far-field. 
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2. VVind Speed and Direction from the Bivane Anemometer 

A Gill-type bivane anemometer, mounted on the tovver mast at 
11.5 m (Figure 28), was used to measure total wind speed and simultaneous 
azimuth and elevation angles of the wind (R.M. Young Co., 1974). The wind 
speed sensor is a four blade polystyrene propellor which drives a tachometer 
generator mounted in the counterweight section of the vane assembly, whose 
output is a d.c. voltage directly proportional to wind speed. The wind speed 
range measureable with this propellor is from a threshold speed of about 
0.2m/s to a maximum of 31m/s. The propellor was designed to provide 
approximately one revolution for each foot of wind for all wind speeds. In 
fact, this criterion is met for all wind speeds above 1.2m/s, which means that 
between the threshold speed of 0.2m/s and 1.2m/s, the propellor is "spinning 
up" and the wind speeds are not reliable. 

The fin portion of the vane assembly is also molded of polystyrene. 
With the passage of a wind gust (threshold speed 0.2m/s), the light vane 
moves horizontally and vertically to maintain equilibrium of wind pressure 
on the vertical and horizontal vane surfaces. Changes in azimuth are 
transmitted through the main vertical shaft to the shaft of a linear conductive 
plastic potentiometer by a set of gears with a one-to-one ratio. Changes in 
elevation angle are transmitted by a pulley to a small bead chain passing 
down the hollow vertical shaft to a pulley mounted on the elevation 
potentiometer. A counterweight on the lower pulley eliminates backlash in 
the system. Elevation angles of +50 are measureable. With a constant 
voltage applied to the elements of the potentiometers, the voltage outputs are 


directly proportional to the position of the vane. 
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Figure 28. The Bivane Anemometer on the Research Tower Mast 11.5m 
above the Water Surface 
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The azimuth potentiometer has a 5% open section which equals Ee 
to 9” rotation ("dead angle"). Therefore, the full scale calibration of the 
azimuth signal is only 352° instead of the full 360°. In the field, the "dead 
angle" is aligned with true South. When the wind is from the South, the 
voltage will oscillate between the minimum and maximum values (OV, 160°; 
5V1527). The average value (2.5V,about 156°) will have a correspondingly 
large standard deviation which can be monitored and corrected for in data 
processing. 

For optimum performance, the propellor and vane assemblies have 
carefully matched dynamic response characteristics. The distance constant of 
the propellor (the wind passage for 63% recovery for a step change in speed) is 
“1.0m. The delay distance of the vane (the wind passage for 50% recovery from 
a step change in direction) is 1.0m and its damped natural wavelength 
(distance in which directional oscillations are reduced by 1/e) is 5.8m. 

3. Wave Spectra from a Six-Element Wave Staff Array 

The six wave staffs deployed in this experiment were arranged at the 
apices and center of a pentagon. The distance between the center staff and the 
outer wave staffs is 0.25m. The azimuthal resolution of this symmetrical 
array depends on processing parameters, here taken to be 5°; the range of 
wave periods resolved is 0.6 to 8.0s. 

Following Der (1976), the sensors are capacitance transduction 
devices which produce output voltage signals proportional to surface 
elevation. The sensors in WAVES87 were comprised of 6m téflon-sheathed 
wire probes of 4.8mm outer diameter and associated signal conditioning 


electronics. Assuming uniformity in the cross-sectional area of the teflon 
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sheath, in its surface and its chemical composition, then the capacitance of the 
proble vvill be linearly related to the length submerged. The mafor source of 
error is stray capacitance-to-ground which is maximized when the probe is 
fully exposed; i.e, is most nearly like a vertical antenna perpendicular to a 
ground plane. Wetting errors due to water films and algae adhering to the 
teflon surface, account for most of the hysteresis exhibited by these sensors. 
According to Tsanis and Donelan (1988), the wave staffs had a very stable 
linear calibration during the WAVES87 experimental period and were 
cleaned almost weekly from any impurities and algae that had accumulated 
on them. 
4. Air Temperature and Relative Humidity from the Fast Response 

Sensors 

The air temperature thermistor at 11.5m had a time constant of 0.2s 
and was shielded to minimize direct solar warming. The relative humidity 
sensor was mounted beneath a plate at the top of the mast on which the base 
of the anemometer rests. The sensing element is a thin film capacitor 
composed of an upper and lower electrode with an organic polymer dielectric 
about One micron thick (WeatherMeasure Corp., 1975). Water vapor is 
absorbed into the polymer after the vapors pass through the upper metal 
electrode thereby changing its capacitance. Comparison of the capacitor 
output to a reference signal provides a DC voltage proportional to humidity. 
Since the upper electrode must be as transparent to water vapor as possible, 
the thin electrode results in quick response to changes in humidity. The 


response time of the fast response sensor element is less than one second at 
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68F to 90% of the final relative humidity value. The accuracy of the sensor is 
better than +3%. 
5. Surface Water Temperature from an in-situ Thermometer 
The water surface temperature sensor is contained in a modified 
"bucket" which is fixed just beneath the water surface against one of the tower 
legs. The bucket has an annular space between its lid and body and a small 
hole in its base. As a wave washes over the bucket it is periodically filled and 
then slowly drained through the hole. In this way, the sensor is always wet 
and measures a temperature corresponding to the upper 10 cm of water 
(Tsanis and Donelan, 1987). 
6. Ancillary Measurements 
Synoptic charts of surface pressure, temperature and wind velocity 
were obtained from Fleet Numerical Oceanography Center on a daily basis 
from 9 October to 11 December to provide a larger context for the local 
weather at the experimental site and to facilitate selection of test cases based 
on similar weather conditions. 
Aside from the radars, information regarding the instruments 
installed on the tower was obtained primarily from Tsanis and Donelan, 1987. 
Calibration of the sensors was performed in the laboratory before and after 


field exposure. 


D. DATA ACQUISITION 

Two alternate modes of data acquisition were employed. In the high 
frequency (continuous) mode, all sensors were sampled and recerded at 20 Hz 
for 95 minutes and in the low frequencv (averaging) mode, samples were 


made at 20Hz and averaged every ten minutes. The continuous mode was 
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used mostly during regular vvork hours only, since this mode required an 
attendant. VVhen significant vveather events occurred, vvatches vvere posted 
and the high frequency mode continued through the night. 

The 20 Hz, 48 channel data vvere digitized on the tovver and transmitted to 
a trailer onshore where they were automatically received by an IBM-AT 
microcomputer. They were then stored on the IBM's 30 MB hard disk and 
plotted in analog form on chart recorders for real-time monitoring. Each 95 
minute run required 7 Mb so that the hard disk is full after four runs. At this 
point the data were copied onto 40 MB magnetic tape, the hard disk erased 
and the logging sequence reinitiated by the operator. 

The ten minute averaged data followed the same procedure as above 


except that they were written directly to floppy disks. 


E. SUMMARY 

The WAVESS7 experiment was performed from a tower in Lake Ontario 
during October to December, 1987. One aspect of the experiment was to 
investigate environmental effects on the NRCS measured by radar 
scatterometers. In support of this investigation, the tower was outfitted with 
a suite of instruments which provided the meteorological and limnological 
parameters thought to influence the 9۰۰1777 
experimental site and the weather conditions which prevailed during the 
experiment. The intrinsic and operating characteristics of the radars, bivane 
anemometer, temperature and humidity sensors, and wave staffs were 
discussed with the intent of identifying the limitations of each of these 


devices. 
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IV. DATA PROCESSING METHODS AND RESULTS 


The 14 data channels pertinent to this study were copied from the raw 
data magnetic tapes to backup tapes at CCIW and subsequently were written to 
the NPS IBM3033 mainframe computer for processing. At this stage, the data 
consisted of voltages to be converted to geophysical units via the laboratory 
calibration equations. Next, the raw environmental data needed averaging in 
time intervals which were sensible for a particular geophysical parameter and 
the resulting averages assigned a time tag for subsequent collocation of 
multiple parameters. Since the rotating Ku-band radar moved in azimuth 
every ten seconds, while the C-band radar was fixed in one position, different 
averaging schemes were applied to the individual radars and their associated 


environmental parameters. 


A. THE KU-BAND SCATTEROMETER AND ASSOCIATED 
ENVIRONMENTAL PARAMETERS 


An example of the voltage outputs of the Ku-band radar channels is 
shown in analog form in Figure 29, which illustrates the first 7.5 minutes of 
Run 87056 separated into three 2.5 minute segments of elevation (incidence 
angle) (EL), azimuth (AZ), Ku-band horizontal polarization (KuH) and Ku- 
band vertical polariztion (KuV) voltages. For approximately 40s of the first 
minute the radar samples the calibration circuit while at its "home" position 
of 60° (incidence) and 250' T North (azimuth). In the next 50s, the radar 
completes its movement to the initial starting point of O° incidence and 95°T 
azimuth. The calibration circuit is switched for the Ku-hpol and Ku-vpol 


channels and typical random signals are observed. 
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Figure 29. Analog Record of the First 7.5 Minutes of Data Run 56 
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The azimuth signal is seen to decrease in small increments, until about 
four minutes have elapsed. At this time, the incidence angle is increased to 
10° and the radar is rotated back to 95° T. This procedure is repeated for the 
remaining incidence angles (recall Table V) except that the radar does not 
return to 95° T each time but simply steps backward from the end of the 
previous aZimuth sweep. 

1. Incidence and Azimuth Angle Position 

The incidence and azimuth angle positions of the Ku-band radar 
were calibrated before and after the experiment using the most precisely 
known physical positions of the radar as anchor points. That is, elevation 
was measured with respect to true vertical as visualized with a Sperry digital 
Bevel, as observed on the antenna rotator azimuth/elevation meter and as 
measured by the rotator output voltage. The azimuth angle was calibrated 
with respect to the known tower orientation of 250° T and the points + 90° of 
250°T. With this calibration method, it is believed that the incidence angle is 
absolutely determined to within 1° and the azimuth angle is absolutely 
determined to within 5” 

The radar incidence and azimuth angles are determined by the 
voltage levels recorded on the appropriate data channels. The voltage values 
at each elevation setting are known, virtually constant and measured for no 
less than four minutes. Therefore, the radar incidence angle movements in 
the analog data record are easily discriminated as changes from the known 
values (within a small tolerance). Spikes in the incidence angle voltage 


measurements, which are misinterpreted as elevation movements, are 
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compensated for by requiring that the apparent movement span at least five 
consecutive data polnts. 

The azimuth angle determination proved much more difficult. 
Tolerance checks based on explicit voltage values are not possible for the 
azimuth angle movements for vvhich the voltage increment betvveen 
azimuth steps is only 71mV and the measurement period is less than ten 
seconds. To discriminate between adjacent steps, a five point running mean 
of the azimuth voltage is computed. A step change is considered to have 
occurred when the difference between two consecutive averages exceeds 
7mV. All of Ku-hpol and Ku-vpol voltage values up to the step change are 
binned together and their statistics (mean, variance, skewness and kurtosis) 
are computed in preparation for the cross-section computations. 

The azimuth direction was computed during the preliminary data 
processing stage using the average azimuth voltage and the calibration 
equation as measured in the field. However, small amounts of noise in the 
voltage measurement occasionally caused the wrong direction to be produced 
by the calibration equation (usually the direction was shifted backwards by 
one five degree bin and the data in the previous bin overwritten). The error 
induced by this bin shifting was cumulative: the position at the end of a ten 
minute azimuth sweep was as much as 15° from the true value. Visual 
examination of the analog plots revealed other anomalies: the total number 
of five degree bins per azimuth sweep was 59 instead of 60, bin durations were 
occasionally shorter or longer than the ideal 10s length, voltage spikes in the 
bin produced anomalous movements, and there was no clear separation 


between the first and second pass at some incidence angles. For these reasons, 
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a less automatic but more reliable, procedure vvas adopted for assigning the 
correct azimuth direction. 

In this procedure, a counter equal to the correct bin number is 
initialized at the start of each azimuth sweep. The counter is then 
incremented (decremented) at each azimuth step of a forward (backward) 
sweep. In this way each real azimuth step is assigned a unique direction. 
Comparison of the bin numbers chosen this way to those computed by the 
actual calibration equation reveal no difference until a "noisy" average 
voltage measurement occurs. Only at this time does the calibration equation 
compute the wrong backwards shifted value, while the counter correctly 
moves to the next step. 

The occasions on which the forward sweep goes beyond the fifty- 
ninth bin, i.e., the lowest voltage is .47xx instead of .5xxx, are accounted for by 
setting the initial value of the bin counter one vnit higher than normal on 
the backward sweep. Voltage spikes are accounted for bv requiring a bin to 
consist of a certain minimum number of points (usually 100 of the possible 
180). Bins which are truly shorter or longer than the normal 10s are still 
correctly differentiated by changes in the running mean values. Finally, 
separation of the two sweeps is accomplished when necessary by splitting the 
second sweep from the first according to a time duration limit of 608s for the 
first sweep. The absence of a systematic bias between the forward and 
backward sweeps was established by plotting the azimuthally averaged 
cross-sections from the first and second passes against each other for several 


“ 
data runs (Figure 30). 
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Figure 30. Azimuthly Averaged NRCS of First Radar Rotation Plotted against 
Second Rotation. Upper: 40° v-pol. Lower: 40° h-pol. 


To summarize the Ku-band radar position processing, the full length 
analog plots of azimuth and incidence angle voltages are visually examined 
to identify anomalous movements or potential problems. These problems 
alvvays fall into the above mentioned categories, hence standard corrections 
are vvritten into the data processing softvvare. The corrections are applied by 
setting the appropriate flags in the data statements of the processing program. 

1. Internal Calibration 

Errors in the cross-section values due to performance characteristics 
of the radar components are eliminated through internal calibration of the 
radar at the beginning of each full 95 minute data run. A portion of the 
transmitted signal is sampled after its passage through the internal 
components of the radar system (amplifier, mixers, and oscillator). The 
measured calibration voltage is recorded on the Ku-hpol channel. In Figure 
29, the calibration period is 51.65 (nominal: 50s) and the calibration voltage is 
0.299 volts (nominal: 0.30 volts). Over the course of the experiment, the 
calibration voltage increases at a rate of approximately 2mV/day (Figure 31a) 
but it shows no dependence on air temperature (Figure 31b) as would be 
expected if the efficiency of these internal components were affected by air 


perature. 
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Figure 31a. Internal Calibration Voltage of Ku-Band Radar as a Function of 
Time from Start of Experiment (November 11, 1987) 
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Figure 31b. Internal Calibration Voltage of Ku-Band Radar as a Function of 
Air Temperature 


Drift of the calibration voltage is compensated for by using the actual 
measured calibration voltage in the NRCS computation except for the runs 
which do not begin at the radar home position. For these incomplete runs, 
the laboratory calibration voltage of 0.30 volts was used.! 

3. External Calibration 

The radar system was calibrated at NRL before and after the 
experiment using two metallic spheres of known radar cross-section to 
establish the absolute calibration of the measured cross-section. The spheres 
were suspended from strings and oscillated to discriminate the radar return 
from background reflections by the Doppler shift of the returns from the 
moving spheres. A 1 KHz filter was used to eliminate returns from 
nonmoving objects (noise). 

The radar equation which expresses the power received at the 
antenna in terms of the antenna properties and the target cross-section is 
written as 


m5 o 


(113) ا نے 
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where p is the ratio of the power received (Pj), to the power transmitted (Py, 
A is the microwavelength, G is the antenna gain, R is the range to the target 


and O is the target cross-section. The cross-sections measured in the field are 


lThe preferred value would have been either the value from the next 
closest data run or the regressed value derived from the time series of the 
calibration voltage. Since the data were not processed in sequence and since 
the calibration voltages were required at the start of the data processing, the 
laboratory value was adopted as the default. 
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absolutely calibrated by forming the ratio of pş of the sea-surface to pç for the 
calibration sphere and deducing a calibration constant. Since the antenna 


parameters are identical for each p, the ratio becomes 


4 
p, _ S, F. (114) 


so that the sea-surface cross-section 1s 


4 
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The area-normalized radar cross-section 1s then 
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With the illuminated area A defined as in equation (107), and R  a/cosö, 


then 
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where Co = c'/cos0 p, accounts for the system and calibration constants and 
incidence angle. This constant was found to be 2.31 / cos@ for the Ku-band 
radar. 
4. System Noise Evaluation 
A continuous wave radar, such as those used in this experiment, is 
always transmitting a signal. Hence, system noise defined as the voltage 


measured by the radar in the absence of any target, is estimated by 
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minimizing the return signal. Typically, this means pointing the antenna 
into free space away from any potentially reflecting targets. In this 
experiment, when under very low wind conditions, the 80° incidence angle 
was considered equivalent to looking into free space. These data were used to 
compute an "equivalent noise variance," N, which was then subtracted from 
the signa] variance in the computation of the absolute cross-section. 

The basis for the equivalent noise variance is that if the computed 
cross-section is actually representative of the system noise only, then Pg is a 


constant and 


(119) 





should be nearly equal for the 60° and 80° incidence angles, where the average 
on the cross-section, denoted by angle brackets, is over all azimuth angles. In 
this noise case, the cross-sections are constant with azimuth and their 
magnitudes increase with incidence angle by virtue of the radar equation. 
From equation (117), the ratio of the average cross-section at 80° to that at 60° 
IS 


577 Gas 80) ` p. 


— = -—” (120) 
GEG CCM 0) 





In log space the ratio is a positive difference of 10log(2.87) or 4.5 dB between 
the 80° and 60° average cross-sections. 
Figure 32 illustrates the last 7.5 minutes of Run 87048 from which a noise 


computation is made. At this stage, the radar is at 80* incidefice and about 
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Figure 32. Final 3.5 Minutes of Run 48, Used for Evaluation of the Ku-Band 
Radar System Noise Level 
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two minutes into the second azimuth sweep. Return from the tower legs 1s 
observed as large sinusoidal fluctuations between times 24:36 and 24:38. 
Because of these erroneous returns, the five azimuth bins at the extreme 
ends of the sweep at 60° and 80° are eliminated from further analysis. Beyond 
24:38 the returns are clearly down to the digitization level of the instrument. 
The average wind speed for this run was 2.6 m/s. 

Using these data and the internal calibration voltage, the NRCS was 
computed and plotted as a function of azimuth (Figure 33) and incidence 
angle (Figure 34). As explained above, the h-pol NRCS at 80° is in fact a 
relatively constant value greater in magnitude than that at 60°. From 
equation (4), the noise variance, N, was computed to be 1.6E-04 volts-squared. 
Other similarly low wind speed runs were used to arrive at an average noise 
value of 1.5E-04 volts-squared. This value is then used to compute the lower 


limit of NRCS values for each incidence angle from 


- 2.31 1 sE-04) (121) 


cos 0 





mun 


which yields o°_. of -34.6 dB (0°), -34.5 dB (10°), -34.3 dB (20°), -33.4 (407), 


-31.5 dB (60°) and -26.9 dB (80°). In subsequent analyses, NRCS values lower 

than these for a particular incidence angle are eliminated from consideration. 
5. Normalized Radar Cross-Section (NRCS) 

In this experiment, the radar return voltage was measured by a 

linear detector. Therefore, the return power is computed from the mean- 
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square ac voltage variance, V; 


ac? 


or the voltage variance, o, since 
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Figure 33. NRCS-Versus Azimuth Angle of Radar for all Incidence Angles. 
Data Run 48, November 18, 2312 GMT. Upper: v-pol. Lower: h-pol 
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Figure 34. NRCS Versus Incidence Angle. Data Run 48, November 18, 
2312 GMT. Upper. v-pol; Lower: h-pol 
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Here, only the mean-square ac component of the return voltage, Vac, 1s 
considered since the voltages are Doppler shifted due to the moving vvater 
surface. The dc component due to microvvave leakage betvveen the antennae 
is removed from the returns vvith a 1Hz high pass filter 

The NRCS is then computed from the voltage values in each 
azimuth bin using the formula 

o ја - N) (123) 
cal 
where o, is the area-normalized (nondimensional) cross-section intensity 
(not dB), 6; is the voltage variance, V7, is the calibration voltage squared, N is 
the noise value, and C, is the calibration factor. 

The cross-sections were then plotted as a function of radar azimuth 
angle. Although the azimuthally averaged values from the first and second 
passes were comparable (recall Figure 30), the individual pass values at a 
particular azimuth bin often exhibited significant scatter (Figure 33). 
Therefore, the two-pass values were averaged to reduce the scatter and reveal 
more of the azimuthal dependence. 

The error bars on Figure 35 are at the 95% confidence level and are 
estimates of the random measurement error for each average assuming that 
the voltage variances, from which the cross-sections are determined, are chi- 
Square distributed. This assumption is valid provided that the measured 
voltages are normally distributed. Figure 36a and 36b illustrate histograms of 


all voltages at incidence angle 20° and both polarizations, from which it is 
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Figure 35. Two-Pass Averaged NRCS as a Function of Radar Azimuth. 
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Dotted Lines Show 95% Confidence Levels Based on Chi-Square Distribution 
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Figure 36. Histograms Showing Gaussian Distribution of Voltage 
Measurements from Ku-Band Radar. (left) 20” v-pol, (right) 20” h-pol. 
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seen that the voltages are in fact very nearly Gaussian as required. Because 
each averaged value has 360 degrees of freedom, the error bars are 
correspondingly small (+1dB).? 

The final step in the processing of the cross-sections was to smooth 
the averaged values across three azimuth bins using a 1:2:1 weighted average 
filter (Figure 37). In the analysis figures to follow, the cross-section values 
will be identified when necessary as individual pass, two-pass averaged or 
smoothed values according to their level of processing as described above. 

6. Field Calibration using the Nadir and 10 Degree Incidence Data 

After smoothing the cross-sections and plotting them with respect to 
azimuth, it was noticed that the nadir incidence angle v-pol and h-pol values 
were offset from each other. Specular scattering theory predicts no difference 
between v-pol and h-pol at nadir; therefore, this bias was taken to be a 
calibration error between the two polarizations. On average, the bias was 
about +1.5dB (v-pol larger than h-pol). Since it was not known which 


polarization was in error, the correction procedure was to subtract 


2 Under most circumstances, the first and second pass values are close 
enough that their individual error bars would overlap, implying that they 
come from the same population. However, there are occasions when the pass 
by pass values are sufficiently separated to suggest that the cross-sections are 
from different populations. These data runs should be examined on a point 
by point basis and correlated with (possible) short term changes in the wind 
feld. 
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half of the bias from the v-pol and add half to the h-pol cross-section values. 
This method produces identical cross-section values for the tvvo polarizations 
at nadir. 

The NRCS is negatively correlated vvith vvind speed at nadir and 
positively correlated at progressively higher incidence angles (recall Figure 2). 
The 10° incidence angle is the "transition" angle between the two regions and 
as such, displays virtually no dependence on windspeed and minimal 
dependence on wind direction. Therefore, this angle is used as an external 
"field" calibration point. This calibration point may be more representative 
of the true calibration since the radar is designed to detect Doppler shifts in 
the radar returns induced by the movement of the ocean surface. This 
27: is only simulated in the laboratory procedure through the use of an 
oscillating or rotating sphere. 

Initial comparisons of the measured NRCS with the model function 
predictions as a function of incidence angle showed a positive 2 to 3 dB bias 
between the measured and predicted values at all incidence angles. Before 
attributing erroneous physical meaning to these biases, a further antenna 
pattern correction was applied to the NRCS data using the 10° incidence angle 
data. 

It was presupposed that the tower radar cross-sections at 10° may be 
higher than the corresponding predictions from the model functions which 
were tuned to airplane or satellite scatterometers, since the tower radar 
antenna has a beamwidth of nearly 7° as compared to 1° or less Yor the higher 


altitude scatterometers. That is, while the tower radar has a boresight 
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direction of 10°, the average power may be coming from an angle less than 
this, thus producing higher cross-section values. 

To estimate the offset which would arise from the wider beamwidth, 
SASS1 and SASS2 predictions at 10° (averaged over all azimuth and 
windspeeds between 3 and 11m/s) were "smeared" over a 7° beamwidth by 
assuming a Gaussian antenna pattern for the tower radar. The SASS1 and 
SASS2 model functions were chosen for this calibration since they are purely 
statistical fits of data from RADSCAT (1.5? beamwidth) and SEASAT SASS 
(0.5 ” beamvvidth). 

The half-povver points of the Gaussian curve are mapped to the radar 
antenna half povver points (of knovvn incidence angles), Figure 38. The 
standard coordinates of the Gaussian curve are then obtained as +/-0.67 or 
1.34 standard deviations; i.e., about 5? per standard deviation. The averaged 
cross-section intensity values of either the SASS1 or SASS2 model over 5° to 
15° are then weighted by the Gaussian ordinates to obtain the average value at 
10° for a 6.7° beamwidth. 

For SASS1, the averaged value was 7.73dB or 0.73dB higher than that 
obtained from the 1.5° HPBW of the aircraft RADSCAT. Similarly, the 
averaged value was 7.44dB for the SASS2 model was about 1.1dB higher. 
Therefore, about 1dB of the bias in the tower radar data can be attributed to 
the wider antenna beamwidth. 

6. Environmental Parameters 
Since the air/water temperatures, relative humidity, surface 


elevation, wind speed and wind direction are generally nearly constant 
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Figure 38. Antenna Pattern Correction to Convert Ku-Band 7.0” Beamvvidth 
to 1? Beamwidth 


during a 95 minute run, one minute running means, evaluated devery 10s are 


computed for these quantities. The first mean is calculated from the first full 


1 


minute of data after the radar reaches its initial starting position. A centered 
time tag of 30s is associated vvith the mean, hence neither the first tvvo radar 
bins (10s, 20s) nor the last two bins (5690s, 5700s) will have associated 
environmental parameters. 

It should be noted that the one minute running mean value of wind 
direction is used to compute the relative angle between the antenna and wind 
directions for each azimuth bin. The wind speed values are handled in 


several different ways as shall be described in Section C of this chapter. 


B. THE C-BAND SCATTEROMETER AND ASSOCIATED 
ENVIRONMENTAL PARAMETERS 


The data processing for the C-band radar is much simpler since it is fixed 
in one position throughout the experiment at 45°incidence and 250°T. 
azimuth. 

1. Internal Calibration 

As with the Ku-band radar, the C-band radar has an internal 
calibration cycle at the beginning of each 95 minute data run which is 
recorded on the hpol channel. The calibration period is also 50s (nominal) but 
the calibration voltage is -1.3V (nominal). Contrary to the Ku-band radar, the 
C-band calibration voltage shows no particular drift with time (Figure 39a) but 
there is clearly a dependence on air temperature (Figure 39b)3. The actual 


measured value of the calibration voltage is used to compute the NRCS. 


* [n retrospect, a multiple regression of the calibration voltage on time 


and air temperature should have been performed for both the Ku- and 
C-band radar calibration voltages. 
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Figure 39. Internal Calibration Voltage of C-Band Radar as a Function of a) 
Time from the Start of the Experiment (November 11, 1987) and b) Air 
Temperature. 
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2. External Calibration 
The C-band radar was calibrated at NRL before and after the run 
using two metallic spheres of known cross-section in the same fashion as the 
Ku-band radar. The calibration constant Cp at @ equal to 45° is 11.09. 
3. System Noise Evaluation 
No estimates of the system noise were made since the radar was 
permanently aimed at 45° to the surface and the values of the returns o... 
well above the system noise. 
4. Normalized Radar Cross-Section 


The NRCS is computed in one-minute data blocks using the formula 





o? 
o,= 1109 — (124) 
Va 
where ovy is the voltage variance within the data block. 
5. Environmental Parameters 

The environmental parameters are also computed in one-minute 
blocks to align in time with the C-band NRCS values. Judging from Figure 
40, it is clear that the radar measures nearly instantaneous changes in the 


wind speed. 


C COMPUTATION OF THE WIND STRESS, HEAT FLUX AND NEUTRAL 
WIND SPEED 


1. Justification of the Method 
The observation that the radar is representative of short-term 
changes in wind speed again confirms the basic premise of scatterometry: 


microwave radars operating at centimeter wavelengths are sensitive to wind 
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Figure 40. One-Minute C-Band NRCS and Associated Environmental 
Parameters as a Function of Time from the Start of Data Run 80, 
November 18, 0501 GMT. 
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generated capillary vvaves vvhich have grovvth and decay time scales of the 
order of seconds. This implies that the proper development of a physically- 
based radar model function must consider the local (microscale) problem of 
momentum transfer to waves which are small compared to the longest 
gravity waves present in the larger, spatially-averaged (macroscale) wave 
field. Since the momentum transfer from the wind to waves is quantified by 
the friction velocity, u*, in the growth rates of the short waves, the wind 
stress should be evaluated at time intervals comparable to the time-scale of 
these waves in microscale studies, i.e., O(seconds). 

On the other hand, the macroscale wave field exerts a drag on the 
wind which modulates the wind stress. The stress is supported by a range of 
wavenumbers whereas the radar is wavenumber selective (within a small 
band about the Bragg wavenumber). Therefore, a different time scale is 
needed to quantify sea-state effects on the wind stress. With the assumption 
of horizontal homogeneity in the surface wave field and ignoring the 
microscale, the typical approach to this problem has been to determine an 
overall roughness scale, zo, or drag coefficient, which parameterizes the wave 
drag on the wind (recall Chapter II, Section E) Wave state and stratification 
have time and space scales much larger than the microscale, therefore the 
drag coefficient has inherently longer time scales O(minutes) and larger 
spatial scales, (O( 100m)). 

The separation of the microscale effects from the macroscale is 


accomplished by writing 


: (125) 


us Olma) iiz - uo) O(sec) 


O(sec) - Cp 
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In this expression, the observed rapid changes in the cross-section with wind 
speed are accounted for by using a short time scale wind speed to compute the 
wind stress. The modulation of the stress over longer time scales is 


accounted for by computing a longer term drag coefficient. The procedure is 


to compute ux and the mean wind, <u>, from the wind components over the 


2 
long time scale and then to form the ratio of ux to <u>? for the long-term drag 


coefficient. This drag coefficient is used in (125) with the short term wind 
speed average to compute a short term stress. 

There are several ways to measure (compute) wind stress, to date the 
most reliable and accurate is the eddy covariance technique as used here 
(Geernaert, et al., 1988a; Fairall, et al., 1988). In this method, the 
instantaneous wind speed is considered to consist of a mean Uj, and 


fluctuating component, uj’ 


Uu te (126) 


where uj is the wind component in the i-th direction. The engineering 
requirements to be satisfied in the eddy correlation method are the necessity 


for reliable fast response wind, temperature and humidity sensors and a stable 


platform so that no artificial vertical motions are introduced; i.e, w must 
approach zero. 
The stress is then estimated directly by forming the covariance 


between the fluctuating horizontal and vertical components of the wind and 


applying Reynold"s rules of averaging?. The stress component in the x- 


direction is then 


u'w' -(u-u)w-W) (127) 


and, in the y-direction, 


v'w'=(v-V)\(w-W) (128) 


where the overbar denotes an average. The stress vector is then 


t= p(-w vi - vj (129) 
with magnitude 
2 27; 
It / p| = (2222) + (ow) İ (130) 
and direction, a, relative to the wind direction, given by 
Oz GEN (131) 
uw 


The velocity components are subjected to a coordinate transformation such 


that the -u'w' component is aligned with the (20 min mean) wind (along 


thesshaft of the bivane) end =~ \ sicko - right of the wind direction. 


Turbulence theory provides the guidance for the separation of the 


time scales required in (125). It has been noted that a "spectral gap" occurs at 


3 For a stationary process, all the fluctuating quantities average to zero, 
the correlation between the fluctuations and the average quantities must 
vanish, and the averaging process applied to an average must reproduce the 
same average (Panofsky and Dutton prot mO- 
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periods of about one hour betvveen the lovv frequency vvind speed maximum 
("synoptic peak”) and the high frequency wind speed maximum ("turbulent 
peak", Figure 41a (Lumley and Panofsky, 1964). From the figure it appears 
that if the averaging interval is between 0.4 and 4 hours then the longer 
period wind speed contributions will be separated from the higher frequency 
fluctuations, with the (apparently) best partition occurring between 0.7 and 
1.0hr (42 to 60min). The figure is for overland data, a similar spectral gap 
occurs over water as obtained from Lake Ontario data, Figure 41b. (Pierson, 


1983). 
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Figure 41a. Schematic Spectrum of Wind Speed Near the Ground Estimated 
from a Study of Van der Hoven (1957) 
(from Lumley and Panofsky, 1964) 


The u-w cospectrum as computed with standard spectral methods 
can also be used to determine the appropriate averaging interval. As shown 
in Figure 42, most of the turbulent energy is contributed at frequencies 


between .001 Hz (17min) and 10 Hz. Therefore, averaging intervals at least 
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17min long are required for adequate representation of the fluctuating 
components. Longer averaging times vvill include energy from the synoptic 


scale features. 


— 
“” 
— 
E 
on 
c 
eo 
c 


| hr. 


ZOmin. from meosured 


- 0 +i +2 +3 
log ,g ^ (cycles /hr) 





Figure 41b. Frequency Spectrum of Wind Speed Over Water from Pierson 
(1983) 


The entire planetary boundary layer (PBL) has a characteristic time 
scale associated with the rotation and advection times of the largest eddies 


contained within the layer. The near-neutral PBL height in mid-latitude, hpi, 


can be modelled by 


ا 
hpi = “zı (132)‏ 


where c is an empirical constant of about 0.3 and f is the Coriolis parameter of 


order 1071 s! (Tennekes, 1973). The largest eddies in the layer will be of height 


148 


equal to hp; and rotate at xhpj/U, where U is the PBL mean wind speed, 


assuming the eddy is circular. Then 


_ "hei 
5 Ü (155) 
Or 
oi 
L S f VCD (134) 


which is about 5min for Cp equal to 103. Therefore, to measure these larger 
scale fluctuations requires time records of at least 10min to sample at the 


Nyquist frequency. 
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Figure 42. Normalized Stress Cospectrum for Varying Stability after Kaimal 
et al., (1972) (from Panofsky and Dutton, 1984) 


The above considerations are summarized as follovvs: the turbulent 
covariance is adequately measured with a sampling interval of about 20min, 
which is twice the Nyquist frequency of the largest eddies in the boundary 
layer, and also still within the spectral gap between the synoptic and turbulent 


peaks in the wind speed spectrum. The fact that 20 minutes is at the upper 
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limit of the spectral gap means the stress values will be slightly 
underestimated. Nevertheless, in this study a 20min averaging period was 
selected for the computation of the wind stress. 

The wind speed averaging interval (microscale) was chosen at first 
for convenience to be 10s to correspond to the Ku-band radar azimuthal 
movements. However, this short time interval introduced significant 
variance in the computed short term stresses so that one minute running 
mean wind speeds, centered on the radar measurements, were subsequently 
used. The selection of a one-minute average was based on a calculation that 
showed the one-minute means to reduce the variance in the wind speed by 
50%. 

The choice of averaging interval for the wind stress determines the 
accuracy of the computed stress. Wyngaard (1973) recast Lumley and 
Panofsky's (1964) results to give an expression for the averaging time, Tu'w' 
needed for a specified accuracy in the covariance. For neutral conditions at a 
measurement height, z, wind speed, u, and specified accuracy, €, the 


averaging time is given by 





T ass = 7 (135) 


For lower wind speeds, the averaging time is longer unless 
compensated for by reducing the measurement height. Conversely, for a 
given wind speed, the averaging interval is longer at higher measurement 
heights. With Tum = 20min, 6 U=7.5mps, the average accuracy of 


any particular covariance computed in the present study is +16%. Under 


very unstable conditions, the numerator in (135) increases to 100z, therefore, 
the accuracy is reduced to +35% of the true value. 

With the acceptance of these accuracies comes the question of how 
exactly to calculate the averages. The averaging scheme must produce 
continuous averages which are differentiable to any order required by the 
equations of motion and they must satisfy Reynold's rules of averaging. Here 
it was decided to use a moving 20min average evaluated every 10s in order to 
maintain continuity of the drag coefficient between independent 20min 
records and to allow correlation of the computed drag coefficient with a 
particular 1min (microscale) wind speed. This averaging method allows each 
Ku-band azimuth bin to be assigned a unique stress value. 

As Panofsky and Dutton (1984) point out, running averages are 
continuous but they do not satisfy Reynold's postulates exactly. They 
maintain that the defect is not serious as long as the averaging period is 
chosen in correspondence with the spectral gap. 

3. Results of the Method 

Using the 20min us and wind speed values, a 20min drag coefficient 

is computed and corresponding one-minute stresses (Figure 43a-c) , as well as 


a 20 min running mean bulk Richardson number (Figure 43e). 
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Figure 43. Wind and Stability Parameters as a Function of Time: (a) Wind 
5tress, (b) Drag Coefficient, (c) Wind Speed, (d) Wind and Stress Direction, (e) 
Bulk Richardson Number and (f) Relative Angle between Wind and Stress 
Directions 
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Variations in the azimuthal dependence of the NRCS may be due to 
relative differences between the wind direction and wind stress vector. 
Therefore, the wind vector and wind stress vector directions were computed 
and plotted (Figure 43d,f). In the example shown, the drag coefficient has a 
mean value of about .0018, no trend, but significant temporal variations, 
which span a range of values from .0012 to .0026. The twenty minute drag 
coefficient follows the variation of the stress as normalized by the mean wind 
speed. Similarly, the one-minute stress has the overall character of the 
longer-term 20 minute drag coefficient, but it also incorporates the effect of 
the short time scale variations in the wind speed. 

The atmospheric stratification was stable for this data run as 
indicated by the positive value of the Richardson number. Under stable 
conditions, Geernaert (1989) found that the stress vector is often directed to 
the left of the wind vector as is shown in Figure 43f. 

4. Conversion to a Neutral Wind Speed 

The wind speed predictions of the existing scatterometer model 
functions are in terms of a neutral wind speed at 19.5m. Therefore, 
comparison of the results of this experiment to model predictions requires 
that the measured wind speeds be converted to their corresponding neutral 
values. 

The procedure is to first convert the measured 20 minute drag 
coefficient to its neutral value using the integrated stability profile, equation 
(87). Then, the one-minute neutral wind speed at the measurement height is 


given by 


1 


EH (136) 


vvith the one-minute friction velocity, u,. At any other level, z, 


15 


| de wd 3 
u,(z)= 2. ( (197) 


The stability parameter, z/L, is required for the integrated stability 
function (83) and in this study it vvas obtained from equation (80) using the 


empirical relationship between z/L and Rig established by Donelan (1974): 


2 7 ine O 
L 


= 16.0Ri, Ri, > 0 (138) 
for -0.05SRips0.02 and wind speeds 4«U(10)«17m/s. Here, the twenty minute 


running mean value of the bulk Richardson number was used in (138). 


D. COMPUTATION OF THE WAVE FIELD PARAMETERS 
1. One-Dimensional Wave Spectra 

One-dimensional wave spectra were computed using measurements 
from the wave staff at the center of the pentagonal array. The laboratory 
calibration equations converted the wave staff output voltage to surface 
elevation and fast Fourier transforms were used to evaluate the spectra 
(Bendat and Piersol, (1986)). Since the data runs varied in length from 
fifteen to ninety-five minutes, the 20Hz data were blocked in 51.2s records 
(1024 points) to ensure sufficient degrees of freedom for the computation of a 
reasonable spectrum for even the shortest runs. 

First order descriptive quantities obtained from the one-dimensional 


spectra were the frequency (period) of the spectral peak, fp, and peak 


wavelength, Lp. The wavelength was computed from the linear wave 
dispersion relation, using the peak frequency and the numerical method of 


Wu and Thornton (1986); 1.e., for a given wave period, T, the wave length is 


= E T^tanh = (139) 
7 


jocum 
л 


Р 


where d is the water depth. The significant wave height was computed 
from the standard deviation of the surface elevation, o 


5 (140) 


2 


Second order quantities computed vvere the significant slope 


o 
§= T (141) 
t 
peak phase speed 
2: fe, (142) 
and wave age 
C, 
W = Tm (143) 


for wind speed, U. These quantities are descriptive of the degree of wave 
development: longer waves have smaller slopes and faster propagation 
speeds. 
2. Two-dimensional (Directional) Wave Spectra X 
A limited number of directional spectra were computed using the 
full array of wave staffs and the maximum likelihood processing method 


outlined in Tsanis and Donelan (1988). The one-dimensional spectra vvere 


frequently noted to be bimodal at low wind speeds (containing low and high 
frequency peaks), and under these circumstances the radar cross-section 
exhibited anomalous azimuthal behavior. Therefore, the purpose of these 
directional spectra was to confirm that the long wave (low frequency peak) 
was traveling down the long axis of the lake and the short wave (high 
frequency peak) was more or less in the wind direction. The hypothesis was 
that if the swell was propagating at a moderate angle to the wind it may 
modulate the short wind waves and Bragg waves, hence the azimuthal 


dependence of the radar backscatter. 


E. SUMMARY 

This section has described the data processing methods used to convert 
the WAVES87 raw data into useable geophysical and radiometric quantities. 
The final data set consists of the Ku-band (10s) and C-band (1min) NRCS 
temporally collocated with wind speed, air and water temperatures, 
significant wave height, bulk Richardson number, stability parameter, 
neutral wind at several heights, wind stress, and measured and neutral drag 
coefficients. Each data run also has associated with it a one-dimensional 
wave spectrum from which peak frequency, wavelength, significant slope, 
phase speed and wave age are computed. A set of twelve directional spectra, 
representative of a range of wave conditions, was obtained for the initial 


examination of long wave effects on the radar cross-section. 


V. DISCUSSION AND RESULTS 


A. INTRODUCTION 

Methods by which the raw data were processed to obtain the radar cross- 
sections and collocated geophysical parameters of interest were described in 
the previous chapter. Random and systematic errors of each variable must be 
kept in mind as the following analyses and interpretations are performed. 
An additional source of variability is that the correlations between the radar 
measurements and the geophysical variables involve different temporal and 
spatial averaging schemes. In spite of these caveats, clear dependencies of the 
radar cross-section on environmental parameters will be demonstrated. 

The organization of this chapter follows the schematic diagram of 
Figure l. The behavior of the NRCS with wind speed will be discussed first, 
followed in turn by discussions of the observed effects of wind direction, long 


waves, atmospheric stratification and water temperature. 


B. NRCS DEPENDENCE ON WIND SPEED PARAMETERS 

As noted in the discussions of the scatterometer model functions, there is 
some controversy regarding the appropriate wind speed parameter to relate 
to the NRCS (i.e; traditionally the wind speed, neutral wind speed or wind 
stress). Ideally, the most appropriate wind parameter should reduce the 
number of dependent variables; in particular, eliminate the stratification 
dependence and reduce the wave dependence. The possibilities are then to 


use what is known as the neutral wind at some level (equation 137) or to use 


directly the wind stress at the surface. As a measurement convenience, 
existing model functions predict neutral wind speed at 19.5m although it 1s 
quite certain that this rather artificial quantity is not responsible for the 
generation and growth of the capillary waves to which the NRCS is sensitive. 
On the other hand, the (neutral) wind speed very near the sea surface is 
included dynamically in some expressions for the short wave growth rate, 
which hold that it is the ratio between the wind speed and phase speed of the 
waves which amplifies the waves. 

The difficulty with the choice of the neutral wind at any level is that its 
determination is based on a particular surface or boundary layer model: i.e., zo 
or drag coefficient. Air-sea interaction specialists have yet to identify an open 
ocean standard which normalizes the drag coefficient for wave state in 
addition to measurement height, wind speed, and stability. Even with a 
satisfactory standard, it is not necessarily possible to extend the assumed 
.logarithmuic profile directly to the surface since the waves influence the wind 
field above them up to about one wave height (Snyder, et al., 1981). 
Therefore, correlation of the wind just above the surface with radar 
backscatter is hindered by the difficulty of measuring the near surface wind. 

The alternate parameter of choice is then the friction velocity or wind 
stress, which by virtue of its definition, is stability independent. However, 
verification of a wind stress algorithm for the scatterometer has the same 
verification problem of the drag coefficient/neutral wind verification 
problem mentioned above. While the use of wind stress should 
hypothetically eliminate a second order dependence on stability, the influence 


of wave state on the resulting model function based on friction velocity may 
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increase. Note that the stress represents the slope spectrum of all surface 
waves while Bragg scattering is due to waves within a small range of 
wavenumbers; i.e., wave modulation of the stress spectrum will differ from 
the degree of Bragg wave modulation by longer waves. At some point a 
major effort to measure the wind stress at sea will have to be conducted in 
order to establish the open ocean drag coefficient to which all other 
measurements can be scaled. 

In the present study, the wind stress was measured using a bivane 
anemometer and processed as described in the previous chapter. The friction 
velocity values were used to determine the drag coefficient and then the 
neutral wind. All of the neutral winds were based on stability corrections 
which used the bulk Richardson number to estimate z/L as in (80). Friction 
velocity values less then 0.12m/s were omitted from the plots and data 
interpretations (wind speed too low for bivane), as were points for which the 
relative humidity was above 90% (probable rain) and two runs which were of 
questionable value (one of short duration and one contaminated by specular 
returns from what appears to be a solid object). The NRCS values were all 
reduced by 1.0dB in the following plots to account for the wider beamwidth of 
the tower radar to allow comparison with the satellite based model functions. 

The data were stratified by upwind, downwind, crosswind near 90° and 
crosswind near 270° relative wind angles with an allowance of 5° about each 
prescribed direction. The vertically and horizontally polarized "smoothed" 
radar cross-sections were then plotted as a function of the measured friction 
velocity (Figures 44, 45) and the neutral wind at 19.5m (Figures 46 and 47 for 


comparison with the model functions) and the neutral wind at 1m (Figures 


48 and 49). The 1m neutral vvind level vvas employed instead of one-half 
Bragg vvavelength in order to be more realistic about the application of the 
logarithmic vvind profile in the presence of the vvaves. 

The figures shovv the dependence of the cross-section on the vvind 
parameters for each of the incidence angle settings of the Ku-band radar. 
There is still considerable scatter in the data, hovvever, the scatter of the data 
points and the slopes of the lines at each angle are consistent vvith previous 
results. At nadir there is a small negative dependence on vvind speed or 
stress, reflective of the increased scatter of the incident radiation as the surface 
becomes roughened by the vvind vvavesl. The ten degree incidence angle is 
nearly wind speed independent although the best fit least squares line does 
generally assign either a small positive or negative slope value. The slopes of 
the lines increase with incidence angles up to 60° (the downwind slopes being 
the exception). At 60° the h-pol data begin to show sensitivity to incidence 
angle by becoming more scattered. The 80° data still display some dependence 
with wind speed in the upwind direction for the vertical polarized data; the 
horizontal polarized cross-sections at 80° generally fall below the noise level 
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1 This data set could be used to examine the wind speed algorithm used 
in altimetry, especially as regards the significant wave height contribution to 
the measured cross-section. It should be noted here, however, that the data 
editing criteria used in the present analysis no doubt "overedit" at nadir 
where the cross-sections are more a function of the long waves than the 
vagaries of the wind stress. In addition, since the backscatter is isotropic in 
direction, all of the data points should be used irrespective of wind direction 
instead of being stratified by relative wind angle. 
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Comparison of the slopes of the neutral vvind speed lines vvith previous 
results shovv that the radar cross-sections from Lake Ontario display stronger 
dependence on the vvind speed. A compilation of previous results for 
upwind v-pol data is reproduced here from Phillips (1988) in Table VI, where 
the Ku-band data are those listed under 13.9 and 14.6GHz. The slopes at 20° 
and 40° differ by about 0.3 to 0.7 and by about 1.0 to 2.0 at 60°. The slope 
coefficients (Figures 44 and 45) on the ustar curves at the lower three 
incidence angles (0°, 10° and 20°) are comparable to the neutral wind speed 
slopes, but are significantly lower at 40° and 60°, e.g., 1.6 v. 2.4 upwind, v-pol. 

The correlations on friction velocity are poorer than the corresponding 
neutral wind speed values. In this analysis the neutral wind speeds are 
derived from the measured wind stress, but they have been smoothed by the 
ratio of the measured and neutral drag coefficient in the conversion process. 


That is, the neutral wind speed is defined by 
(1, (11.5) = Coin, (144) 
Here, 


u. = C;U(11.5) (145) 


(62 


TABLE VI. COMPARISON OF NRCS VERSUS WIND SPEED SLOPE 
COEFFICIENTS WITH PREVIOUS STUDIES 
Partial table reproduced from Phillips (1988). 


Radar frequency 


(GHZ) 139 139 13.9 
Rada: wavelength 
icm} 2 16 2 16 2 16 
Incidence angie 
Tw - 047 — 
30° 139 140 1 
161 (32°) 
40° 149 204 218 
1 3 (43°) 
50" 169 217 23 
132 
60° — 072 — 


Summary of Upwind Slope Coefficients from Present Study 


Incidence V-POL H-POL 
Angle Scaled Scaled 
u, Ux(19.5) un(1) Ux(19.5) uis un(19.5) UNÜ) un(19.5) 
20: (0 LI 1 .98 Ri 1.1 1.1 .97 
40 1.7 2.4 7 1.6 1.7 2.4 22 1 
60 1.6 27 27 13 1.3 20 2.4 bi 
& 50 1.4 1.6 11 -- — -- -- 
Summary of Upwind Correlation Coefficients from Present Study 
Incidence V-POL H-POL 
Ane Scaled Scaled 
u, u.(19.5) ux(1) Ux,(19:5) u, ux(19.5) un(1) UN(19.5) 
° .76 0.85 .80 7 74 .82 .78 74 
4 ./9 93 SCH .70 .78 .92 .90 .65 
6C 76 .85 .80 55 7 .82 57 .42 
& 70 73 .66 ور‎ -- -- — -- 
Summary of Upvvind Exponent Standard Errors from Present Study 
Incidence V-POL H-POL 
Anrle Scaled Scaled 
u. uN(19.5) ux (1) un(19.5) u, ux(19.5) uN(1) ux(19.5) 
2 .08 0.08 109 09 08 08 10 07 
1 16 dü E .19 .16 12 13 10 
60 217 .20 .25 .23 .16 .20 .24 21 
5 2 27 .26 T3 -- -- - -- 
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Therefore, the neutral vvind is obtained from 


U, (11.5) = E U(11.5) (146) 


DN 


This expression shows that the variability of the neutral wind is reduced, 
since the ratio of the measured and neutral drag coefficients is nearly equal to 
one except for very stable conditicns for which the ratio decreases to about 0.5. 
Therefore, instead of being susceptible to variability in both the drag 
coefficient and the wind speed, the neutral wind variability is due mostly to 
changes in the wind speed, since stability is accounted for in the neutral drag 
coefficient. 

If the friction velocity slope coefficients are accepted as valid, then the 
differences in slopes between the upper and lower incidence angles and the 
differences from the neutral wind values can be interpreted in terms of the 
sea state. At nadir and 10°, the cross-section values are due to scattering from 
a rough surface; 1.e., basically the wind waves, not the Bragg-length capillary 
waves. However, as the incidence angle steps into the Bragg domain (say 20° 
to 60°), the capillary waves become the critical feature of the sea surface 
responsible for the observed returns. These waves of slow phase speed and 
large steepness are extremely responsive to changes in wind speed or stress 
and to the background state of the sea. One quantity which links all of the sea 
States is the drag coefficient, and the observed variations in the slopes of the 
wind stress/wind speed lines imply that this coefficient has a sea state 


dependence. 
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Recall the relationship betvveen the vvind stress, friction velocity and 


neutral vvind speed, 


әти. (147) 


о | 


If the neutral drag coefficient vvere only a function of vvind speed, then 
according to this expression, the stress vvould be proportional to a povver lavv 
dependence on Un with an exponent greater than 2 and approaching 3. 
Similarly, if the NRCS is proportional to the wind stress then it should also 
have the same power law dependence on the neutral wind speed. The NRCS 
dependence on neutral wind is observed to be generally greater than two and 
approaching three at intermediate incidence angles which strongly suggests a 
stress dependence. However, the fact that both the friction. velocity canela E 
neutral wind speed slope coefficients are less than what 1s-expecteslgsibo 
suggests a different functional dependence for the drag coefficient; i.e., a 
dependence on sea state. 

The possibility of a sea state dependent drag coefficient can be examined 
by comparing the lake NRCS versus wind speed slopes to those predicted by 
the model functions discussed earlier. The models were run at 40°, v-pol for 
windspeeds ranging from 1 to 15m/s. The cross-sections from each of the 
models and the lake observations are shown in Figure 50, where the line 
labelled "Lake normal" is the best fit line shown in Figure 46. The vertical 
error bar at UN(19.5) = 5.5m/s is the +1 standard deviation of the absolute 
calibration error of the radar as determined from the ten degree data; i.e., the 


level of either of the lines could go up or down by about 1.9dB. 
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The best fit line of the lake observations (Lake normal) is considerably steeper 
than the corresponding model function predictions. A possible explanation 
for this result is found by noting that all of the model functions displayed use 
the Large and Pond open ocean drag coefficient to convert measured vvind 
speed to stress and/or neutral wind speed. However, in the lake the neutral 
drag coefficient has a greater dependence on wind speed than the Large and 
Pond results for the open ocean. This is due to the fact that the lake waves are 
often in an active growth stage and even the longest waves are quite steep by 
comparison to open ocean swell (recall Figure 12). Therefore, for a given 
wind stress (and hypothetically, also for a given NRCS) the neutral wind for 
open ocean conditions is higher than the corresponding lake value through 


the scaling 


i 
C ۱ “ 
CON -| -— | Oy Late (148) 


Con Ocean 


where the ratio of the Large and Pond (88) and Donelan (89) drag coefficients 
has been used as a first guess scaling factor. The line labelled "Lake scaled" 
illustrates the result of this scaling at 40° incidence and Figures 51 and 52 
show the results for all of the incidence angles. The slope of the scaled line is 
now parallel to the open ocean results implying that the sea state difference 
between the lake and ocean has an effect and should be incorporated in the 
drag coefficient. This also implies that the wave influence on both the stress 
and NRCS is similar but not the same; i.e., for a given Uy, the NRCS and u. 
both increase with increasing long wave slope. The scatter in the converted 


points in Figures 51 and 52 is a ramification of the fact that the Large and 
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Pond and Donelan vvind speed dependent drag coefficients explain only a part 
of the variance in the vvind stress. 

In summary, the radar cross-sections have been analyzed as a function of 
neutral wind speed and wind stress. The overall behavior of the cross- 
section vs. wind parameter curves are consistent with previous results which 
show the dependence to be negative at nadir, transitioning through zero 
slope at 10°, and becoming positively correlated above 10° incidence. The 
wind stress and neutral wind speed exponents are similar at lower incidence 
angles but the wind stress exponent is smaller than the wind speed exponent 
at the intermediate, Bragg-scatter incidence angles. This observation, and the 
steepness of the lake wind speed curve relative to the open ocean model 
function curves, suggests that the cross-section in the Bragg domain is 


affected by the spectral characteristics of the non-Bragg waves. 


C NRCS DEPENDENCE ON WIND DIRECTION 
1. General Observations 

The second most characteristic feature of the scatterometer signal is 
its cosine-like dependence on the relative angle between the radar antenna 
and the wind direction (recall Figure 2b). This behavior was first examined 
empirically with a series of aircraft flights in which circular flight tracks were 
used to generate the relative angles between the wind and radar (Jones, et al., 
1977). Since that time, aircraft scatterometers with rotating or phased array 
antennas have been developed so that the airplane can fly a standard straight 
line "flight plan" thereby allowing more control over the GER of ocean 


features to be examined (e.g., Fronts and Air-Sea Interaction Experiment, Li, 


et al., 1989). 


In the WAVES87 experiment, the Ku-band scatterometer was 
mounted on an antenna rotator and extended away from the tower about 8m 
on a boom. The antenna rotator scanned about 300° in azimuth every ten 
minutes at set incidence angle settings as summarized in Table V (Chapter 
Ill). A full 360° rotation was desirable but impractical since the angles behind 
the radar would have intersected the tower and automatic operation of the 
radar would have required an additional level of sophistication to allow for a 
certain number electrical cable wraps before the antenna was reversed and the 
cabling unwound. Even with this angular gap of about 60°, the sinusoidal 
behavior of the cross-section is evident. 

The observed "smoothed" cross-sections at each incidence angle and 
polarization were used to create three-dimensional surfaces of the cross- 
section as a function of wind stress and relative angle (Figures 53 and 54). The 
nadir surface is characteristically flat with azimuth; the cross-section scale 
used here is too large to reveal the small negative dependence on wind stress. 
At 10° incidence angle there are frequently large azimuthal modulations of 5 
dB or more, which are not predicted by existing model functions. The cosine- 
like dependence becomes apparent for both polarizations at 20° and 
increasingly obvious at the higher incidence angles, especially for v-pol. The 
h-pol surfaces are rougher and the 60° and 80° data for this polarization are 
basically noise. The surfaces at 20°, 40° and 60° show the surfaces to have a 


positive slope with respect to the friction velocity. 
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SIGMA-O V. USTAR AND CHI FOR VARYING ELEVATION 
Ku-band Vertical Polarization 
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Figure 53. Three-dimensional Surfaces of Vertical-Polarized NRCS asa 
Function of Incidence Angles, Friction Velocity and Relative Angle, Chi, 
(Wind Minus Antenna Angle). Incidence Angle Increases Clockwise from 
Nadir at Top Left, the Relative Angle is Given in Radians with r 
Corresponding to the Downwind Look Direction 


177 


WAVES 87 
SIGMA-0 V. USTAR AND CHI FOR VARYING ELEVATION 


Ku-band Horizontal Polarization 





Figure 54. Three-dimensional Surfaces of Horizontal-Polarized NRCS asa 
Function of Incidence Angle, Friction Velocity and Relative Angle, Chi 
(Wind Minus Antenna Angle). Incidence Angle Increases Clockwise from 
Nadir at Top Left; the Relative Angle, Chi, is in Radians with x 
Corresponding to the Downwind Look Direction 
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After this holistic vievv of the vvind and azimuth dependence, the 
smoothed cross-sections for each individual data run, incidence angle other 
than nadir and polarization were fit with a three term truncated Fourier 
series as described by equations (91). These observations of the azimuthal 
dependence of the scatterometer are unique as tower-based radar 
measurements; therefore, all of the available azimuthal plots at the various 
incidence angles for both polarizations are provided in Appendix B. Some 
examples of these fits are shown in Figure 55 (v-pol, 20°) in which it is seen 
that azimuthal variations at moderate wind speeds are apparent, but the 
variation is frequently asymmetric and the crosswind minima are not 
necessarily at 90° and 270°. The upwind/downwind asymmetry is more 
obvious for h-pol and increases with incidence angle (see Appendix B). Since 
the pure cosine fit cannot simultaneously fit both the upwind/downwind 
asymmetry and the position of the minima, the largest errors in minima 
placement occur for greater differences in the upwind/downwind cross- 
section values. Some of this asymmetry was hypothesized to be induced by 
the directionality of the wave field. As an attempt to quantify these phase 
shifts for comparison to wave spectrum quantities, a second order sine 
component was added to the three-term cosine series. 

f oil CUrvesmincluding themsine component improve the 
placement of the minima (cf., runs 80 to 92, Figures 55 and 56), but there are 
still problems with the overall fits due to data gaps and spikes and to 
relatively long sweep rate of the antenna; i.e., the field evolves in the twenty 


minute sampling period. The SASS1 model function team tested fits to 
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similar data by computing 36 terms in the series only to find that the first 
three explained most of the variance in their data (Schroeder et al., 1982). 
Recently in an examination of the azimuthal dependence observed with an 
aircraft scatterometer, Li et al. (1989) used an iterative method to locate the 
minima in which a truncated Fourier series of three sines and a cosine was 
used to estimate the position of the minima, followed by a polynomial fit to 
all points 45° on either side of the estimated minimum. Donelan and Pierson 
(1987) point out the fallacy of fitting the curves with a truncated series of any 
degree: inasmuch as the data points are not distributed evenly over 360° 
adding additional terms will not improve the fits in a general sense. 

The azimuthal behavior shown in the WAVES87 data indicates the 
model functions which use the cosine series approach to predict the wind 
direction may be in error by ten degrees or more. The baseline configuration 
of all the models discussed previously results in predictions of symmetric 
dependence of the cross-section with relative wind angle; i.e., the minima are 
near 90° and 270° and separated by 180°. The physically based model functions 
contain coefficients which can be adjusted to examine wave effects such as 
wave slope and hydrodynamic modulation of the placement of the minima. 
Rather than perform this type of sensitivity test, it was decided to examine the 
behavior of the radar under different wave conditions as indicated by the 
directional wave spectra. These analyses will be discussed in the Section D. 

2. Rain Effects on the NRCS 

Entries in the experimental log indicate high rain rates during data 

runs 143 to 146. Prior to these runs, the wind had been blowing at moderate 


to high speeds from the long fetch direction for about two days. This relative 


constancy of speed and direction allovved the vvave spectrum to become fuliy 
developed with a peak frequency of 0.16Hz and peak wavelength of about 
60m. The direction of the primary wave was nearly in alignment with the 
wind, down the long axis of the lake. 

Under these circumstances the 10? and 20? cross-sections (Figure 57) 
show a broad peak centered on the downwind direction (180?), indicative of a 
low amplitude swell propagating in the wind direction. The scenario 
suggested by the 40° data is that between runs 140 (wind speed of only 1.3m/s) 
and 143 (7.8m/s), a wind sea had begun to develop as evidenced by the near 
perfect sinusoidal shape of 143. The rain then started toward the end of 143 
and continued through run 146 (12 hours). It appears that the rain drops 
augment the small, wind-generated capillary waves (of Bragg wavelength 
approximately 1cm) by creating new rain-generated capillary waves? The 
resulting surface is then isotropic in wind direction (Stoker, 1977; LeMehaute, 
1988) and the distribution and number of capillary waves is radically altered. 
As a result, the azimuthal dependence of the cross-section is obliterated 
(Bliven, et al., 1989) and the level of the cross-section 1s increased above that 
which is normal for the particular wind speed (cf., runs 134 and 144) (Moore 
et al., 1979) suggesting that there were additional Bragg-length scattering 
waves on the surface.. 

As the rain continues, it is possible for a thin rain induced mixed 


layer of relatively high eddy viscosity to form which damps the longer waves 


2 Path-loss of the Ku- and C-band signals through rain is negligible at the 
radar ranges used in this experiment (Battan, 1973); therefore, this 
interpretation of the rain effects in terms of the scattering surface is valid. 
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Figure 57. Azimuthal Dependence of the NRCS in the Presence of Rain: a) 
10”, b) 20”, (Runs 143 to 145) and c) 40? v-pol for Runs 143 to 146. Note 
Attenuation of Azimuthal Modulation at 40° between Runs 143 and 144. 
Rain Started near the End of Run 143. 
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Figure 58. VVave Amplitude Spectra Shovving Decrease of Spectral Density of 
VVavenumbers above the Spectral Peak from Runs 143 to 146. 
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present (Nystuen, 1989). The one-dimensional frequeney spectra for runs 143 
to 146 are overlaid in Figure 58, from vvhich it is seen that the frequencies 
above the spectral peak decay as the rain continues. Therefore, runs 145 and 
146 shovv no azimuthal dependence and the level of the curve falls as the 
longer vvaves are damped and the vvind decreases. This scenario is 
appropriate at the C-band vvavelength. İn this case the rain damps the longer 
Bragg vvaves (3cm) but the rain drops are not large enough or distributed 
appropriately to create 3cm vvaves in their place. The corresponding C-band 
data for runs 144 to 146 shovv the cross-section to decrease steadily vvith the 
continuing rain and decreasing winds, rather than showing an initial increase 
as with the Ku-band NRCS. 

This section has demonstrated that although the overall behavior of 
the cross-section as a function of relative wind angle is a cosine-like 
dependence, there appear to be circumstances under which the sinusoid 
becomes asymmetric or even isotropic. The observations do not support the 
| predictions of existing model functions which show cross-section minima to 
be exactly at crosswind under all conditions. The polarization of the incident 
radiation is also a critical parameter: the v-pol data often showed an 
azimuthal dependence at 80* whereas the h-pol data were more scattered at 


60° and mostly noise at 80°. 


D. NRCS DEPENDENCE ON SPECTRAL WAVE PARAMETERS 

The azimuthal behavior of the cross-sections shown in the previous 
section displayed significant asymmetry. In this data set, several related 
physical mechanisms could be working together to produce this asymmetry. 


In some instances, the wind speed changed in the twenty minute period over 
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vvhich the measurements vvere made so that the level of the NRCS vvas 
altered. Since the radar scanned back and forth instead of always starting 
from one side, the averaged cross-section at the start/finish azimuth bin 
consists of two NRCS measurements twenty minutes apart, whereas the 
averaged cross-section measurements near the center of each run consists of 
NRCS measurements taken within minutes of each other. This asynoptic 
measuring scheme usually results in an "odd" function of the NRCS with 
respect to the upwind or downwind maximum (see annotated Run 80, 40°, v- 
pol, Appendix B). 

Alternatively, the azimuthal asymmetry may be due to different wave 
conditions; 1.e., varying spectral forms from one run to another. To express 
this in terms of the two-scale model, the anisotropy of the cross-section with 
direction is due to changes in the large scale slope of the tilting surface and to 
the directional spread of the small scale waves about the wind direction. In 
this section, the azimuthal behavior of several data runs will be discussed 
within the framework of the spectral composition of the waves present as 
determined from frequency and directional spectra and descriptive 
parameters of the sea state. The observed effects are those due to long waves 
relative to the Bragg waves. 

Three classes of wave spectra, having progressively more complex wave 
fields, are chosen for comparison to the azimuthal variations of the radar 
cross-sections: 1) a moderate wind speed (7.5m/s) wind wave sea (4m waves 
aligned with the wind, Runs 49 and 89); 2) a combined sea /swell in which the 


swell component opposes the sea (Runs 79 and 80); and 3) a moderate wind 
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speed, "confused" sea of broad frequency and directional bandwidths (Run 
— 
1. Class 1: Moderate Wind Speed and a Wind-driven Sea 

This test case is presented as the "norm" which satisfies the basic 
assumptions of scatterometry upon which the model functions are based. 
That is, the ind and waves are aligned and the waves are locally generated 
wind waves. Two runs of similar wave conditions will be examined as a 
reliability measure. 

The averaged environmental conditions during the 70 minutes of 
Run 49 (Figure 59) were 8.8m/s wind speed from 231°T. The significant wave 
height was about 25cm and the atmospheric stratification as inferred from 
the air-water temperature difference (+5.5°C) was stable “Under cce 
conditions, the NRCS exhibited the azimuthal behavior relative to the wind 
direction shown in Figure 60a. The curve is flat at nadir, the 10° curve is 
different on either side of 180° due to a decrease in wind speed in the first five 
minutes of the 10° scan and sinusoidal curves are obvious at all other 
incidence angles. Even though the downwind maximum is in the tower gap 
and therefore not measured, it appears from the symmetry of the curves that 
the upwind/downwind difference is greater for h-pol than v-pol. This 
difference as well as the basic difference in NRCS magnitudes between h-pol 
and v-pol increases with increasing incidence angle (recall the reflection 
coefficient of horizontal polarized radiation has a stronger incidence angle 
dependence, (11a) and Figure 5). 

A perspective diagram of the directional wave spectrum associated 


with this run (Figure 60b) shows a4m wind wave sea propagating with the 
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Figure 59. One-Minute Block-Averaged C-band NRCS and Environmental 
Data from Run 87049. Time Associated with Particular Incidence Angles o? 
Ku-band Radar are Annotated on Wind Speed Plot 
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wind at about 225°T. The one-dimensional frequency spectrum for Run 49 
shows a bifurcated peak with two nearly equal frequency (.69Hz, .0.62Hz; 4m 
peak wavelength), equal energy wave components. The directional spectrum 
assigns directions of 220°T and 230°T to these waves; therefore, on the polar 
plot of the wind, wave and radar angles (Figure 60c), the direction of the peak 
wave is shown as a 10° arc rather than a single angle. The upwind maximum 
of the NRCS is within 5° of the wind/wave direction which is the best 
alignment one can expect in this experiment, since 5° is the minimum 
resolution of both the radar and directional spectra. The NRCS minima are 
separated by 160° and are 75° and 85° on either side of the maximum. The 
positions of the minima approximate the directional spread of the wind 
waves but the correlation cannot be made exact since 1) the directional 
spectrum can be altered by a change in the wave staff data processing 
parameters, and 2) the spectrum has no absolute minima. 

The pure wind sea case allows the most legitimate comparison of the 
measurements with the model function predictions. All of the model 
functions were run for an average neutral wind speed of 8.8m/s and wind 
direction of 230°, corresponding to the conditions for Run 49, 40? v-pol. The 
NRCS values computed from the model functions are illustrated in 
Figure 61 from which it is seen that there is significant variation among the 
model functions with the SASS1, D&P and D&V models forming a set, 
separate from the SASS2 and Plant models. The D&P and D&V models were 
tuned to the same in-situ data set so it is not surprising that they produce 
similar values; the SASS2 was constrained by SEASAT data, SASS1 by other 


in-situ. data and the Plant model was not tuned to any particular data set. The 
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Figure 61. Azimuthal Modulations of Model Functions (o ~ Plant, x — 
SASS2, + - D&P, -.- line SASS1, solid line D&V) and Measurements from 
Run 49 ( * ). 
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amount of upwind/crosswind modulation is least for SASS2 (5dB) and about 
7dB for the other models; the upwind/downwind modulation varies from 
about 1 to 2dB (Plant, SASS2; D&P, D&V, SASS1). The overall shapes of the 
azimuthal modulation are all similar except for the Plant model which has 
broad, flat-topped maxima and sharper nulls at the crosswind minima. 

The NRCS data points from Run 49 are superimposed on the model 
function predictions after subtracting a constant 5dB bias from the 
measurements to account for the lake/ocean wind speed scaling as discussed 
in the previous section. The sinusoidal pattern of the data points follows 
that of the model functions except the minima are located at relative wind 
angles of about 90° and 285°, instead of being located at 90°, 270° (SASSI, 
SASS2, Plant) or 95°, 265° (D&P, D&V). The cluster of three points of 
unusually high magnitudes near 110° are due to a data spike of unknown 
Origin and do not change the location of the minimum at 90°. The distance 
between the relative wind angle minima in the measurements is 195° which 
corresponds to the angular separation of the minima in the wind wave 
directional spectrum. Therefore, there is a difference of 15° to 20° between 
the directional dependence of the model functions and the radar field 
measurements. 

Run 89, which has a similar wind wave spectrum as Run 49, had an 
averaged wind speed of 7.5m/s from the short fetch direction, 230°T (Figure 
62); however, the atmospheric stratification was unstable. The smoothed 
NRCS for all incidence angles and both polarizations are plotted‘ against radar 


azimuth in degrees True North (Figures 63 and 64) from which it is seen that 


1. 


all of the incidence angles above nadir are azimuthally modulated in 
increasing amounts. 

The maxima of the modulated curves fall in the vicinity of the 
upvvind direction, vvith the 60? v-pol curve having the best conformance in 
shape vvith respect to the vvind direction. The 40? v-pol curve is somevvhat 
skewed and the 80° v-pol curve has a broad peak of relatively constant NRCS. 
As before, the difference between the h-pol and v-pol NRCS increases with 
incidence angle. The decrease in h-pol is dramatic enough that even at 
7.5m/s the 60° h-pol values are frequently below the noise level of the radar. 
The polar plot of the radar wind and wave angles, and the perspective plot of 
the directional spectrum (Figure 65) confirm that in the case of a pure wind 
sea the NRCS is maximum in the upwind direction, perpendicular os 
wind wave crests which carry the Bragg scatterers. The minima are not 
necessarily 180° apart nor 90° from the wind direction; rather, the angular 
distance between the minima approxima’ the angular spread of the wind 


wave directional spectrum. 
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WAVES 87 --- November 19 Run 87089 Time 1301 
AvES : | inin blocks Overull : Moderate Wind Short Fetch Unstable 
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Figure 62. One-Minute Block Averaged C-band NRCS and Environmental 
Data From Run 89. Times Associated with Particular Incidence Angles of Ku- 
band Radar are Annotated on Wind Speed Plot 
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2. Class 2: Moderate VVind Speed, Combined Sea and Svvell 

A bimodal spectrum in which a nonlocally generated long wave 
propagates through a local wind sea is a frequent possibility in the open 
ocean. Due to the fetch geometry of Lake Ontario, a small change in wind 
direction produces a similar bimodal spectrum in which the peak frequencies 
and directions of the component waves are quite distinct. In this test case, 
two data runs for which the wave field consisted of two primary wave 
components are examined. 

The southern half of a semi-stationary, low pressure system was 
positioned over the lake for 19hrs on 17 and 18 November such that the wind 
velocity varied from about 7 to 10m/s and 190° to 170°T from the western to 
eastern end of the lake. The primary wave frequency and direction which 
developed were approximately 0.2s! and 60°T. The weather system 
progressed eastward on 18 November, and the wind direction and speed 
changed to about 230°T, 7.5m/s. 

Three hours after the shift in wind direction, radar and 
environmental measurements were acquired, at which time the 
environmental conditions were as shown in Figure 66. The wind speed was 
somewhat erratic, varying between 5 and 9m/s, but the wind direction was 
uniformly from 230° to 240°T. The significant wave height was about 30cm 
and the mean water level indicated that a seiche may have been induced by 


the movement of the pressure system across the lake. 
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Figure 66. 


WAVES 87 -—-- November 18 Run 87079 Time 304 
AvgS : 1 inin blocks Overall : Moderute Wind Short | Fetch SLable 
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One-Minute Block-Averaged C-Band NRCS and Environmental 
Data from Run 79. 


Analysis of the associated directional spectrum shovvs that a local 
vvind sea developed in addition to the long vvave in vvhich the total povver at 
the wind wave peak frequency is 76% of the power at the lower swell 
frequency. These two wave components propagated in nearly opposite 
directions: 66°T (44m long wave) and 221°T (3.2m short wave). The 
directional spread of the long wave is not well resolved by the maximum 
likelihood method for low spectral density (Donelan, pers. comm); therefore, 
only the peak direction will be used here. 

The azimuthal modulations of the NRCS for both polarizations and 
all incidence angles except nadir (Figures 67 and 68) still exhibit sinusoidal 
dependence with wind direction even in the presence of the long wave. The 
magnitudes of the NRCS values appear nearly equal to the wind sea case 
(Run 89), after allowing for wind speed changes at each incidence angle (the 
average speed is marked on each incidence angle curve). However, the 
position of the NRCS maximum, located at 245°T, could be considered to be 
aligned with the either the wind or long wave as illustrated in Figure 69. The 
minima are not symmetric with any one of the wind, short wave or long 
wave directions, instead they are located at the midpoints between the long 
and short wave directions and are separated by 180°%. This effect may be due 


to the increased spreading of the wind wave directional spectrum that occurs 


3 A fourth directional choice, not examined here, is the direction of the 
stress vector. In the presence of waves and for differing atmospheric 
stratification, the wind and stress vectors are not necessarily aligned 
(Geernaert, 1988b). 
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vvhen a vvind vvave sea runs against an opposing current, in this case the 
orbital velocity of the longer vvave component (Sugimori, et al., 1985). 

The vvind sea continued to grovv during Run 79 and by the time Run 80 
vvas completed, about three hours later, the total povver at the vvind vvave 
frequency was twice that of the long wave, Figure 70. With the wind wave 
dominating the spectrum, the maximum of the NRCS is in the wind 
direction and the minima (polar plot, Figure 70) are again only 165° apart. 

3. Class 3: Moderate Wind Speed, "Confused" Sea 

When the low pressure system described above moved across the 
lake, the wind changed from a light, northerly wind, that had been 
prevailing, to a southerly direction, then from the southwest. Runs 79 and 
80 described above were performed after the winds had settled into the final 
southwest direction. Between the initial northerly, and the final, 
southwesterly winds, the wave field was in a period of active readjustment to 
the changing winds. 

The environmental conditions of Run 75 were steady winds of about 
2— s from the south for the first 55 minutes of the data run; at 55 minutes, 
the wind speed started to increase while the direction remained southerly; the 
atmospheric stratification was stable (Figure 71). The azimuthal modulations 
of the NRCS with respect to True North (Figures 72 and 73) for the various 
incidence angles differ quite markedly from the previous cases. At 10* and 
20°, the modulation with wind direction is gone in both polarizations. The 
intermediate angles for the v-pol data (40° and 60°, Figure 72) show the peak- 


to-trough modulation to be smaller than the wind sea (only about 5dB) 


205 


compared to 10dB) and overall, the curves are more noisy.4 The maxima 
shift progressively toward the wind direction with increasing incidence angle 
while the cross-wind minimum stays near 270°. The h-pol curves at 40° and 
60° display even less modulation (Figure 73). 

The complexity of the wave/wind field is illustrated in Figure 74. 
The directional spectrum has three peaks spread out over 0.3Hz in frequency 
(0.27, 0.43, 0.58) and 211° in direction (56°, 101° and 136°). The polar diagram 
in Figure 74 shows the maxima and minimum positions for the 40° and 80° 
v-pol curves. No apparent symmetry of these positions with respect to any of 
the wave directions is obvious, but it is clear that away from the upwind 


direction the placement of NRCS minima is not a function of the wind alone. 


* The near equal magnitudes of the 60? data with the 40° NRCS values is 
due to the lower wind speed during the 40? portion of the data run. 


206 


зәјдиҹ aAeM pue purM ‘Jepey Buoury Ayjawoay aarjeyayy 
yo 3014 Tefoq pur urnnəəds feuonəəriq yo qq əAnəədsiəq “08 una “02 əm3iq 


COMO SON. N — 
CIPS ə ə ° 
op) XA аа, t d 
s ə — XN pO oO SOS رخ‎ 
——— o 
/ 4/42? کک کک‎ ә PEESSSSSSSSSSSSSESSSS SS SSE, 
ESS PCS OL ELSES SSS LGD PSPS SE SESS SS SSL SESS DE 
Һәбси .. — o —— IS o — — نے‎ 
Eet, eh EE 
COS A TSS SSS SPE کے‎ 
eO OS ASS SS SISOS SOS > کک‎ 
,o£r-SYY I I SELES SSS SS SSSSSESSS DSS DUDS ST SUS SSS SSSSSDSDSSS SES IS 
9 Std نے‎ ٠ ooo o rri Ssss ٠ 
XV/v/17) ooo پٹ‎ a — e eg e e e e e e 
٢ o OS a OS >L $$$ > i 2 + +X 2 ++ +; O S 
Dore, OO S D ES a əə oo ooo o سے‎ ə ii :کے‎ 
کک کک کک کک و ى‎ ə 
ə əə“ o کک کو ےو ےو کہ مد ك‎ 
с̧ор әт әәә 
: سد سو سوہ‎ oo Ән 2 $$ 25 $o کک کک‎ 
e یی کک‎ ON SOA SSS Se “ 
ə” ےک ےک ے٠ ہے دی سح سح سد ہی سر‎ S ہکوہ ےہ پر ید وش‎ e e + — edd 
ə “ox ə ə کے‎ H OS YESS OSS әс̧әсәс̧әәс̧әс̧әс̧әсә-сә-ә- 
/ یک ےہ ےک ےا‎ o o” 3 PSII SIS IOS PS Sc ےه ے7‎ 
0 ees ea e ee e 
SUMI əN D o... 
-t bar ب‫‎ 
سے رع ویک یک وک و کک یک‎ 0s کیہ ےہ‎ ə ə SS SSS SS 


w ə |‏ 
06 کیک یک یی یک یسیک یلص ا 





o “o o ns Y 
O 22: SPS PS PSO Sg M oo əc. ove, 
a. Ov I 3۸۶۸ 1J0US  gpz 
LSS IS OE OSS SSS OSS OPER LS BET ESOS i +; $ 
ooo” o ےك ے٠ ہے‎ ٣ کے ہے‎ tg Ké 
CC gt rt Eh Et EE و و و‎ 
Ly سک پر سر 7ر‎ <. 
5 EU X 
کک کک کک ی و کی یک یک و کیک‎ S Y 


2. 3” 
—— o. ek E rods, 


9 3 o o ors Z 
dE PM 
әд Z 
"ә لا‎ 23 
SS zer 
1983 wa +0 
NL sad 19 


)V222 325 әләм Du0? GOU 


U} JON 


207 


WAVES IY --- Novcinber 17 Run 87075 Time 1530 
Avgs : 1 niin blocks Overall :— ligi Wind. Short FYC Stable 
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Figure 71. One-Minute Block-Averaged C-Band NRCS and Environmental 
Data for Run 75 
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4. Dependence of the NRCS on Descriptive Spectral Parameters 

The azimuthal modulation at 10° and 20° was frequently observed to 
be 5 to 10dB which, for 10°, is significantly more than that predicted by the 
model functions. An initial visual examination of the data plots suggested a 
correlation between the amount of modulation and sea state, perhaps 
significant wave height. The modulation did not appear to be related directly 
to wind speed since both high and low wind speed cases were observed. A 
closer inspection of the data revealed that the modulation was greatest when 
the peak frequency of the wave amplitude spectrum increased to frequencies 
corresponding to wavelengths of 5m or less. This suggested a comparison of 
the modulation with wave steepness parameters, since shorter, steeper 
waves increase the mean-square slope of the surface, contributing to greater 
specular scattering at lower incidence angles. 

Two spectral wave parameters which are indicators of wave steepness are 
the significant slope, § (eqn 141), and wave age, W (143). These quantities 
were objectively computed using the peak frequency from the wave 
amplitude spectrum. It will become clear in the ensuing discussion that 
misleading calculations of quantities based on the spectral peak can occur for 
doubly peaked spectra. The significant slope and wave age are "integral" 
representations of the wave field during the measurement period; therefore, 
a similar "integral" representation of the radar measurement is required for 
proper comparison of the NRCS to these wave parameters. 

The coefficients of the cosine curves which approximate the NRCS 
azimuthal dependence (Section C) are used for comparison with the spectral 


wave parameters. In particular, the ratio of the A to Ao coefficients reflects 


the amount of cosine modulation, as normalized by the wind speed. Initial 
plots of the regression of the A2?/ Ay ratio with the wave parameters were 
scattered but much of the variance was associated with data runs which were 
taken when there was rain or low friction velocity values (< 0.12 - 0.17cm/s). 
The rainy days were immediately excluded from all incidence angles. The 
low friction velocity runs were at first eliminated only from the 40° data since 
the lower angles are less sensitive to wind speed. However, these low 
friction velocity values seemed always to be associated with outliers in the 
regression plots, so they were eventually eliminated from consideration at all 
angles?. 

The ratio of the A» to Ao coefficients are plotted against significant 
slope in Figure 75 and coded as unstable, stable or neutral according to 
negative, positive, or zero Richardson number. An increase of the ratio with 
the significant slope is observed and the unstable and stable runs are clearly 
separated at the low and high ends of the curve. This is not interpreted as a 
stability effect, rather it seems to be a combined effect of stability and fetch. A 
plot of Richardson number with direction around the lake (not shown) 


indicates that all of the stable runs occur for (warm) winds from the west 


? Based on the NRCS vs. friction velocity plots and these analyses, a 
lower limit of 0.12 - 0.17cm/s for friction velocity became a subjective editing 
criterion in order for the radar to behave "sensibly". This could indicate a 
threshold wind speed value for the radar as suggested by Donelan and Pierson 
(1987) however, there was no systematic dropoff of the NRCS values with 
wind speed observed here, only an increase in the scatter of the NRCS values 
when the low power returns extended beyond the dynamic range of the radar. 


ҺӘ 
نیم‎ 
Q3 


(short fetch), unstable runs occur mostly from the North (polar air and short 
fetch), neutral runs occur from the longer fetch directions (a marine layer). 

A cluster of points (numbered and inside the dashed line, Figure 75) 
degrades the correlation of NRCS on the significant slope. These runs can be 
eliminated as follows: Runs 44, 52, 53 had low friction velocity values as 
explained in the footnote; Run 147 was at the end of a rain event and the 
remaining outliers all had low energy, bimodal wave spectra. For these cases, 
the spectral density of the lower frequency peaks was greater than, but nearly 
equal to that of the higher frequency peak; therefore, the lower frequency 
peak was chosen for the computation of the significant slope. However, the 
total variance of the spectrum was increased by the presence of the higher 
frequency waves so that in (141) the significant slope based on the peak of the 
spectrum appears "misplaced". The removal of these anomalous points 
improved the correlation from 0.62 to 0.91 (lines 1 and 2, Figure 76). This 
correlation indicates that the azimuthal modulation increases as the 
steepness of the dominant wave increases. 

The significant slope is related to the wind in the sense that the total 
variance of the spectrum increases with the wind speed, and the wavelength 
corresponding to the peak frequency increases. However, the radar cross- 
section is only affected by the wave slope at second order: the wind speed or 
Stress 1s still the primary factor in the variation of the NRCS. Therefore, the 
wave age 1S a more appropriate spectral parameter to compare to the 
NRCS.Utilizing the same editing criteria as described previously, A2/Ao at 20° 


h-pol was plotted against inverse wave age in order to yield a positive slope 
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on the regression line (Figure 77). The correlation of A2/ Ao on inverse wave 
age is 0.91, after removal of the anomalous points listed above. Since the Ao 
coefficient is highly correlated with wind speed, a multiple regression of 
A2/Ao on both wind speed and inverse wave age was performed to eliminate 
the possibility of any spurious correlations with wind speed. The partial 
correlation coefficient, adjusted for wind speed dependence, is 0.83. 

The positive dependence of Az/Apgon inverse wave age is 
interpreted to mean that the azimuthal modulation at 20° increases for a 
growing sea composed of shorter, steeper waves. This result holds for 10° and 
20°, both polarizations; at 40°, there is still a positive, albeit reduced 
correlation of A»/ Ao on inverse wave age of 0.75 (Figure 78). For these higher 
incidence angles, the quasispecular return is reduced for increasing slopes, 
and the Bragg scattering term dominates (12). 

The two-scale model by Plant (1986) expresses the NRCS in terms of 
the upwind and crosswind slopes. An adjustable parameter, €, of the model 
represents the ratio of the crosswind to upwind slopes to be adjusted, where € 
as set by Plant is 0.4. The Plant model was run to examine the effect of 
changing slopes on the azimuthal modulation at 20°, this parameter was 
varied from 0 to 1 (extreme range of 0 to 100% crosswind slope). The results, 
summarized in Table VII, show the expected inverse relationship between 
the azimuthal modulation and the crosswind slope. However, the 
modulation never disappears as was frequently observed in . E Ontario 
NRCS data (recall Test Case 3 and see Appendix B). Since these "flat" curves 


occurred when the peak frequency of the wave spectrum was low, it was 
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hypothesized that the reduction in the modulation vvas due to a decrease in 
the total mean-square slope. The Plant model vvas run again vvith the mean 
square slope halved and & equal to 0.4. The results show that the the 
upwind/downwind (UW/DW) difference increases while the crosswind 
minima maintain about the same value. Therefore, in this model the total 
slope determines the UW/DW values and the ratio of the upwind/ crosswind 
slopes determines the amount of azimuthal modulation. The observed 
"flat" curves which prompted this sensitivity could not be reproduced with 


any realistic values of the slopes or crosswind/upwind ratio. 


TABLE VII. PREDICTED EFFECTS OF SLOPE CHANGES ON NRCS FROM 
PLANT MODEL (20°, WIND SPEED 6M/S) 
UW/CW Difference UW/DW Difference 


vpol hpol vpol hpol 
12-5 13.4dB 0.5 


10.8 05 


9.4 ۱ 0.6 

8.1 : 07 
Maximum no longer at upwind. Assumption 
of collinear wind and wave angles is violated. 
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E. DEPENDENCE OF NRCS ON ATMOSPHERIC STRATIFICATION 

The previous sections have shown the NRCS to be functionally 
dependent on wind speed, wind and wave directions and the state of wave 
development. If the atmospheric stability is now entered as an influencing 


variable, then the NRCS may be written as follows 


NRCS - NRCS(t ut 9, y, z/L, 52) (149) 


where uj is either friction velocity or wind speed; q, y are the wind and wave 
directions; z/L is the stability parameter; and 5215 (һе mean square slope. 

In the effort to define an operational algorithm for the NRCS, the wind 
speed, wind direction and stability dependencies are accounted for in the 
SASS1 and SASS2 algorithms, which parameterize the NRCS as a power law 
of neutral wind and a Fourier series in relative wind direction. Studies 
subsequent to SEASAT (Keller et al., 1985; 1989) revealed that the NRCS was 
also related to atmospheric stratification and long wave slope (recall Figure 2), 
although a unique relationship of NRCS on friction velocity, wave slope and 
stability could not be established for all stability conditions. The conclusion of 
Keller et al.(1985) was that the surface layer theory used to estimate the 
friction velocity did not properly account for stability, or that the NRCS was 
not determined by the stress in any simple fashion. 

The problem lies in the assumption of independence of the variables 
contributing to the NRCS measurement. The Lake Ontario data provide an 
extreme example of how tangled the web can become when the 
interdependence of the physical parameters is considered. It was noted earlier 


in the discussion of the NRCS azimuthal variations that the stability varied 
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systematically vvith direction around the lake: the most unstable data vvere 
from the North-Northvvest (cold polar air), near-neutral, from the North- 
Northeast (marine winds), and most stable from the South-Southwest. 
Inasmuch as the fetch gradient of the lake has a similar directionality, stability 
and fetch, hence stability and sea state, are highly correlated in this study. 
This correlation was displayed in Figure 75, in which all but one of the stable 
values of the ratio of Az to Ag corresponded to higher values of significant 
wave slope or inverse wave age. 

The comparison of the NRCS values with any particular geophysical 
parameter then requires that al] other influencing variables be held constant. 
This filtering of the data quickly reduces many thousands of points to a few 
dozen. Therefore, to keep as many of the Ku-band radar NRCS values as 
possible in this examination of stability effects, and to eliminate the influence 
on the NRCS from the primary factors of wind speed and direction, the data 
were normalized using the fitted sinusoidal curves which approximate the 
dependence of NRCS on relative wind angle as shown in Section C. The 
difference of the measured NRCS (in linear units) and the fitted value at a 
particular azimuth was computed at each data point to account for the wind 
speed dependence of each curve. All of the residuals were then normalized 
to the "upwind" value, for a particular run, by adding to it the sum of the 
three cosine coefficients (eqn(91) with x = 0); i.e., all points were translated to 
the same relative angle position. 

The normalized "upwind" points were stratified by ranges of a particular 
wind parameter and plotted against air-sea temperature difference, bulk 


Richardson number and z/L as determined from the bulk Richardson 


bə 
bə 
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number and Donelan's ratio of z/L to Richardson number. The residuals for 
20? incidence (v-pol) stratified on measured wind speed at 11.5m, Un(19.5), 
UN(1), and friction velocity are shown in Figures 79 to 82, in which the 
measured wind speed range is from 6.5 to 7.5m/s and the other wind 
parameters were scaled to correspond with this wind speed range. The NRCS 
values show a steady decline amounting to about 5dB in 25°C. The obvious 
decrease of the cross-section with increasing stability is not removed by any of 
the wind speed parameters. 

A similar decrease of the NRCS is observed in the 40° v-pol data for the 
same measured wind speed range of 6.5 to 7.5m/s (Figure 83). However, 
stratification of the NRCS residuals on the corresponding friction velocity 
(which is independent of stability) or neutral wind (which has been corrected 
for stability) removes most of the dependence of the NRCS on stability 
(Figures 84 to 86). Comparison of the plots based on the three 
stability-independent wind speed parameters used, shows the neutral wind at 
lm appears to be more effective than the friction velocity at removing the 


stability dependence. 
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Ku-band Vertical Polarization 20.0 Deg 
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Figure 79. Normalized 20°, V-Pol Ku-Band NRCS for Measured Wind Speed 
of 6.5 to 7.5m/s as a Function of a) Air-Water Temperature Difference, b) Bulk 
Richardson Number, and c) Stability Parameter. 
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Ku-band Vertical Polarization 20.0 Deg 
Neuti9m Windspeed Range: 6.00 to 7.00 
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Figure 80. Normalized 20°, V-Pol, Ku-Band NRCS for 7.0<Un(19.5)<8.0m/s as 
a Function of a) Air-Water Temperature Difference, b) Bulk Richardson 
Number and c) Stability Parameter 
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Figure 81. Normalized 20°, V-Pol NRCS for 4.5<Un(1) < 5.5 m/s as a 
Function of a) Air-Water Temperature Difference, b) Bulk Richardson 
Number and c) Stability Parameter 
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Ku-band Vertical Polarization 20.0 Deg 
Friction Velocity Range: 0.23 to 0.26 
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Figure 82. Normalized 20°, V-Pol NRCS for 0.23<u,< 0.26m/s as a Function 


of a) Air-Water Temperature Difference, b) Bulk Richardson Number and c) 
Stability Parameter 
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Ku-band Vertical Polarization 40.0 Deg 
Windspeed Range: 6.50 to 7.50 
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Figure 83. Normalized 40°, V-Pol Ku-Band NRCS for Measured Wind Speed 
of 6.5 to 7.5m/s as a Function of a) Air-Water Temperature Difference, b) Bulk 
Richardson Number, and c) Stability Parameter. 
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Ku-band Vertical Polarization 40.0 Deg 
Neuti9m Windspeed Range: 7.00 to 8.00 
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с 
Figure 84. Normalized 40°, V-Pol, Ku-Band NRCS for 7.0<UN(19.5)<8.0nV/s as 
a Function of a) Air-Water Temperature Difference, b) Bulk Richardson 
Number and c) Stability Parameter 
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Ku-band Vertical Polarization 40.0 Deg 
Neutlm Windspeed Range: 4.50 to 5.50 
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Figure 85. Normalized 40°, V-Pol NRCS for 4.5<UN(1)< 5.5 m/s as a Function 
of a) Air-Water Temperature Difference, b) Bulk Richardson Number and c) 
Stability Parameter 
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Ku-band Vertical Polarization 40.0 Deg 
Friction Velocity Range : 0.23 to 0.26 
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Figure 86. Normalized 40°, V-Pol NRCS for 0.23<u,< 0.26m/s as a Function 


of a) Air-Water Temperature Difference, b) Bulk Richardson Number and c) 
Stability Parameter 
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The simple functional expression for the NRCS stated at the beginning of 
this section (149) can be used in the interpretation of these observations. 
Given that the NRCS decreases with increasing stability for a particular 
measured wind speed mandates that a stability independent wind parameter 
must be used in any parameterization of NRCS on wind speed. This would 
account for the wind parameters and stability in (149). However, the 
observations at 20° which show the stability dependence to remain, even after 
Stratification of the data with a stability independent wind speed parameter, 
implies that the accounting for stability is not complete. Since the azimuthal 
asymmetry in the NRCS at 20° is highly correlated with long wave slope and 
the stability is similarly correlated with slope (Figure 75), it is suggested here 
that the "stability dependence" observed at 20° is actually a guise for a 
stability /slope dependence (Janssen and Komen, 1985). 

At 40°, where Bragg-resonance dominates and the slope is not as 
important, the friction velocity or neutral wind speed fully accounts for 
atmospheric stratification. The observation that the neutral wind speed at 1m 
seems to remove the dependence most completely, negates Keller et al.'s 
(1985) suggestion that the stability correction to the friction velocity for 
neutral wind speed (equation 81) is inadequate since the same correction was 
applied here. Instead the results affirm their second conclusion, which was 
that the NRCS is not a function of wind stress alone. 

These findings have bearing on the surface layer theories incorporated in 
the model functions. Evidently, the choice of a neutral wind 
parameterization at 19.5m was appropriate for the statistical models (SASSI 


and SASS2) in spite of being physically unsatisfying. Assuming a logarithmic 
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wind speed profile enables the neutral wind speed at a higher (measurable) 
level to be lowered to a near-surface level; in this analysis the neutral wind at 
1m was used with some success. This supports the contention by Donelan 
and Pierson (1987) that the wind speed nearer the surface is the controlling 
parameter, however, it does not support their assumption that the neutral 
profile can be extended to within 1cm of the surface. Since the neutral wind 
speed in this study was computed from a direct measure of the wind stress 
and the stability parameter, and by virtue of the observed slope effect, Plant's 
model which incorporates friction velocity and long wave slope could be 
exercised to examine the combined friction velocity/stability/slope effects. 
However, in its current form, neutral conditions are assumed so that 
E ration of a stability correction to the friction velocity is required before 
the effect of stability can be evaluated. 

Finally, the Durden and Vesecky (1986) model does include a stability 
correction to the neutral wind speed. Based on the observations here that the 
neutral wind is effective at removing the dependence of the NRCS on 
stability implies that the model NRCS predictions should similarly have no 
dependence on air-sea temperature difference, at a constant neutral wind 
speed. To test this, the D&V model was run for air-sea temperature 
differences from -20°C to 9°C for 40° incidence, v-pol, and Un(19.5) = 7m/s. 
The NRCS predictions show a 2.8dB decrease from -14.1dB at -20°C to -16.9dB 
at +10°C; i.e., there is still a residual stability dependence predicted by this 
model. As stated previously, the expressions used for the drag coefficient and 


Zoin this model are inconsistent. Correcting the numerical coefficients in the 


— 


Zo expression may eliminate altogether the small dependence of the NRCS on 


stability predicted by this model. 


F. NRCS DEPENDENCE ON VVATER TEMPERATURE 

The Donelan and Pierson model function specifies the waves in the Bragg 
domain to be the result of a balance between wind input and viscous 
dissipation. Since viscosity is inversely proportional to water temperature, 
the NRCS predictions from their model are water temperature dependent. 
The differences are small, e.g., only 1dB between the NRCS values for water 
temperatures of 0°C and 30°C and wind speed of 7m/s. In the WAVES87 
experiment, the water temperature variation was only about 4°C; therefore, 
the predicted effect of water temperature on the NRCS could not be verified. 
Examination of the WAVES86 data, which spanned 10°C, also showed no 


discernible dependence of NRCS on water temperature. 


G. SUMMARY 

In this chapter the dependence of the NRCS on the individual and 
combined effects of wind speed, wind direction, wave direction and spectral 
parameters, atmospheric stratification and water temperature were discussed. 
Where indicated, model function experiments were performed to provide 
additional insight into the observations. 

The NRCS was shown to have a greater dependence on neutral wind 
speed than that of the model predictions based on open ocean NRCS data. It 
was found that more favorable comparisons between the NRCS versus wind 
speed curves were obtained after re-scaling the data to account for the steeper 


waves, hence larger drag coefficient of the lake. The wind direction and 
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friction velocity dependence vvas examined in three-dimensions to confirm 
the characteristic sinusoidal modulation of the NRCS and the positive slope 
of the surface vvith increasing friction velocity. Individual curves of the vvind 
speed as a function of the relative angle vvere examined vvithin the context of 
vvave effects. Placement of the minima of these curves appeared to be related 
to the angular spread of the vvind vvaves, taken here to be indicative of the 
spreading of the Bragg e Increased upwind/crosswind modulation of 
the cross-sections at lower incidence angles was determined to be a function 
of significant wave slope. It was suggested that the effect of atmospheric 
Stratification at 20° was not completely accounted for by simply using a 
stability independent wind speed parameter since there was a correlation of 
stability and slope in this experiment. However, it was found at 40°, where 
Bragg resonance dominates, and scatterometers are designed to operate, the 
selection of friction velocity or neutral wind speed removed the stratification 
dependence with the neutral wind speed at 1m providing the best 
normalization. A parametric check of the effect of water temperature on the 
NRCS could not be performed due to the limited range of water temperatures 


encountered during the experimental period. 
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VI. SUMMARY AND CONCLUSIONS 


A. SYNOPSIS OF THE STUDY AND MAJOR RESULTS 

Microwave radar scatterometers transmit radiation of centimetric 
wavelength to the sea surface and receive the averaged power reflected from 
the sea surface by ocean waves of approximately one-half the electromagnetic 
wavelength. The number and distribution of ocean waves of the appropriate 
centimeter wavelength are primarily due to the speed and direction of the 
wind. Therefore, it is possible to relate the intensity of the power received by a 
radar scatterometer to the wind velocity at the sea surface in the form of a 


"transfer function," commonly referred to as a "model function." 

A satellite scatterometer was flown on a oceanographic satellite, SEASAT, 
in 1978, to demonstrate the concept of measuring global oceanic winds from 
“space. The data from this mission, processed according to a model function 
known as SASS1, were analyzed to produce detailed maps of the synoptic 
oceanic wind field with resolution previously unimaginable. Unfortunately, 
systematic errors in the SASS1 algorithm were discovered in subsequent | 
analyses and additional empirical studies implicated effects on the measured 
radar returns due to geophysical variables other than the wind. The 
documentation of the SASS1 weaknesses, the appearance of new model 
functions in the literature, and the suggestions of additional environmental 
effects on the model function provided the motivation for the present study. 


The development of a scatterometer model function requires 


consideration of electromagnetic, ocean wave and air-sea surface layer 
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theories. Concepts from each of these areas pertinent to scatterometry vvere 
discussed to provide the background for the description of existing model 
functions. Based on these model functions, environmental parameters were 
identified which may influence the fundamental measurement of the 
scatterometer, the Normalized Radar Cross-Section (NRCS). 

The Water-Air Vertical Exchanges 1987 experiment (WAVES'87), 
designed in part to evaluate the effects of environmental parameters on the 
NRCS, was performed from a research tower located in Lake Ontario, on 
which two radar scatterometers operating at Ku- and C-band (14.0 and 5.0GHz) 
were installed for six weeks in the autumn of 1987. A novel aspect of the 
radar portion of this experiment was that the Ku-band radar was installed on 
a rotating, incidence angle adjustable, antenna platform. Therefore, both the 
azimuth and incidence angle variations in the received power were 
examined within the context of environmental effects. The noteworthy 
aspects of the in-situ measurements were the simultaneous acquisition of 
wind stress from a bivane anemometer and high resolution directional wave 
spectra. The attributes of the scatterometers and the supporting 
meteorological/limnological instrumentation were described, as well as the 
data processing required to reduce the uncalibrated voltage measurements to 
a set of collocated radar cross-section and geophysical variables. Using this 
collocated data set, the dependence of the NRCS on each of the individual 
effects of wind speed, wind direction, wave spectra and spectral parameters, 
and atmospheric stability was examined. The limited (4°) range of water 
temperature encountered during this experiment prevented a similar 


parametric analysis of the effect of water temperature on the NRCS. 
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The NRCS measured by the multiple incidence angle, Ku-Band 
scatterometer in Lake Ontario exhibits the characteristic incidence angle 
behavior as a function of wind speed at 19.5m measured in previous studies: 
ie. a negative dependence on wind speed at nadir, transitioning through 
almost no wind speed dependence at 10? and a positive dependence at 40? (2.4, 
+0.11, upwind, vpol) to 60° (2.7, +0.20, upwind v-pol). The slopes of ‘the 
NRCS curves as a function of neutral wind speed are steeper than the 
corresponding curves obtained from existing model functions, which are 
based on open ocean conditions (cf., 2.4 and 1.7, 40° v-pol, upwind). Scaling 
the lake data to simulate open ocean conditions by using the ratio of a lake 
drag coefficient to an ocean drag coefficient reduces the slopes of the NRCS 
versus neutral wind speed curves and makes them parallel to those predicted 
by the model functions (cf., 1.6 and 1.7, 40°, v-pol). The successful adjustment 
of the curves using this lake/ocean drag coefficient scaling indicates the drag 
coefficient should include a sea state dependence. If the drag coefficient is 
considered to be a function of wind speed alone, then the nearly cubic 
dependence of the NRCS on neutral wind also implies the NRCS is 
proportional to the wind stress. 

However, the correlation coefficients for the regressions of NRCS 
against neutral wind speeds at 19.5m and 1m and friction velocity were 
higher for the neutral wind speed than for the friction velocity (cf., 0.91 and 
0.79, 40°, upwind, v-pol). This statistic alone is insufficient to declare the 
neutral wind speed to be a better predictor of NRCS. Although the neutral 
wind speeds, at levels higher and lower than the anemometer height, were 


deduced from the friction velocity in this experiment, they were smoothed by 
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the ratio of the measured and neutral drag coefficients in the conversion 
process. Therefore, the neutral wind speeds have less inherent variability 
and the correlations are improved. 

The two primary characteristics of the NRCS, that is, the azimuthal 
modulation of the NRCS as a function of both relative wind angle and 
friction velocity were viewed in three dimensions for the six incidence angle 
settings of the Ku-band radar. The sinusoidal modulation becomes evident at 
20° and more obvious at 40°. The 60° and 80° surfaces at v-pol continue to 
display the sinusoidal modulation whereas the sinusoidal appearance of the 
60° h-pol NRCS surface is much less apparent, and not at all apparent at 80° 
where the NRCS measurement is generally at the noise level of the radar. 
| The minima in the azimuthal modulation of the NRCS as a function 
of relative wind angle are often shifted by as much as 45? such that the 
minima do not always occur at 90? and 270? cross-wind angles. The addition 
of a second order sine term to the standard three term cosine fit of the 
azimuthal modulation often places the minima better. It cannot be said that 
the fits with the additional sine term improve the fits in a general sense since 
the NRCS values were not distributed evenly over 360°. The effect of rain on 
the NRCS azimuthal modulation is to increase initially, then decrease, the 
level of the Ku-band NRCS, and to flatten out the azimuthal variability in the 
NRCS; i.e., the NRCS is isotropic in the presence of rain. The C-band NRCS 
does not appear to be similarly affected by the rain, indicative of the greater 
influence of rain on 1cm Bragg waves compared with 3cm Bragg waves. 
Further comparisons of rain effects on the Ku- and C-band NRCS values will 


be performed with additional test cases identified in the WAVES'87 data. 
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The azimuthal modulation as a function of relative wave angles indicate 
that the shape of the azimuthal dependence conforms best to the sinusoidal 
model when the wave state is dominated by wind waves in the direction of 
the wind. However, the angular spread between the NRCS minima is 
generally rather less than 180°; therefore, the minima in the NRCS do not 
occur at 90° angles to the wind. In the case of a long swell component 
opposing the wind sea, the minima become symmetric to the relative angles 
between the primary wave components and the separation of the minima 
widens to 180°. The azimuthal modulation disappeared altogether in the case 
of a "confused" sea consisting of multiple peaks and broad directionality. 

The azimuthal modulation as a function of wave spectrum parameters 
is highly correlated with significant slope and inverse wave age at 20° 
(correlation coefficient of approximately 0.90) and moderately correlated with 
these same parameters at 40? (0.75). Therefore, increased azimuthal 
modulation at 20° is associated with a steeper wave field, which increases the 
specular contribution to the NRCS at this angle. At higher incidence angles, 
Bragg-resonance with the capillary waves is the dominant scattering 
mechanism. 

The dependence of the NRCS on atmospheric stability shows the NRCS 
to decrease by about 5dB between air-sea temperature differences of about -16° 
to +10°C. This means that the NRCS in an unstable environment indicates a 
wind speed at least 2m/s higher than for a stable environment, if a square-law 
dependence between the NRCS and wind speed is assumed. Since the design 
error for the scatterometer wind speed is 2m/s, improper accounting for this 


effect alone will prohibit the design criterion from being met. It was not 


240 


possible to remove the stability dependence at 20? for stratification of the 
NRCS on any wind speed parameter or friction velocity. The conclusion is 
that the assumption of independence among the variables influencing the 
NRCS was violated in this experiment since the stability and fetch (hence 
wave state) are correlated. Therefore, at 20° the failure of a wind speed 
parameter to remove the stability dependence is an indication of an 
additional correlation between siability and wave slope. At 40°, where the 
NRCS is more dependent on wind speed than mean-square slope, stratifying 
the data on either the friction velocity or neutral wind speed removes the 
Stability dependence. Therefore, model function parameterization of the 
NRCS in terms of the friction velocity or neutral wind speed accounts for the 
observed decrease in NRCS with atmospheric stability. 

In light of the interdependence of the factors influencing the NRCS, a 
diagnostic model for the NRCS must either use a physical approach to model 
the evolution of the sea surface through the wavenumber spectrum, or be 
derived from an enormous Statistical data base which represents all possible 
environmental conditions. Quantification of the relative importance of a 
particular environmental parameter can only be established as a function of 
scale and application; i.e., increased resolution and accuracy requirements of a 
Satellite sensor demand physical understanding and measurements at 
correspondingly finer scales. 

In this study, existing model functions of both types (physical and 
Statistical) have been examined to elucidate how and why certain 
environmental parameters are included in the models. The attributes and 


basic predictive characteristics of each of the models (SASS1, SASS2, Durden 
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and Vesecky, Plant, and Donelan and Pierson) were described. All of the 
model functions, with the exception of SASS1, are prototypes which have not 
been evaluated for operational use. Therefore, this analysis focussed on the 
validity of the basic assumptions regarding environmental parameters within 
the models, rather than establishing the skill of the different models. Further, 
tħe electromagnetic scattering model per se was not critically analyzed; 
however, the greater importance of the significant slope at lower incidence 
angles as compared to the intermediate incidence angles does suggest a 
composite surface approach. 

The empirical results described above indicate that two of the major 
aspects of the model functions which need to be reconsidered are the 
methods by which the wind waves are spread about the wind direction and 
the assumption of fully developed wind seas of one primary wave direction. 
Since the minima in the azimuthal modulation of the NRCS appear to 
approximate the angular spread of the wave spectrum rather than the actual 
crosswind directions, the manner in which the wave spectrum models spread 
the Bragg waves about the wind direction is important. The D&V and D&P 
models do allow for this effect by factorizing the wavenumber spectrum into 
a product of wavenumber amplitude and directional dependence terms. 
However, the adjustable parameters in these models are still tuned to 
produce minima near the cross-wind directions. In the model function by 
Plant, the directionality of the NRCS is obtained by expanding the cosine 
dependence of the wave growth expression and by the mean-square slopes of 
the longer waves. This model also produces cross-wind minima at 90° and 


270° and rather flat-topped, broad maxima. None of the models described 
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currently consider additional wave components such as long waves 
propagating at an angle to the wind. Both of these wave effects shift by 15° to 
20° the positions of the NRCS minima away from a simple +90° relationship 
with the wind; therefore, forcing the minima to be at 90° and 270° through 
the use of a three-term cosine fit will not produce an accurate representation 
of the directionality of the scatterometer measurements. 

The azimuthal modulation was also shown to be correlated with the 
significant slope of the wave spectrum. In the physical models of Plant (1986) 
and Durden and Vesecky (1985), the Cox and Munk (1954) mean-square slope 
expressions are used; however, there have been recent studies of the short 
wave slope which suggest different parameterizations. The Donelan and 
Pierson (1987) model integrates the long-wave portion of their wavenumber 
spectrum to obtain mean-square slopes, but they could not reconcile their 
numerical values with those of the Cox and Munk values. Therefore, a re- 
evaluation of the mean-square slopes used in the model functions is 
warranted. 

The third area of the models which needs updating is in the basic surface 
layer theory that connects the measured wind speed to the surface. The 
apparent effect of atmospheric stability at intermediate Bragg angles is 
properly accounted for if a stability independent wind speed parameter; i.e., 
either friction velocity or neutral wind at a level near the surface, is used in 
the parameterization of the NRCS on wind speed. If the friction velocity is 
measured at any substantial distance above the surface, model functions 
which then parameterize the NRCS in terms of a neutral wind speed at a 


lower level, need air-sea temperature differences to compute the stability 
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correction. Except for the D&V model, the models presently assume 
averaged neutral conditions over the sea, therefore, the air-sea temperature is 
not required. In either case of neutral wind speed or friction velocity, 
improvement of the model functions will depend on improved 
measurements of the friction velocity at sea, either for direct comparisons to 
the NRCS or for the establishment of an open ocean drag coefficient which is 


normalized for wind speed, stability and sea State. 


B. CLOSURE 

This study began by establishing the perimeter of the scatterometry model 
function problem. Contained within this perimeter were the scientific issues 
facing model function developers, the special operational requirements of the 
Navy for surface winds and the programmatic status of scatterometry in the 
US and abroad. 

This study has shown that radar scatterometers are an especially sensitive 
means by which to study the air-sea interface. The magnitudes of the Bragg- 
resonant NRCS measurements from radar scatterometers respond nearly 
instantaneously to changes in the near-surface winds, but the azimuthal 
modulation of the (Ku-band) NRCS appears to be a strong function of the 
angular spread and slopes of the dominant waves on the sea surface, rather 
than a simple sinusoidal dependence on the wind direction. In this analysis, 
the neutral wind speed at 1m, obtained from the measured friction velocity 
and corrected for atmospheric stability, was a better predictor of NRCS 
magnitude and was a better normalization quantity for removing the 
observed stability dependence of the NRCS. Taken together, these 


observations state that the radar scatterometer can measure surface winds to 
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the accuracy desired by the Navy, if the interdependence of the influencing 
variables is properly taken into account. That is, by modeling the NRCS in 
terms of 1) a neutral wind speed near the surface determined from a neutral 
drag coefficient formulated in terms of sea state, friction velocity and stability, 
2) updated mean-square slopes and spreading functions of the wind waves 
and 3) the background wave field which includes wave components traveling 
obliquely to the wind vector. 

The launch of the European Space Agency satellite in the near future will 
carry a combination of instruments designed to study wind and waves and 
thus, will provide a means by which to validate the above conclusions and 


examine model function improvements. 
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APPENDIX A. DEFINITION OF THE NORMALIZED RADAR 
CROSS SECTION 


When microwave radiation strikes the sea surface it is scattered in many 
directions. Some part of this scattered radiation is eventually received by the 
radar. The mathematical relationship which describes that portion of the 
transmitted power received at the radar after having been reflected from a 
surface scatterer is called the radar equation. The definition of the radar 
cross-section arises in the derivation of the radar equation and can be thought 
of as the equivalent surface area of a target as observed at microwave 
frequencies (Long, 1989). The Normalized Radar Cross-Section (NRCS) is the 
radar cross-section divided by the target area. 

Suppose that the surface consists of a dominant scatterer, then the 
geometry and quantities related to the scattering problem are as in Figure Al. 
Following Ulaby, et al., (Chapter 7, v. 2, 1982), the final received power at the 
radar may be thought of as the result of a sequence of events starting with a 
radar transmission of power, Py, in the direction of the scatterer. The radar 
antenna pattern is nonisotropic so that the transmitted power must be 
multiplied by the antenna gain, GT, in a particular direction. Also, the 
transmitted power is spread over a sphere surrounding the antenna of radius, 


kr, equal to the distance between the antenna and the scatterer. 
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Therefore, at the scatterer, 


Power per Unit Area Incident = Power Transmitted * Spreading Loss 





on Scatterer toward Scatterer 
5 = F.Ə, 
° 4xR” (A1) 


Gain G, FSG, 
Power/Unit Solid Effective Area A, 
Angle P, G, 


Effective m 
Area A. Power Through 


Spreading Loss Lm S. 
] / nR Spreading Loss 
1/ 4۲۲ 
d BN Solid 
Power Through Angle, Pis Gis 
1 m” s S, Gain бу 
Fraction Absorbed, f, 





Figure A1. Geometry of and quantities involved in the radar equation (from 
Ulaby, et al., 1981). 


If the scatterer is considered to be a receiving antenna, then it removes 

power from the effective area of the incidence beam which it intercepts. The 
d 

effective area is that area of the incident beam from which all pogver would be 


removed if one assumed that the power going through all the rest of the 
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beam continued uninterrupted. The total povver intercepted at the scatterer is 
then 


Power Intercepted by Scatterer = Incident Power * Effective Area 
per Unit Area of the Scatterer 


Pas = SsAgs (A2) 


The actual value of Ars depends upon the effectiveness of the scatterer as a 
receiving antenna. 

The power received by the scatterer is not all reflected: some is absorbed, 
fa , the remainder 1 - fa, is reflected in various directions. That is, 


Power Reflected by Scatter = Power Intercepted by Scatterer * Fraction of 
Power not absorbed 


Prs = Pes- fa) (A3) 

Inasmuch as the scatterer receives radiation as an antenna, it also radiates 

as one. Therefore, the power leaving the surface, Prs, must be multiplied by 

the gain of the receiver, Gr. Again, the power spreads out along a spherical 
wavefront and is thus reduced. Then 


Power per Unit Area = Power radiated by Scatterer * Gain toward * Spreading 
at Receiver receiver Loss 


— Pr G; 


S, = A4 
T Ap ә 


Finally, the total power at the receiver is 


Total Power at the Receiver = Power per Unit Area at Receiver *EffectiveArea 
of Receiver 


Р, = 5,А, (A5) 
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Substituting for Sp from the previous four equations gives 


n F Grâr Ñ 
P, f Ary RR} Jac fa) Grs] an 


The quantity in square brackets describes the radiating characteristics of 
the surface scatterer. Individually, the variables are difficult to measure and 
they are not particularly relevant if the measurement of interest is the 
magnitude of the received power. Taken together, the variables define the 


radar scattering cross-section, 


O = Ass (] — f.)Grs (A7) 


The cross-section is a function of the incident and reflected wave directions, 
the shape of the surface scatterer and the scatterer's dielectric properties. 

In the case of a monostatic radar in which a single antenna is used for 
transmitting and receiving, the following equalities hold: Rr = Кр = R, Gr = 
Gr =G, Af = Ar = A. Further, the effective area of the antenna is related to 
the gain of the antenna by A=12?G/4r. The radar equation appropriate to 
scattering from a single point target using monostatic radar may then be 


written as 


"+0 


n (4 m)° R° (A8) 


However, the sea surface does not represent a single point target, rather the 


scatterers are distributed over some area À on the surface. 
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In this case, it is typical to think of the average received power rather than 
the power at a particular instant, and the scattering cross-section per unit area, 


o°, such that 


o”dA (A9) 





fe arena 


Area 


This is the area extensive form of the radar equation used for remote sensing 
radars. There are restrictions on the application of this form of the equation 


which guarantee that no one scatterer dominates: 


(1) In a differential area, AA, in which P, G, and R are nearly constant, there 
must be enough scattering centers for a reasonable average to be 
computed and allow substitution of GAA for Gj, and 


(2) In the total illuminated area there must exist many more scatterers tO 
allow replacement of the differential form with an integral. 


In a real system, to solve for o? from (A8) then requires inversion of the 
integral. Since the antenna gain, cross-section and range are not constant 
over the target area, Bracalente, et al., (1980) approximated 6? by including 
correction factors in the SASS algorithm to allow the equation for o? to be 


written as 


с DÄI PR 
CT === e 
PAL COI YUA 
T 5 IS 


o 


(A10) 


where Ls is an electronics loss term, G is the peak antenna gain, G/Go is the 
gain loss ratio (both as functions of the antenna look direction) and A is the 


illuminated area on the earth. Values from this approximating equation 
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Over a parametric range of conditions were compared with the exact equation 
to show that the error in o? due to the use of (A10) was less than 0.1dB. 

A particular o? can be evaluated from (A10) given measured values of Pg, 
PT, calculated values of R and A, and laboratory determined values of Ls, G 
and Go. Noise contamination is eliminated by sampling the noise-only 
power and subtracting it from the received power (signal) measurement. The 
problem is then how to relate the measured NRCS to the physical 


characteristics of the scattering surface. 
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